The anharmonic force field of HCN and HCP

The quadratic, cubic, and quartic force field of HCN has been calculated by a least squares refinement to fit the most recent observed data on the vibration-rotation constants of HCN, DCN and H13CN. All of the observed parameters are fitted within their standard errors of observation. The corresponding parameters for other isotopic species are calculated. For HCP and DCP the more limited data available have been fitted to an anharmonic force field using constraints based on comparison with HCN. Using this force field the zero-point rotational constants B0 have been corrected to obtain the equilibrium constants Be, and hence the equilibrium structure has been determined to be re(CH) = 1•0692(7)A, and re(CP) = 1•5398(2)A.

Vibration-rotation spectra and the harmonic force field of HOCl

Vibration-rotation spectra of HOCl have been measured at a resolution of 0.05 cm−1 to determine vibration rotation constants, and 35–37 Cl isotope shifts in the vibration frequencies. The spectrum of DOCl has also been recorded, and a preliminary analysis for the band origins has been made. The vibrational frequency data and centrifugal distortion constants have been used to determine the harmonic force field in a least-squares refinement; the force field obtained also gives a good fit to data on the vibrational contributions to the inertial defect. The equilibrium rotational constants of HOCl have been obtained, and an equilibrium structure has been estimated.

An anharmonic force field for SO3

An anharmonic force field for SO3 based on the valence force model has been investigated. The results of extending the model to include some further estimated cubic interaction potential constants have also been investigated. The phenomenological parameters calculated from both model force fields agree with those few values which have been experimentally determined. A calculation of the inertia defect has been made, and thus the value of C0 has been determined. The equilibrium structure has been determined to be: re = 1.4184 ± 0.0010 Å.

Nitrous acid: vibrational frequency shifts due to 15N substitution, and the harmonic force field

Vibration rotation spectra of HO15 NO and DO15 NO have been measured at a resolution of 0•04 cm-1 to determine the isotopic shifts in the vibrational band origins. These have been used together with recently determined data on the vibrational band origins, Coriolis constants, and centrifugal distorition constants, to determine the harmonic force field of both cis and trans nitrous acid in least squares refinement calculations. The results are discussed in relation to recent ab initio calculations, the inertia defects, and the torsional potential function.

The harmonic force field and rz structure of HNCO

The presently available microwave, millimeter wave, and far-infrared data of five isotopic species of isocyanic acid, namely, HNCO, H15NCO, HN13CO, HNC18O, and DNCO, have been used to obtain improved values of the ground-state rotational constants, the five quartic distortion constants, and some higher-order distortion constants in the Ir S reduced Hamiltonian of Watson. The appropriate planarity relation among the quartic centrifugal distortion constants has been imposed in the fitting procedure. The general harmonic force field of isocyanic acid has been determined using all existing data, and assuming a trans bent equilibrium geometry of the molecule with an NCO angle of 170°. Finally an rz structure has been obtained using the Az, Bz, and Cz rotational constants of five isotopic species. The bending of the NCO chain is found to be 8° in the trans configuration.

The infrared spectrum, equilibrium structure and harmonic and anharmonic force field of thioborine, HBS

Infrared spectra of the two stretching fundamentals of both HBS and DBS have been observed, using a continuous flow system through a multiple reflection long path cell at a pressure around 1 Torr and a Nicolet Fourier Transform spectrometer with a resolution of about 0•1 cm-1. The v3 BS stretching fundamental of DBS, near 1140 cm-1, is observed in strong Fermi resonance with the overtone of the bend 2v2. The bending fundamental v2 has not been observed and must be a very weak band. The analysis of the results in conjunction with earlier work gives the equilibrium structure (re(BH) = 1•1698(12) , re(BS) = 1•5978(3) ) and the harmonic and anharmonic force field.

Intensities in the v7, v8 and v9 bands of CH2NH, and the harmonic force field of methyleneimine

A high resolution Fourier transform infrared spectrum of methyleneimine, HN=CH2, has been obtained in the gas phase in the region 700 to 1300 cm−1. The rovibrational line intensities of the three lowest fundamentals ν7 (A′), ν8 (A″), and ν9 (A″) have been simulated including all Coriolis interactions between the three bands, and by fitting the observed spectrum the relative signs and magnitudes of the vibrational transition moments have been determined. All of the available spectroscopic data have been used to determine the harmonic force field of methyleneimine.

The general harmonic force field of methyl fluoride

The complete general harmonic force field of methyl flouride was recalculated using the most recent literature frequency, Coriolis ζ, and centrifugal distortion data for 12CH3F, 13CH3F, 12CD3F, 12CHD2F and 12CH2DF. The anharmonic corrections applied to the observed frequency data and the adopted molecular geometry are considered to be more realistic than those used hitherto. There is excellent overall agreement between the fitted force constants and the highest quality ab initio force field currently available.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Optimal space manoeuvres in a non-Keplerian force field

In this work the problem of a spacecraft bi-impulsive transfer between two given non coplanar elliptical orbits, with minimum fuel consumption, is solved considering a non-Keplerian force field (the perturbing forces include Earth gravity harmonics and atmospheric drag). The problem is transformed in the Two Point Boundary Value Problem. It is developed and implemented a new algorithm, that uses the analytical expressions developed here. A dynamics that considered a Keplerian force field was used to produce an initial guess to solve the Two Point Boundary Value Problem. Several simulations were performed to observe the spacecraft orbital behaviour by different kind of perturbations and constraints, on a fuel consumption optimization point of view. (C) 2002 COSPAR. Published by Elsevier B.V. Ltd. All rights reserved.

Ramachandran-type Plots for Glycosidic Linkages: Examples from Molecular Dynamics Simulations using the Glycam06 Force Field

The goals of this article are to (1) provide further validation of the Glycam06 force field, specifically for its use in implicit solvent molecular dynamic (MD) simulations, and (2) to present the extension of G.N. Ramachandran's idea of plotting amino acid phi and psi angles to the glycosidic phi, psi, and omega angles formed between carbohydrates. As in traditional Ramachandran plots, these carbohydrate Ramachandran-type (carb-Rama) plots reveal the coupling between the glycosidic angles by displaying the allowed and disallowed conformational space. Considering two-bond glycosidic linkages, there are 18 possible conformational regions that can be defined by (α, ϕ, ψ) and (β, ϕ, ψ), whereas for three-bond linkages, there are 54 possible regions that can be defined by (α, ϕ, ψ, ω) and (β, ϕ, ψ, ω). Illustrating these ideas are molecular dynamic simulations on an implicitly hydrated oligosaccharide (700 ns) and its eight constituent disaccharides (50 ns/disaccharide). For each linkage, we compare and contrast the oligosaccharide and respective disaccharide carb-Rama plots, validate the simulations and the Glycam06 force field through comparison to experimental data, and discuss the general trends observed in the plots.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Cost-Effective Force Field Tailored for Solid-Phase Simulations of OLED Materials

A united atom force field is empirically derived by minimizing the difference between experimental and simulated crystal cells and melting temperatures for eight compounds representative of organic electronic materials used in OLEDs and other devices: biphenyl, carbazole, fluorene, 9,9′-(1,3-phenylene)bis(9H-carbazole)-1,3-bis(N-carbazolyl)benzene (mCP), 4,4′-bis(N-carbazolyl)-1,1′-biphenyl (pCBP), phenazine, phenylcarbazole, and triphenylamine. The force field is verified against dispersion-corrected DFT calculations and shown to also successfully reproduce the crystal structure for two larger compounds employed as hosts in phosphorescent and thermally activated delayed fluorescence OLEDs: N,N′-di(1-naphthyl)-N,N′-diphenyl-(1,1′-biphenyl)-4,4′-diamine (NPD), and 1,3,5-tri(1-phenyl-1H-benzo[d]imidazol-2-yl)phenyl (TPBI). The good performances of the force field coupled to the large computational savings granted by the united atom approximation make it an ideal choice for the simulation of the morphology of emissive layers for OLED materials in crystalline or glassy phases.