Trusted, Bureaucrats, Administrators
1,660
edits
Changes
Jump to navigation
Jump to search
#hRGS7:13@N 0.8142253659 0.0290794774
#hRGS7:14@N 0.7922308188 0.0801096067
</source>
==== Define the magnetic dipole-dipole relaxation interaction ====
In the sample script, we see that we should use the command "interatom.define".
'''Use the help command to find out more.'''
Add this
<source lang="python">
# Define the magnetic dipole-dipole relaxation interaction.
interatom.define(spin_id1='@N', spin_id2='@H', direct_bond=True)
interatom.define(spin_id1='@NE1', spin_id2='@HE1', direct_bond=True)
interatom.set_dist(spin_id1='@N*', spin_id2='@H*', ave_dist=1.02 * 1e-10)
interatom.unit_vectors()
# Define the chemical shift relaxation interaction.
value.set(-172 * 1e-6, 'csa', spin_id='@N*')
</source>
Now we are just going to try this.
<source lang="bash">
relax sam_script.py
</source>
We see that relax calculates a lot of stuff, for spins we have no data for.
'''Let us practice to select and deselect spins.'''
So instead we do this:
<source lang="python">
# Deselect
deselect.all()
# Select
#select.spin(':12@N|:12@H|:13@N|:13@H|:14@N|:14@H|:15@N|:15@H')
select.spin(':12|:13|:14|:15')
for spin, mol_name, res_num, res_name, spin_id in spin_loop(full_info=True, return_id=True, skip_desel=True):
print spin, mol_name, res_num, res_name, spin_id
# Define the magnetic dipole-dipole relaxation interaction.
interatom.define(spin_id1='@N', spin_id2='@H', direct_bond=True)
interatom.define(spin_id1='@NE1', spin_id2='@HE1', direct_bond=True)
interatom.set_dist(spin_id1='@N*', spin_id2='@H*', ave_dist=1.02 * 1e-10)
interatom.unit_vectors()
# Define the chemical shift relaxation interaction.
value.set(-172 * 1e-6, 'csa', spin_id='@N*')
</source>
