Parameterizing a Simple Potential Energy Function Appropriate for Vibrational Frequencies Using Molecular Dynamics

A Maple scheme for quickly parameterizing vibrational potential energy functions is presented. As an example, the potential energy function's parameters for the vibrational motions in H_2O_2 are obtained assuming the simplest potential energy function. This paper was originally written as a research paper, but rejected by the referees. It is therefore being edited into an ``educational'' paper for student usage.

Universality of the self gravitational potential energy of any fundamental particle.

Using the relation proposed by Weinberg in 1972, combining quantum and cosmological parameters, we prove that the self gravitational potential energy of any fundamental particle is a quantum, with physical properties independent of the mass of the particle. It is a universal quantum of gravitational energy, and its physical properties depend only on the cosmological scale factor R and the physical constants ℏ and c. We propose a modification of the Weinberg’s relation, keeping the same numerical value, but substituting the cosmological parameter H/c by 1/R.

Raman and infrared spectra of dimethyl ether 13C-isotopologue (CH3O13CH3) from a CCSD(T) potential energy surface

So far, no experimental data of the infrared and Raman spectra of 13C isotopologue of dimethyl ether are available. With the aim of providing some clues of its low-lying vibrational bands and with the hope of contributing in a next spectral analysis, a number of vibrational transition frequencies below 300 cm−1 of the infrared spectrum and around 400 cm−1 of the Raman spectrum have been predicted and their assignments were proposed. Calculations were carried out through an ab initio three dimensional potential energy surface based on a previously reported one for the most abundant dimethyl ether isotopologue (M. Villa et al., J. Phys. Chem. A 115 (2011) 13573). The potential function was vibrationally corrected and computed with a highly correlated CCSD(T) method involving the COC bending angle and the two large amplitude CH3 internal rotation degrees of freedom. Also, the Hamiltonian parameters could represent a support for the spectral characterization of this species. Although the computed vibrational term values are expected to be very accurate, an empirical adjustment of the Hamiltonian has been performed with the purpose of anticipating some workable corrections to any possible divergence of the vibrational frequencies. Also, the symmetry breaking derived from the isotopic substitution of 13C in the dimethyl ether was taken into account when the symmetrization procedure was applied.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Total differential scattering cross section of {Ar^+}-Ar at 15 to 400 keV

We report on the measurement of the total differential scattering cross section of {Ar^+}-Ar at laboratory energies between 15 and 400 keV. Using an ab initio relativistic molecular program which calculates the interatomic potential energy curve with high accuracy, we are able to reproduce the detailed structure found in the experiment.

Theoretical evidence for quasi-molecular structure at small internuclear distances in elastic ion-atom scattering

The potential energy curve of the system Ne-Ne is calculated for small internuclear distances from 0.005 to 3.0 au using a newly developed relativistic molecular Dirac-Fock-Slater code. A significant structure in the potential energy curve is found which leads to a nearly complete agreement with experimental differential elastic scattering cross sections. This demonstrates the presence of quasi-molecular effects in elastic ion-atom collisions at keV energies.

Radiative association of Ti and O atoms

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Ab Initio correlated all electron dirac-fock calculations for Eka-Francium Fluoride (E119F)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Ab Initio correlated all electron Dirac-Fock calculations for eka-francium fluoride (E119F)

Resultados obtidos por cálculos Dirac-Fock correlacionados de 4 componentes para o fluoreto do elemento E119 (Eka-Frâncio) com base estável e precisa, livre de prolapso variacional, são reportados neste trabalho. No nível CCSD(T), a distância de equilíbrio R_e, frequência harmônica ω_e e energia de dissociação D_e são 2,432 Å, 354,97 cm_-1 e 116,92 kcal mol_-1, respectivamente. Também são reportados base livre de prolapso variacional de 4 componentes para o elemento 119, uma curva analítica de energia potencial precisa e o espectro vibracional a partir dos dados obtidos no nível CCSD(T). Nossos resultados sugerem que a molécula E119F deva ser menos iônica que seus fluoretos alcalinos homólogos mais leves, em contraste com o senso químico comum baseado nas propriedades periódicas - era de se esperar nesta molécula a ligação química mais iônica possível. Também encontramos que a correção do tipo modelo de carga para negligenciar as integrais do tipo SS resulta em erros insignificantes e acelera os cálculos cerca de 3 vezes no nível CCSD(T) e cerca de 4 vezes no nível DFT/B3LYP.

The formation of AlF by radiative association

The formation of the aluminium monofluoride molecule AlF by radiative association of the Al and F atoms is estimated. The radiative association of Al(P-2) and F(P-2) atoms is found to be dominated by the approach along the A(1) potential energy curve accompanied by spontaneous emission into the X-1 Sigma(+) ground state of the AlF. For temperatures ranging from 300 to 14 000 K, the rate coefficients are found to vary from 1.35 x 10(-17) to 9.31 x 10(-16) cm(3) s(-1), respectively. These values indicate that only a small amount of AlF molecules can be formed by radiative association in the inner envelope of carbon-rich stars and other hostile environments.

Radiative association of Ti and O atoms

Rate coefficients for the radiative association of titanium and oxygen atoms to form the titanium monoxide (TiO) molecule are estimated. The radiative association of Ti(F-3) and O(P-3) atoms is dominated by an approach along the C-3 Delta potential energy curve, accompanied by spontaneous emission into the X-3 Delta ground state of TiO. For temperatures ranging from 300-14 000 K, the total rate coefficients are found to vary from 4.76 x 10(-17) to 9.96 x 10(-17) cm(3) s(-1), respectively.

X-ray diffraction and theoretical investigation of the Gedunin crystal structure

The Gedunin compound (C28H34O6) is a natural product extracted from Trichilia pallida that has shown a wide activity. The crystallographic structure shows two conformers in the asymmetric unit, which differ in a rotation of the furan group. To understand this molecular arrangement, the density functional calculations. Molecular Electrostatic Potential (MEP) and thermodynamic function calculation have been performed at the B3LYP/6-311++g(d,p) level. Both conformers were optimized and the agreement with the experimental structure was very good, making possible further theoretical analysis of the structure. The inter-conversion between two conformers depends on the energy barrier. This process is studied in the vacuum and shows two transition states with a low energetic barrier for a potential energy curve scanning rigid around furan group: 4.37 kcal/mol and 16.52 kcal/mol. As the first transition state has a notably lower energetic barrier, the preferred inter-conversion pathway between the conformers involves the first rather than the second transition state. Understanding this transition state in detail led us to perform its optimization, showing an energetic barrier around 3.66 kcal/mol. The negative free energy and low enthalpy confirm that the process is spontaneous and exothermic. The results show that this requirement makes the existence of the two conformers in the asymmetric unit possible. The structure of molecules in the asymmetric unit is better understood when the MEP is used on the interaction between molecules. For Gedunin, both molecules have shown MEP with well-defined regions, and this behavior contributes to the observed link between molecules and for the negative regions complementing positive regions of another molecule. (C) 2011 Elsevier B.V. All rights reserved.

Dissociative electron attachment to chloromethane.

A low energy electron may attach to a molecule, forming a metastable resonance, which may dissociate into a stable anion and a neutral radical. Chloromethane has been a good target for dissociative electron attachment studies, since it is a small molecule with a clear dissociative ‘sigma*’ shape resonance. We present potential energy curves for CH3Cl and its anion, as a function of the C-Cl distance. Due to the resonant nature of the anion, a correct description requires a treatment based on scattering calculations. In order to compute elastic cross sections and phase shifts we employed the Schwinger multichannel method, implemented with pseudopotentials of Bachelet, Hamann and Schlüter, at the static-exchange plus polarization approximation. At the equilibrium geometry, the resonance was found arround 3.3 eV, in accordance to experience. The incoming electron is captured by a ‘sigma*’ orbital located at the C-Cl bond, which will relax in the presence of this extra electron. We took this bond as the reaction coordinate, and performed several scattering calculations for a series of nuclear conformations. The phase shift obtained in each calculation was fitted by a two component function, consisting in the usual Breit-Wigner profile, which captures the resonant character, and a second order polynomial in the wave number, which accounts for the background contribution. That way, we obtained position and width of the resonance, which allowed us to build the potential energy curve. For larger distances, the anion becomes stable and usual electronic structure calculations suffice. Furthermore, the existence of a dipole-bound anion state is revealed when we employed a set of very diffuse functions. The knowledge on the behaviour of the neutral and anionic electronic states helps us in elucidating how the dissociation takes place.

The formation of AlF by radiative association

The formation of the aluminium monofluoride molecule AlF by radiative association of the Al and F atoms is estimated. The radiative association of Al(2P) and F(2P) atoms is found to be dominated by the approach along theA1 potential energy curve accompanied by spontaneous emission into theX1 + ground state of the AlF. For temperatures ranging from 300 to 14 000 K, the rate coefficients are found to vary from 1.35×10−17 to 9.31×10−16 cm3 s−1, respectively.These values indicate that only a small amount of AlF molecules can be formed by radiative association in the inner envelope of carbon-rich stars and other hostile environments.