Skip to main content
Sandbox Physics
Nonadiabatic molecular dynamics

L4 · Chapter 9

Reading interference and molecular-benchmark papers

Sandbox Physics ·

Open Curchod, Agostini and Gross, arXiv:1604.02351v2. Identify its model and preparation before reading Figure 4, then follow its potential representation. Keep the 2025 molecular-benchmark roadmap v3 nearby: one source examines interference, the other asks how comparable calculations should be organized.

Begin with Figure 3 and equation (15)

The 2016 paper uses a Shin–Metiu system with two fixed ions, a moving ion and an electron on one line. Equation (15) retains electron and nuclear coordinates. Write down the moving coordinates, masses and units before inspecting the energy curves.

The browser’s Tully matrix is a different model. Similar avoided-crossing shapes do not establish identical dynamics. Reproducing a source result requires its Hamiltonian, electronic basis and preparation.

Reading: Curchod

Read the preparations in equations (16)–(18)

The paper compares coherent preparation with an independent-wavefunction approximation, or IWA. Distinguish evolving an amplitude sum from adding probabilities after separate propagation; the latter removes interference between histories reaching a common channel.

For another paper, record the nuclear amplitude, relative phase, width and momentum on each electronic state. If preparation is mixed, record the component weights. An average energy alone rarely reconstructs the state.

Reading: Curchod

Read Figure 4 from its axes

The left panel scans initial momentum and compares final lower-state population from exact quantum dynamics and trajectory surface hopping (TSH). The right panel compares exact time evolution with the independent-wavefunction approximation (IWA) for two selected momenta. These panels test different approximations.

In the right panel, hide the legend and predict which calculation retains the relative phase between histories. Then check your prediction. Ask whether the same approximation would change if the two packets never reached the coupling region together. When reading an oscillation, distinguish coherent structure from sampling noise and numerical error.

Reading: Curchod

Exact factorization reorganizes the same wavefunction

Equations (5)–(6) factor the full state into nuclear amplitude χ(R,t)\chi(R,t) and normalized conditional electronic state ΦR(r,t)\Phi_R(r,t). Its nuclear density equals the electronic integral of the full density.

Verify this relation independently in an orthonormal channel expansion: electron integration removes cross terms and leaves a sum of nuclear component densities. Exactness here concerns an information-preserving representation. Approximate propagation equations for its factors still need validation.

Figures 5 and 7 examine the associated time-dependent surface; Figure 6 tracks an overlap indicator. First identify regions carrying appreciable nuclear probability and the time of recombination. Low-density regions can be numerically sensitive. Relate a feature to the probability passing through it before assigning dynamical importance.

Ψ(r,R,t)=χ(R,t)ΦR(r,t),∫∣ΦR(r,t)∣2 dr=1\Psi(r,R,t)=\chi(R,t)\Phi_R(r,t),\qquad \int |\Phi_R(r,t)|^2\,dr=1
The conditional electronic state depends on both nuclear position and time.
∣χ(R,t)∣2=∫∣Ψ(r,R,t)∣2 dr=∑n∣χn(R,t)∣2|\chi(R,t)|^2=\int|\Psi(r,R,t)|^2\,dr=\sum_n|\chi_n(R,t)|^2
The product and orthonormal-channel forms describe the same nuclear density.

Reading: Curchod

A geometric-phase result that can be checked completely

Read Berry 1984 equation (6) and section 3, especially equation (19). A circuit enclosing the real two-state degeneracy gives returning factor −1-1; a nonenclosing one gives +1+1. Derive the centered-loop sign with chapter six’s half-angle eigenvector before inspecting the overlap product.

This limited reproduction checks electronic loop phase and attaches it to a prescribed two-arm output. It supplies an analytically testable research step. Predicting a molecular product yield additionally requires that molecule’s nuclear paths, weights and transitions.

Reading: M. V. Berry

Turn the 2025 roadmap into a calculation plan

Fix arXiv:2502.14569v3. Read sections II C, II D and II E in order: electronic-structure input, matched preparation, then selected properties and observables. Section II F addresses comparison with experiment.

Write four questions in your own words. Are you comparing dynamics or electronic-structure input? Do different representations prepare the same physical state? Which outputs describe electrons and which describe nuclei? Could numerical and sampling errors explain the discrepancy? Each answer should point to a saved input or measurement.

Open research still concerns representative molecular systems, matched preparations across method families and tractable observables. Low-dimensional problems isolate mechanisms; higher-dimensional molecular calculations test how strongly those mechanisms affect an application. Connecting the two develops transferable evidence.

Reading: Cigrang et al.

Try the prediction

Supply what a paper reader needs to repeat a comparison

  1. Export one preset and state its model, preparation, exits and cost definition in one sentence each.
  2. Identify grid-reference and trajectory-statistical results and their distinct error sources.
  3. When comparing source figure points, also record figure number, momentum, width conversion and digitization error.

Check the Berry loop sign step by step

  1. Predict the returning sign from the lower-state eigenvector around one centered circuit.
  2. Compare enclosing and nonenclosing paths and inspect the same port center.
  3. Separate the electronic phase result from the ideal recombiner assumptions.

Exercises and explanations

  1. How can a single product in exact factorization retain electron–nuclear correlation?

    Read the explanation

    The electronic factor depends on nuclear position and time. It is not a fixed electronic state independent of the nuclei; local normalization and coordinate dependence retain the full wavefunction’s information.

  2. What should follow an improved result for one avoided crossing?

    Read the explanation

    First verify matched preparation, observables and convergence. Then choose a distinct physical challenge, such as return passage, coherent recombination or multidimensional intersection, and report changes for each observable.

Sources and calculation scope

  1. M. V. Berry, Quantal phase factors accompanying adiabatic changes, Proc. R. Soc. Lond. A 392, 45–57 (1984)

    Published primary derivation: equations (5)–(6) separate dynamical and geometrical phase; section 3 treats a degeneracy. The path experiment is a controlled two-state model, with the reduction and phase comparison specified on its page.

  2. Curchod, Agostini & Gross, An Exact Factorization Perspective on Quantum Interferences in Nonadiabatic Dynamics, arXiv:1604.02351v2

    Fixed 30 June 2016 author version; published in J. Chem. Phys. 145, 034103 (2016), DOI 10.1063/1.4958637. Reading scope: equations (5)–(6), (14)–(18), and Figures 3–7. Its Shin–Metiu electron–ion calculation differs from the Tully and conical teaching models.

  3. Cigrang et al., Roadmap for Molecular Benchmarks in Nonadiabatic Dynamics, arXiv:2502.14569v3

    Fixed 4 July 2025 author version associated with the published Perspective, DOI 10.1021/acs.jpca.5c02171. Sections II A–F discuss systems, dynamics methods, electronic structure, initial conditions, observables and experimental comparisons. This is a benchmark roadmap, not a completed universal ranking.

FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy