Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Molecules in high pressure environments #467

Draft
wants to merge 8 commits into
base: master
Choose a base branch
from

Conversation

hugoastrom
Copy link

@hugoastrom hugoastrom commented Oct 11, 2023

  • Confinement potential class
  • Incorporate cavity function
  • Test for atoms
  • there seems to be problems with convergence


double ConfinementFunction::evalf(const mrcpp::Coord<3> &r) const {
double f = 0.0;
const mrcpp::Coord<3> center = {0.0, 0.0, 0.0};
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is the center hard-wired to the origin? Maybe it would be wise to be able to pass this from input.

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The plan is to obtain the centers as it is implemented in the Cavity function but it is under construction

Copy link
Contributor

@ilfreddy ilfreddy left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I left a few comments for you to consider.

@@ -268,6 +269,7 @@ json driver::scf::run(const json &json_scf, Molecule &mol) {
auto relativity = json_scf["scf_solver"]["relativity"];
auto environment = json_scf["scf_solver"]["environment"];
auto external_field = json_scf["scf_solver"]["external_field"];
auto confinement_potential = json_scf["scf_solver"]["confinement_potential"];
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Seems you have some indentation issues which should be fixed

src/driver.cpp Outdated
/////////////// Confinement Potential /////////////////
///////////////////////////////////////////////////////////
if (json_fock.contains("confinement_operator")) {
auto r_0 = json_fock["confinement_operator"]["r_0"];
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

indentation should be 4 spaces :-)

@Gabrielgerez
Copy link
Contributor

About the indentation, I think you can set git so that it can automatically format your files before committing changes

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants