Package schrodinger :: Package application :: Package desmond :: Module remd

[frames] | no frames]

Module remd

Utilities for REMD.

Functions

[hide private]

erf(x)

is_water(site, constraint)

get_nonduplicated_atom(ct)

get_num_wateratom(model, selected_atom)

get_num_nonwateratom(model, selected_atom)

get_num_constraint(model, selected_atom)

get_degrees_of_freedom(model, selected_atom)

predict_temperature_ladder(temp, exchange_probability, model, asl, should_fix=True, floaty=False)
`temp' should be a two-tuple, specifying the temperature range.

predict_with_temp_and_exch(temp, exchange_probability, model, asl)
Given temperature range and exchange_probability, predicts the number of replica and returns the temperature ladder.

predict_with_nreplica_and_exch(n_replica, exchange_probability, base_temp, model, asl)
Given the base temperature, number of replicas, and exchange_probability, predicts the top temperature and returns the temperature ladder.

predict_with_temp_and_nreplica(temp, n_replica, model, asl)
Given the temperature range and number of replicas, predicts the exchange probability and returns the temperature ladder.

split_ct(ct, selected_atom)

set_freezing_atommass(model, mass_scale)

freeze_atom(model, asl, frozen_atom_mass_threshold=FROZEN_ATOM_MASS_THRESHOLD)

write_ff_solute_tempering(model, out_fname, asl, scaling_factor, should_scale_torsion=True)
write cms.Cms to a file with scaled force field terms

cms.Cms

rescale_ff_solute_tempering(model, asl, scaling_factor, should_scale_torsion=True)
return cms.Cms with scaled force field terms

rescale_ff(ct, rest_atoms, scaling_factor, should_scale_torsion=True)
scales force field terms for rest atoms

_rescale_charge(ct, rest_atoms, scaling_factor)
scales rest atom charge

_rescale_vdw(ct, rest_atoms, scaling_factor)
scales vdwtype and creates new vdwtype by appending _S for rest atoms

_rescale_dihedral(ct, rest_atoms, scaling_factor)
scales dihedral terms

test_get_num_wateratom(model)

test_get_num_nonwateratom(model)

test_get_num_constraint(model)

test_get_degrees_of_freedom(model)

test_predict_temperature(model)

test_freeze_atom(model)

test_predict_with_nreplica_and_exch(model)

test_predict_with_temp_and_nreplica(model)

Variables

[hide private]

FROZEN_ATOM_MASS_THRESHOLD = 1E9

model = cms.Cms(file= "test.cms")

Function Details

[hide private]

predict_temperature_ladder(temp, exchange_probability, model, asl, should_fix=True, floaty=False)

`temp' should be a two-tuple, specifying the temperature range. `model' should be a `Cms' object.

write_ff_solute_tempering(model, out_fname, asl, scaling_factor, should_scale_torsion=True)

write cms.Cms to a file with scaled force field terms

Parameters:

model (cms.Cms) - input cms.Cms
out_fname (str) - output file name
asl (str) - ASL expression to define rest atoms
scaling_factor (float) - T_ref/T_hot
should_scale_torsion (Boolean) - whether to scale dihedral terms

rescale_ff_solute_tempering(model, asl, scaling_factor, should_scale_torsion=True)

return cms.Cms with scaled force field terms

Parameters:

model (cms.Cms) - input cms.Cms
asl (str) - ASL expression to define rest atoms
scaling_factor (float) - T_ref/T_hot
should_scale_torsion (Boolean) - whether to scale dihedral terms

Returns: cms.Cms

rescaled cms.Cms

rescale_ff(ct, rest_atoms, scaling_factor, should_scale_torsion=True)

scales force field terms for rest atoms

Parameters:

ct (ffiostructure.FFIOStructure) - input/output structure handle
rest_atoms (list of int) - sequence of atom indices in rest region
scaling_factor (float) - T_ref/T_hot
should_scale_torsion (Boolean) - whether to scale dihedral terms

_rescale_charge(ct, rest_atoms, scaling_factor)


scales rest atom charge

@type ct: ffiostructure.FFIOStructure
@param ct: input/output structure handle
@type rest_atoms: list of int
@param rest_atoms: sequence of atom indices in rest region
@type scaling_factor: float
@param scaling_factor: T_ref/T_hot

Charge scaling in place
q_i *= sf_i
where
    sf_i = sqrt(scaling_factor) when i is a rest atom
    otherwise sf_i = 1.0

_rescale_vdw(ct, rest_atoms, scaling_factor)

scales vdwtype and creates new vdwtype by appending _S for rest atoms

Parameters:

ct (ffiostructure.FFIOStructure) - input/output structure handle
rest_atoms (list of int) - sequence of atom indices in rest region
scaling_factor (float) - T_ref/T_hot
vdw.t1 = scaling_factor

_rescale_dihedral(ct, rest_atoms, scaling_factor)


scales dihedral terms

@type ct: ffiostructure.FFIOStructure
@param ct: input/output structure handle
@type rest_atoms: list of int
@param rest_atoms: sequence of atom indices in rest region
@type scaling_factor: float
@param scaling_factor: T_ref/T_hot

Dihedral scaling
d.t1 = sf_i * sf_l
where
    d.t1 = scaling factor for a dihedral term with i, j, k, l atoms
    sf_i = sqrt(scaling_factor) when i is a rest atom
    sf_l = sqrt(scaling_factor) when l is a rest atom
    otherwise sf_i = sf_l = 1.0