pyemma.thermo.estimate_multi_temperature

pyemma.thermo.estimate_multi_temperature(energy_trajs, temp_trajs, dtrajs, energy_unit='kcal/mol', temp_unit='K', reference_temperature=None, maxiter=10000, maxerr=1e-15, save_convergence_info=0, estimator='wham', lag=1, dt_traj='1 step', init=None)

Wraps multi-temperature data.

Parameters:

• energy_trajs (list of N arrays, each of shape (T_i,)) – List of arrays, each having T_i rows, one for each time step, containing the potential energies time series in units of kT, kcal/mol or kJ/mol.
• temp_trajs (list of N int arrays, each of shape (T_i,)) – List of arrays, each having T_i rows, one for each time step, containing the heat bath temperature time series (at which temperatures the frames were created) in units of K or C. Alternatively, these trajectories may contain kT values instead of temperatures.
• dtrajs (list of N int arrays, each of shape (T_i,)) – The integers are indexes in 0,...,n-1 enumerating the n discrete states or the bins the trajectory is in at any time.
• energy_unit (str, optional, default='kcal/mol') – The physical unit used for energies. Current options: kcal/mol, kJ/mol, kT.
• temp_unit (str, optional, default='K') – The physical unit used for the temperature. Current options: K, C, kT
• reference_temperature (float or None, optional, default=None) – Reference temperature against which the bias energies are computed. If not given, the lowest temperature or kT value is used. If given, this parameter must have the same unit as the temp_trajs.
• maxiter (int, optional, default=10000) – The maximum number of self-consistent iterations before the estimator exits unsuccessfully.
• maxerr (float, optional, default=1E-15) – Convergence criterion based on the maximal free energy change in a self-consistent iteration step.
• save_convergence_info (int, optional, default=0) – Every save_convergence_info iteration steps, store the actual increment and the actual loglikelihood; 0 means no storage.
• estimator (str, optional, default='wham') –

  Specify one of the available estimators

  ‘wham’: use WHAM

  ‘dtram’: use the discrete version of TRAM
• lag (int, optional, default=1) – Integer lag time at which transitions are counted.
• dt_traj (str, optional, default='1 step') –

  Description of the physical time corresponding to the lag. May be used by analysis algorithms such as plotting tools to pretty-print the axes. By default ‘1 step’, i.e. there is no physical time unit. Specify by a number, whitespace and unit. Permitted units are (* is an arbitrary string):

  ‘fs’, ‘femtosecond*’

  ‘ps’, ‘picosecond*’

  ‘ns’, ‘nanosecond*’

  ‘us’, ‘microsecond*’

  ‘ms’, ‘millisecond*’

  ‘s’, ‘second*’
• init (str, optional, default=None) –

  Use a specific initialization for self-consistent iteration:

  None: use a hard-coded guess for free energies and Lagrangian multipliers

  ‘wham’: perform a short WHAM estimate to initialize the free energies

Returns:

_estimator – The requested estimator/model object, i.e., WHAM or DTRAM.

Return type:

MEMM or StationaryModel