BioSimSpace.Process.Amber

class BioSimSpace.Process.Amber(system, protocol, exe=None, name='amber', work_dir=None, seed=None, property_map={})

A class for running simulations using AMBER.

__init__(system, protocol, exe=None, name='amber', work_dir=None, seed=None, property_map={})

Constructor.

Parameters:
  • system (System) – The molecular system.
  • protocol (Protocol) – The protocol for the AMBER process.
  • exe (str) – The full path to the AMBER executable.
  • name (str) – The name of the process.
  • work_dir – The working directory for the process.
  • seed (int) – A random number seed.
  • property_map (dict) – A dictionary that maps system “properties” to their user defined values. This allows the user to refer to properties with their own naming scheme, e.g. { “charge” : “my-charge” }

Methods

__init__(system, protocol[, exe, name, …]) Constructor.
addArgs(args) Append additional command-line arguments.
addToConfig(config) Add a string to the configuration list.
clearArgs() Clear all of the command-line arguments.
command() Return the command-line string used to run the process.
deleteArg(arg) Delete an argument from the dictionary.
exe() Return the executable.
getAngleEnergy([time_series, block]) Get the angle energy.
getArgString() Get the command-line arguments string.
getArgStringList() Convert the argument dictionary into a list of strings.
getArgs() Get the dictionary of command-line arguments.
getBondEnergy([time_series, block]) Get the bond energy.
getCentreOfMassKineticEnergy([time_series, …]) Get the kinetic energy of the centre of mass in translation.
getConfig() Get the list of configuration file strings.
getCurrentAngleEnergy([time_series]) Get the current angle energy.
getCurrentBondEnergy([time_series]) Get the current bond energy.
getCurrentCentreOfMassKineticEnergy([…]) Get the current kinetic energy of the centre of mass in translation.
getCurrentDensity([time_series]) Get the current density.
getCurrentDihedralEnergy([time_series]) Get the current dihedral energy.
getCurrentElectrostaticEnergy([time_series]) Get the current dihedral energy.
getCurrentElectrostaticEnergy14([time_series]) Get the current electrostatic energy between atoms 1 and 4.
getCurrentHydrogenBondEnergy([time_series]) Get the current hydrogen bond energy.
getCurrentKineticEnergy([time_series]) Get the current kinetic energy.
getCurrentNonBondedEnergy14([time_series]) Get the current non-bonded energy between atoms 1 and 4.
getCurrentPotentialEnergy([time_series]) Get the current potential energy.
getCurrentPressure([time_series]) Get the current pressure.
getCurrentRecord(record[, time_series, unit]) Get a current record from the stdout dictionary.
getCurrentRecords() Return the current dictionary of stdout time-series records.
getCurrentRestraintEnergy([time_series]) Get the current restraint energy.
getCurrentStep([time_series]) Get the current number of integration steps.
getCurrentSystem() Get the latest molecular system.
getCurrentTemperature([time_series]) Get the current temperature.
getCurrentTime([time_series]) Get the current time (in nanoseconds).
getCurrentTotalEnergy([time_series]) Get the current total energy.
getCurrentVanDerWaalsEnergy([time_series]) Get the current Van der Vaals energy.
getCurrentVirial([time_series]) Get the current virial.
getCurrentVolume([time_series]) Get the current volume.
getDensity([time_series, block]) Get the density.
getDihedralEnergy([time_series, block]) Get the dihedral energy.
getElectrostaticEnergy([time_series, block]) Get the electrostatic energy.
getElectrostaticEnergy14([time_series, block]) Get the electrostatic energy between atoms 1 and 4.
getHydrogenBondEnergy([time_series, block]) Get the hydrogen bond energy.
getKineticEnergy([time_series, block]) Get the kinetic energy.
getName() Return the process name.
getNonBondedEnergy14([time_series, block]) Get the non-bonded energy between atoms 1 and 4.
getOutput([name, block, file_link]) Return a link to a zip file of the working directory.
getPackageName() Return the package name.
getPotentialEnergy([time_series, block]) Get the potential energy.
getPressure([time_series, block]) Get the pressure.
getRecord(record[, time_series, unit, block]) Get a record from the stdout dictionary.
getRecords([block]) Return the dictionary of stdout time-series records.
getRestraintEnergy([time_series, block]) Get the restraint energy.
getSeed() Return the random number seed.
getStderr([block]) Return the entire stderr for the process as a list of strings.
getStdout([block]) Return the entire stdout for the process as a list of strings.
getStep([time_series, block]) Get the number of integration steps.
getSystem([block]) Get the latest molecular system.
getTemperature([time_series, block]) Get the temperature.
getTime([time_series, block]) Get the time (in nanoseconds).
getTotalEnergy([time_series, block]) Get the total energy.
getTrajectory([block]) Return a trajectory object.
getVanDerWaalsEnergy([time_series, block]) Get the Van der Vaals energy.
getVirial([time_series, block]) Get the virial.
getVolume([time_series, block]) Get the volume.
inputFiles() Return the list of input files.
insertArg(arg, value, index) Insert a command-line argument at a specific index.
isError() Return whether the process errored.
isQueued() Return whether the process is queued.
isRunning() Return whether the process is running.
kill() Kill the running process.
resetArgs() Reset the command-line arguments.
resetConfig() Reset the configuration parameters.
run([system, protocol, autostart, restart]) Create and run a new process.
runTime() Return the running time for the process (in minutes).
setArg(arg, value) Set a specific command-line argument.
setArgs(args) Set the dictionary of command-line arguments.
setConfig(config) Set the list of configuration file strings.
setName(name) Set the process name.
setSeed(seed) Set the random number seed.
start() Start the AMBER process.
stderr([n]) Print the last n lines of the stderr buffer.
stdout([n]) Print the last n lines of the stdout buffer.
wait([max_time]) Wait for the process to finish.
workDir() Return the working directory.
writeConfig(file) Write the configuration to file.
addArgs(args)

Append additional command-line arguments.

Parameters:args (dict, collections.OrderedDict) – A dictionary of command line arguments.
addToConfig(config)

Add a string to the configuration list.

Parameters:config (str, [ str ]) – A configuration string, a list of configuration strings, or a path to a configuration file.
clearArgs()

Clear all of the command-line arguments.

command()

Return the command-line string used to run the process.

Returns:command – The command string.
Return type:str
deleteArg(arg)

Delete an argument from the dictionary.

Parameters:arg (str) – The argument to delete.
exe()

Return the executable.

Returns:exe – The path to the executable.
Return type:str
getAngleEnergy(time_series=False, block='AUTO')

Get the angle energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The angle energy.

Return type:

Energy

getArgString()

Get the command-line arguments string.

Returns:arg_string – The command-line argument string.
Return type:str
getArgStringList()

Convert the argument dictionary into a list of strings.

Returns:arg_string_list – The list of command-line arguments.
Return type:[str]
getArgs()

Get the dictionary of command-line arguments.

Returns:args – The dictionary of command-line arguments.
Return type:collections.OrderedDict
getBondEnergy(time_series=False, block='AUTO')

Get the bond energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The bond energy.

Return type:

Energy

getCentreOfMassKineticEnergy(time_series=False, block='AUTO')

Get the kinetic energy of the centre of mass in translation.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The centre of mass kinetic energy.

Return type:

Energy

getConfig()

Get the list of configuration file strings.

Returns:config – The list of configuration strings.
Return type:[str]
getCurrentAngleEnergy(time_series=False)

Get the current angle energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The angle energy.
Return type:Energy
getCurrentBondEnergy(time_series=False)

Get the current bond energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The bond energy.
Return type:Energy
getCurrentCentreOfMassKineticEnergy(time_series=False)

Get the current kinetic energy of the centre of mass in translation.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The centre of mass kinetic energy.
Return type:Energy
getCurrentDensity(time_series=False)

Get the current density.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:density – The density.
Return type:float
getCurrentDihedralEnergy(time_series=False)

Get the current dihedral energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The dihedral energy.
Return type:Energy
getCurrentElectrostaticEnergy(time_series=False)

Get the current dihedral energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The electrostatic energy.
Return type:Energy
getCurrentElectrostaticEnergy14(time_series=False)

Get the current electrostatic energy between atoms 1 and 4.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The electrostatic energy between atoms 1 and 4.
Return type:Energy
getCurrentHydrogenBondEnergy(time_series=False)

Get the current hydrogen bond energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The hydrogen bond energy.
Return type:Energy
getCurrentKineticEnergy(time_series=False)

Get the current kinetic energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The kinetic energy.
Return type:Energy
getCurrentNonBondedEnergy14(time_series=False)

Get the current non-bonded energy between atoms 1 and 4.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The non-bonded energy between atoms 1 and 4.
Return type:Energy
getCurrentPotentialEnergy(time_series=False)

Get the current potential energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The potential energy.
Return type:Energy
getCurrentPressure(time_series=False)

Get the current pressure.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:pressure – The pressure.
Return type:Pressure
getCurrentRecord(record, time_series=False, unit=None)

Get a current record from the stdout dictionary.

Parameters:
  • record (str) – The record key.
  • time_series (bool) – Whether to return a list of time series records.
  • unit (Unit) – The unit to convert the record to.
Returns:

record – The matching record.

Return type:

Type

getCurrentRecords()

Return the current dictionary of stdout time-series records.

Returns:records – The dictionary of time-series records.
Return type:MultiDict
getCurrentRestraintEnergy(time_series=False)

Get the current restraint energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The restraint energy.

Return type:

Energy

getCurrentStep(time_series=False)

Get the current number of integration steps.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:step – The current number of integration steps.
Return type:int
getCurrentSystem()

Get the latest molecular system.

Returns:system – The latest molecular system.
Return type:System
getCurrentTemperature(time_series=False)

Get the current temperature.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:temperature – The temperature.
Return type:Temperature
getCurrentTime(time_series=False)

Get the current time (in nanoseconds).

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:time – The current simulation time in nanoseconds.
Return type:Time
getCurrentTotalEnergy(time_series=False)

Get the current total energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The total energy.
Return type:Energy
getCurrentVanDerWaalsEnergy(time_series=False)

Get the current Van der Vaals energy.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:energy – The Van der Vaals energy.
Return type:Energy
getCurrentVirial(time_series=False)

Get the current virial.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:virial – The virial.
Return type:float
getCurrentVolume(time_series=False)

Get the current volume.

Parameters:time_series (bool) – Whether to return a list of time series records.
Returns:volume – The volume.
Return type:Volume
getDensity(time_series=False, block='AUTO')

Get the density.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

density – The density.

Return type:

float

getDihedralEnergy(time_series=False, block='AUTO')

Get the dihedral energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The dihedral energy.

Return type:

Energy

getElectrostaticEnergy(time_series=False, block='AUTO')

Get the electrostatic energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The electrostatic energy.

Return type:

Energy

getElectrostaticEnergy14(time_series=False, block='AUTO')

Get the electrostatic energy between atoms 1 and 4.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The electrostatic energy between atoms 1 and 4.

Return type:

Energy

getHydrogenBondEnergy(time_series=False, block='AUTO')

Get the hydrogen bond energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The hydrogen bond energy.

Return type:

Energy

getKineticEnergy(time_series=False, block='AUTO')

Get the kinetic energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The kinetic energy.

Return type:

Energy

getName()

Return the process name.

Returns:

name : str
The name of the process.
getNonBondedEnergy14(time_series=False, block='AUTO')

Get the non-bonded energy between atoms 1 and 4.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The non-bonded energy between atoms 1 and 4.

Return type:

Energy

getOutput(name=None, block='AUTO', file_link=False)

Return a link to a zip file of the working directory.

Parameters:
  • name (str) – The name of the zip file.
  • block (bool) – Whether to block until the process has finished running.
  • file_link (bool) – Whether to return a FileLink when working in Jupyter.
Returns:

ouput – A path, or file link, to an archive of the process output.

Return type:

str, IPython.display.FileLink

getPackageName()

Return the package name.

Returns:name – The name of the package.
Return type:str
getPotentialEnergy(time_series=False, block='AUTO')

Get the potential energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The potential energy.

Return type:

Energy

getPressure(time_series=False, block='AUTO')

Get the pressure.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

pressure – The pressure.

Return type:

Pressure

getRecord(record, time_series=False, unit=None, block='AUTO')

Get a record from the stdout dictionary.

Parameters:
  • record (str) – The record key.
  • time_series (bool) – Whether to return a list of time series records.
  • unit (Unit) – The unit to convert the record to.
  • block (bool) – Whether to block until the process has finished running.
Returns:

record – The matching record.

Return type:

Type

getRecords(block='AUTO')

Return the dictionary of stdout time-series records.

Parameters:block (bool) – Whether to block until the process has finished running.
Returns:records – The dictionary of time-series records.
Return type:MultiDict
getRestraintEnergy(time_series=False, block='AUTO')

Get the restraint energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The restraint energy.

Return type:

Energy

getSeed()

Return the random number seed.

Returns:seed – The random number seed.
Return type:int
getStderr(block='AUTO')

Return the entire stderr for the process as a list of strings.

Parameters:block (bool) – Whether to block until the process has finished running.
Returns:error – The list of stderror strings.
Return type:[str]
getStdout(block='AUTO')

Return the entire stdout for the process as a list of strings.

Parameters:block (bool) – Whether to block until the process has finished running.
Returns:output – The list of stdout strings.
Return type:[str]
getStep(time_series=False, block='AUTO')

Get the number of integration steps.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

step – The current number of integration steps.

Return type:

int

getSystem(block='AUTO')

Get the latest molecular system.

Parameters:block (bool) – Whether to block until the process has finished running.
Returns:system – The latest molecular system.
Return type:System
getTemperature(time_series=False, block='AUTO')

Get the temperature.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

temperature – The temperature.

Return type:

Temperature

getTime(time_series=False, block='AUTO')

Get the time (in nanoseconds).

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

time – The current simulation time in nanoseconds.

Return type:

Time

getTotalEnergy(time_series=False, block='AUTO')

Get the total energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The total energy.

Return type:

Energy

getTrajectory(block='AUTO')

Return a trajectory object.

Parameters:block (bool) – Whether to block until the process has finished running.
Returns:trajectory – The latest trajectory object.
Return type:Trajectory
getVanDerWaalsEnergy(time_series=False, block='AUTO')

Get the Van der Vaals energy.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

energy – The Van der Vaals energy.

Return type:

Energy

getVirial(time_series=False, block='AUTO')

Get the virial.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

virial – The virial.

Return type:

float

getVolume(time_series=False, block='AUTO')

Get the volume.

Parameters:
  • time_series (bool) – Whether to return a list of time series records.
  • block (bool) – Whether to block until the process has finished running.
Returns:

volume – The volume.

Return type:

Volume

inputFiles()

Return the list of input files.

Returns:input_files – The list of autogenerated input files.
Return type:[str]
insertArg(arg, value, index)

Insert a command-line argument at a specific index.

Parameters:
  • arg (str) – The argument to set.
  • value – The value of the argument.
  • index (int) – The index in the dictionary.
isError()

Return whether the process errored.

Returns:is_error – Whether the process errored.
Return type:bool
isQueued()

Return whether the process is queued.

Returns:is_queued – Whether the process is queued.
Return type:bool
isRunning()

Return whether the process is running.

Returns:is_running – Whether the process is running.
Return type:bool
kill()

Kill the running process.

resetArgs()

Reset the command-line arguments.

resetConfig()

Reset the configuration parameters.

run(system=None, protocol=None, autostart=True, restart=False)

Create and run a new process.

Parameters:
  • system (System) – The molecular system.
  • protocol (Protocol) – The simulation protocol.
  • autostart (bool) – Whether to start the process automatically.
  • restart (bool) – Whether to restart the simulation, i.e. use the original system.
Returns:

process – The new process object.

Return type:

Procees

runTime()

Return the running time for the process (in minutes).

Returns:runtime – The runtime in minutes.
Return type:Time
setArg(arg, value)

Set a specific command-line argument.

For command-line flags, i.e. boolean arguments, the key should specify whether the flag is enabled (True) or not (False).

Parameters:
  • arg (str) – The argument to set.
  • value (bool, str) – The value of the argument.
setArgs(args)

Set the dictionary of command-line arguments.

Parameters:args (dict, collections.OrderedDict) – A dictionary of command line arguments.
setConfig(config)

Set the list of configuration file strings.

Parameters:config (str, [str]) – The list of configuration strings, or a path to a configuration file.
setName(name)

Set the process name.

Parameters:name (str) – The process name.
setSeed(seed)

Set the random number seed.

Parameters:seed (int) – The random number seed.
start()

Start the AMBER process.

Returns:process – The process object.
Return type:Process.Amber
stderr(n=10)

Print the last n lines of the stderr buffer.

Parameters:n (int) – The number of lines to print.
stdout(n=10)

Print the last n lines of the stdout buffer.

Parameters:n (int) – The number of lines to print.
wait(max_time=None)

Wait for the process to finish.

Parameters:max_time (Time, int, float) – The maximimum time to wait (in minutes).
workDir()

Return the working directory.

Returns:work_dir – The path of the working directory.
Return type:str
writeConfig(file)

Write the configuration to file.

Parameters:file (str) – The path to a file.