DynamicsOutputs

DynamicsOutputs

Defines a set of functions that can be used to calculate outputs for dynamics simulations.

Outputs the desorption angle in degrees (relative to the surface normal) if a desorption event is detected.

(output::OutputDesorptionAngle)(sol, i)

Outputs the desorption angle in degrees (relative to the surface normal) if a desorption event was detected.

OutputDesorptionAngle(indices; surface_normal = [0,0,1], surface_distance_threshold = 5.0u"Å")

Outputs the desorption angle in degrees (relative to the surface normal) if a desorption event is detected. Use surface_normal to define the direction "away" from the surface. Most commonly, this would be in positive z direction.

A desorption is detected if the centre of mass of the molecule defined with indices is above surface_distance_threshold from the closest surface atom. This is calculated with respect to surface_normal and will take into account periodic boundary conditions.

(output::OutputDesorptionTrajectory)(sol, i)

Only output parts of the trajectory where desorption is occurring.

`OutputDesorptionTrajectory(indices; surface_normal = [0,0,1], surface_distance_threshold = 5.0u"Å", extra_frames = 0)`

Like OutputDynamicsVariables, but only saves parts of the trajectory where desorption is occurring.

Use surface_normal to define the direction "away" from the surface. Most commonly, this would be in positive z direction.

Use extra_frames to save additional steps before the desorption event begins.

A desorption is detected if the centre of mass of the molecule defined with indices is above surface_distance_threshold from the closest surface atom. This is calculated with respect to surface_normal and will take into account periodic boundary conditions.

Output a 1 if the molecule has dissociated, 0 otherwise.

Output the vibrational and rotational quantum numbers of the final image.

Output a ComponentVector with fields reflection and transmission containing the probability of the outcome. Each index in the arrays refers to the adiabatic state.

OutputSubsetKineticEnergy(sol, i)

Evaluate the classical kinetic energy of a subset of the entire system at each save step.

The subset is defined by OutputSubsetKineticEnergy(indices).

OutputAdiabaticPopulation(sol, i)

Output the adiabatic population at each timestep during the trajectory.

OutputCentroidPosition(sol, i)

Output the position of the ring polymer centroid at each timestep during the trajectory.

OutputCentroidVelocity(sol, i)

Output the velocity of the ring polymer centroid at each timestep during the trajectory.

OutputDiabaticPopulation(sol, i)

Output the diabatic population at each timestep during the trajectory.

OutputDiscreteState(sol, i)

Output the discrete state variable at each timestep during the trajectory. This is used for surface hopping simulations and returns the variable that determines the currently occupied adiabatic state.

Requires that the dynamics variable has a field state.

Use OutputDiabaticPopulation or OutputAdiabaticPopulation to get the population estimators.

OutputDynamicsVariables(sol, i)

Output all of the dynamics variables at each timestep during the trajectory.

Output the end point of each trajectory.

Output the first point of each trajectory in DynamicsVariables format. (Useful when using distributions for initial conditions.)

OutputKineticEnergy(sol, i)

Evaluate the classical kinetic energy at each timestep during the trajectory.

OutputMappingMomentum(sol, i)

Output the momentum mapping variable at each timestep during the trajectory.

OutputMappingPosition(sol, i)

Output the position mapping variables at each timestep during the trajectory.

OutputPosition(sol, i)

Output the position at each timestep during the trajectory.

OutputPotentialEnergy(sol, i)

Output the adiabatic potential energy at each timestep during the trajectory.

OutputQuantumSubsystem(sol, i)

Output the quantum subsystem at each timestep during the trajectory. Usually this will refer to a wavefunction or density matrix but will depend on the particular dynamics method.

Output the total number of surface hops during the trajectory

OutputTotalAdiabaticPopulation(sol, i)

Output the total adiabatic population at each timestep during the trajectory.

OutputTotalDiabaticPopulation(sol, i)

Output the total diabatic population at eah timestep during the trajectory.

OutputTotalEnergy(sol, i)

Evaluate the classical Hamiltonian at each timestep during the trajectory.

OutputVelocity(sol, i)

Output the velocity at each timestep during the trajectory.