3.4 Comparing gas phase and SMD solution consistently

These are unexecuted teaching inputs and starting models. Original diagrams are schematics, not calculated results. Validate version-specific syntax, licensed or authorized data, numerical convergence and the scientific model before using this workflow.

3.4.1 Model, units and provenance

Geometry in Å; electronic energy in hartree; vibrational wavenumbers in cm⁻¹. Check each printed field and keep thermal and standard-state terms distinct.

Shared inputs, conventions and evidence

Original schematic: Comparing gas phase and SMD solution consistently. No numerical results are claimed.

3.4.2 Unexecuted inputs and explicit deltas

Use the accompanying instructions to identify the parent calculation and placement of every delta; a snippet is not automatically a standalone input. Preserve all blank-line and file-provenance requirements.

3.4.2.1 Input block 1

%oldchk=case10_water_min.chk
%chk=case12_water_smd_vertical.chk
%mem=4GB
%nprocshared=2
#p B3LYP/6-31G(d) SP SCRF=(SMD,Solvent=Water)
   Geom=AllCheck Guess=Read SCF=Tight Integral=UltraFine

3.4.2.2 Input block 2

%oldchk=case10_water_min.chk
%chk=case12_water_smd_min.chk
%mem=4GB
%nprocshared=2
#p B3LYP/6-31G(d) Opt=(Tight,CalcFC) Freq
   SCRF=(SMD,Solvent=Water) Geom=AllCheck Guess=Read
   SCF=Tight Integral=UltraFine Temperature=298.15 Pressure=1.0

3.4.3 Worked investigation

3.4.3.1 Question and intuition

What changes when a polar molecule is placed in a polarizable environment? An implicit solvent surrounds the solute with a continuum rather than thousands of explicit solvent molecules. SMD includes an electrostatic response and empirically parameterized non-electrostatic contributions. It is useful for controlled comparisons, but it does not show a particular hydrogen-bonded water shell, proton transfer through a solvent chain, or solvent exchange dynamics. We use the Case 10 water monomer as a deliberately small solute to learn the workflow, not as a microscopic simulation of bulk liquid water.

Prerequisite: retain the validated case10_water_min.chk and its gas-phase output. Two different questions are separated. A fixed-geometry comparison asks how the model energy changes without nuclear rearrangement. An independently optimized solution calculation permits structural relaxation. These differences should not be interchanged, and neither automatically includes a complete concentration-dependent experimental free energy.

3.4.3.2 Complete checkpoint inputs

First run a fixed-geometry SMD single point at the gas-phase minimum.

See input block 1 above.

Then independently optimize and compute frequencies in the same SMD environment, starting again from the preserved gas-phase checkpoint.

See input block 2 above.

Geom=AllCheck means neither file has a title line or molecule specification. Guess=Read provides a starting orbital guess; the solution SCF still converges self-consistently in the requested solvent. These are distinct output checkpoints, so the gas minimum and fixed-geometry result are not overwritten. Use the predefined Water solvent rather than manually re-entering a dielectric constant and assuming all SMD parameters remain equivalent.

3.4.3.3 Workflow

Copy the gas-phase method, basis, grid, SCF settings, charge, and multiplicity into a comparison manifest. The solvent term should be the intended change in the fixed-geometry test.
Run the fixed-geometry job, verify that the output names SMD and water, and compare its nuclear coordinates with the gas minimum. A change in molecular orientation is not itself a structural change; compare distances and angles if coordinates are reoriented.
Record the final total model energy. For the ordinary ground-state SMD calculation here, inspect the output's explicit note about non-electrostatic CDS terms being included in the reported total. Do not add the printed CDS contribution again. Keep its kcal mol−1 diagnostic unit separate from the hartree total.
Run solution Opt Freq and apply Case 10's minimum checks. Compare O–H lengths, H–O–H angle, and dipole moment with the gas calculation. State that trends are computed model responses; no particular numerical shift is promised beforehand.
Report separately ΔEfixed = ESMD(Rgas) − Egas(Rgas) and ΔErelaxed = ESMD(Rsolv) − Egas(Rgas). If constructing a solution Gibbs-energy estimate, specify exactly which thermal correction and standard-state conversion are used, following Case 11.

3.4.3.4 Checks pitfalls and interpretation

The SMD total is a model quantity that incorporates solvation contributions; calling its difference a bare electronic energy difference can conceal that meaning. Equally, calling the output's RRHO Gibbs sum a fully rigorous liquid-phase free energy overstates the treatment of translation, rotation, concentration, and solvent structure. Choose one explicit thermodynamic cycle. For example, a hybrid cycle may combine gas-phase RRHO corrections with an SMD solvation contribution, but its reference states and geometry convention must be stated and solvation must not be added twice.

Adding SCRF only to the frequency step after a gas optimization does not verify a solution minimum. Changing to a larger basis only for the solution also confounds basis and solvent effects. Charged species often require particular care with diffuse functions, cavity definitions, and the experimental convention for single-ion solvation; this neutral demonstration does not settle those questions. Explicit solvent molecules may be needed when specific coordination or hydrogen-bond networks determine the chemistry.

EN exercise: Repeat only the fixed-geometry calculation with Solvent=Ethanol, preserving all other settings. Explain which comparison isolates the solvent-model change, and which further jobs would be required to compare fully relaxed minima and consistently defined Gibbs energies. Do not infer that dielectric constant alone must determine the ordering.