Role of the range of the dipole function in the classical dynamics of molecular dissociation

The dissociation dynamics of heteronuclear diatomic molecules induced by infrared laser pulses is investigated within the framework of the classical driven Morse oscillator. The interaction between the molecule and the laser field described in the dipole formulation is given by the product of a time-dependent external field with a position-dependent permanent dipole function. The effects of changing the spatial range of the dipole function in the classical dissociation dynamics of large ensembles of trajectories are studied. Numerical calculations have been performed for distinct amplitudes and carrier frequencies of the external pulses and also for ensembles with different initial energies. It is found that there exist a set of values of the dipole range for which the dissociation probability can be completely suppressed. The dependence of the dissociation on the dipole range is explained through the examination of the Fourier series coefficients of the dipole function in the angle variable of the free system. In particular, the suppression of dissociation corresponds to dipole ranges for which the Fourier coefficients associated with nonlinear resonances are null and the chaotic region in the phase space is reduced to thin layers. In this context, it is shown that the suppression of dissociation of heteronuclear molecules for certain frequencies of the external field is a consequence of the finite range of the corresponding permanent dipole. © 2013 American Physical Society.

Quantum-classical correspondence and the role of the dipole function in molecular dissociation

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

Quantum atom optics with fermions from molecular dissociation

We study a fermionic atom optics counterpart of parametric down-conversion with photons. This can be realized through dissociation of a Bose-Einstein condensate of molecular dimers consisting of fermionic atoms. We present a theoretical model describing the quantum dynamics of dissociation and find analytic solutions for mode occupancies and atomic pair correlations, valid in the short time limit. The solutions are used to identify upper bounds for the correlation functions, which are applicable to any fermionic system and correspond to ideal particle number-difference squeezing.

Controle de dissociação molecular com ferramentas de dinâmica não linear

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

A Computational Analysis of the Collision-Induced Dissociation Mechanisms of the Dipeptide Glycine-Serine Under Mass Spectrometry

Collision-induced dissociation (CID) of peptides using tandem mass spectrometry (MS) has been used to determine the identity of peptides and other large biological molecules. Mass spectrometry (MS) is a useful tool for determining the identity of molecules based on their interaction with electromagnetic fields. If coupled with another method like infrared (IR) vibrational spectroscopy, MS can provide structural information, but in its own right, MS can only provide the mass-to-charge (m/z) ratio of the fragments produced, which may not be enough information to determine the mechanism of the collision-induced dissociation (CID) of the molecule. In this case, theoretical calculations provide a useful companion for MS data and yield clues about the energetics of the dissociation. In this study, negative ion electrospray tandem MS was used to study the CID of the deprotonated dipeptide glycine-serine (Gly-Ser). Though negative ion MS is not as popular a choice as positive ion MS, studies by Bowie et al. show that it yields unique clues about molecular structure which complement positive ion spectroscopy, such as characteristic fragmentations like the loss of formaldehyde from the serine residue.2 The increase in the collision energy in the mass spectrometer alters the flexibility of the dipeptide backbone, enabling isomerizations (reactions not resulting in a fragment loss) and dissociations to take place. The mechanism of the CID of Gly-Ser was studied using two computational methods, B3LYP/6-311+G* and M06-2X/6-311++G**. The main pathway for molecular dissociation was analyzed in 5 conformers in an attempt to verify the initial mechanism proposed by Dr. James Swan after examination of the MS data. The results suggest that the loss of formaldehyde from serine, which Bowie et al. indicates is a characteristic of the presence of serine in a protein residue, is an endothermic reaction that is made possible by the conversion of the translational energy of the ion into internal energy as the ion collides with the inert collision gas. It has also been determined that the M06-2X functional¿s improved description of medium and long-range correlation makes it more effective than the B3LYP functional at finding elusive transition states. M06-2X also more accurately predicts the energy of those transition states than does B3LYP. A second CID mechanism, which passes through intermediates with the same m/z ratio as the main pathway for molecular dissociation, but different structures, including a diketopiperazine intermediate, was also studied. This pathway for molecular dissociation was analyzed with 3 conformers and the M06-2X functional, due to its previously determined effectiveness. The results suggest that the latter pathway, which meets the same intermediate masses as the first mechanism, is lower in overall energy and therefore a more likely pathway of dissociation than the first mechanism.

Electron collisions with the HCOOH...(H2O)n (n=1, 2 and 3) complexes in liquid phase: the influence on the π* shape resonance of formic acid

In this conference we report cross sections for elastic collisions of low-energy electrons with the HCOOH…(H2O)n complexes, with n = 1, 2 and 3. The scattering cross sections were computed with the Schwinger multichannel method [K. Takatsuka and V. McKoy, Phys. Rev. A 24 , 2473 (1981); Phys. Rev. A 30 , 1734 (1984)] with pseudopotentials [M. H. F. Bettega, L. G. Ferreira, and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993)] in the static-exchange and static-exchange plus polarization approximations, for energies from 0.5 eV to 6 eV. We considered some diÆerent hydrogen-bonded structures for the complexes that were generated with classical Monte Carlo simulations [K. Coutinho and S. Canuto, J. Chem. Phys. 113, 9132, (2000)]. The aim of this work is to investigate the effect of the surrounding water molecules on the π* shape resonance of the solute. Previous theoretical and experimental studies carried out in the gas phase reported a π* state for HCOOH at around 1.9 eV. For the n = 1 case and for all complexes, the stabilization of the resonance was observed (it appears at lower energy compared to the value obtained in the gas phase), as reported previously for the CH2O…H2O complexes [T. C. Freitas, M. A. P. Lima, S. Canuto, and M. H. F. Bettega, Phys. Rev. A 80, 062710 (2009)]. This result indicates that the presence of the solvent may affect the processes related to the π* state, such as the molecular dissociation by electron impact. For the n = 2 case we have observed both stabilization and destabilization of the π* resonance, that is associated with the hydrogen bond donor or acceptor role of the water molecules in the complexes. For the n = 3 case, preliminary static-exchange results show the stabilization of the π* state. We propose an explanation of the stabilization/destabilization of the π* state in terms of the polarization of the solute due to the surrounding water molecules and the net charge in the solute.

Numerical modelling of radiating superorbital flows

Experimental aerodynamic studies of the flows around new aerocapture spacecraft configurations are presently being done in the superorbital expansion tubes at The University of Queensland. Short duration flows at speeds of 10--13 km/s are produced in the expansion tube facility and are then applied to the model spacecraft. Although high-temperature effects, such as molecular dissociation, have long been a part of the computational modelling of the expansion tube flows for speeds below 10 km/s, radiation may now be a significant mechanism of energy transfer within the shock layer on the model. This paper will study the coupling of radiation energy transport for an optically thin gas to the flow dynamics in order to obtain accurate predictions of thermal loads on the spacecraft. The results show that the effect of radiation on the flowfields of subscale models for expansion tube experiments can be assessed by measurements of total heat transfer and radiative heat transfer.

Gaussian quantum Monte Carlo methods for fermions and bosons

We introduce a new class of quantum Monte Carlo methods, based on a Gaussian quantum operator representation of fermionic states. The methods enable first-principles dynamical or equilibrium calculations in many-body Fermi systems, and, combined with the existing Gaussian representation for bosons, provide a unified method of simulating Bose-Fermi systems. As an application relevant to the Fermi sign problem, we calculate finite-temperature properties of the two dimensional Hubbard model and the dynamics in a simple model of coherent molecular dissociation.

Water adsorption on O-covered Ru{0001}: coverage-dependent change from dissociation to molecular adsorption

When water is coadsorbed with oxygen at coverages above 0.25ML an intact water species is observed in high resolution X-ray photoelectron spectroscopy up to 220 K, which is significantly more stable than intact water on the clean surface. The presence of this species causes a shift in the O 1s binding energy of the pre-adsorbed oxygen, which indicates the formation of hydrogen bonds between the two adsorbates. Low coverages of oxygen induce partial dissociation and recombinative desorption in the same temperature range, which illustrates that desorption temperatures alone cannot be used to determine whether water is molecularly intact or not.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

Molecular masses and sedimentation coefficients of extracellular hemoglobin of Glossoscolex paulistus: Alkaline oligomeric dissociation

The giant extracellular hemoglobin of Glossoscolex paulistus (HbGp) has a molecular mass (M) of 3600 +/- 100 kDa and a standard sedimentation coefficient (s(20.w)(0)) of 58 S. estimated by analytical ultracentrifugation (AUC). In the present work, further AUC studies were developed for HbGp, at pH 10.0, which favors oligomeric dissociation into lower M species. The HbGp oligomer is formed by globin chains a, b, c and d plus the linker chains. The pure monomeric fraction, subunit d, and HbGp at pH 10.0, in the presence of beta-mercaptoethanol, were also studied. Our results indicate that for samples of pure subunit d, besides the monomeric species with s(20.w)(0) of 2.0 S, formation of dimer of subunit d is observed with s(20.w)(0) of around 2.9 S. For the whole HbGp at pH 10.0 contributions from monomers, trimers and linkers are observed. No contribution from 58 S species was observed for the sample of oxy-HbGp at pH 10.0, showing its complete dissociation. For cyanomet-HbGp form a contribution of 17% is observed for the un-dissociated oligomer, consistent with data from other techniques that show the cyanomet-form is more stable as compared to oxy-HbGp. Masses of HbGp subunits, especially trimer abc and monomeric chains a, b, c and d, were also estimated from sedimentation equilibrium data, and are in agreement with the results from MALDI-TOF-MS. (C) 2010 Elsevier B.V. All rights reserved.

Einstein-Podolsky-Rosen correlations via dissociation of a molecular Bose-Einstein condensate

Recent experimental measurements of atomic intensity correlations through atom shot noise suggest that atomic quadrature phase correlations may soon be measured with a similar precision. We propose a test of local realism with mesoscopic numbers of massive particles based on such measurements. Using dissociation of a Bose-Einstein condensate of diatomic molecules into bosonic atoms, we demonstrate that strongly entangled atomic beams may be produced which possess Einstein-Podolsky-Rosen (EPR) correlations in field quadratures in direct analogy to the position and momentum correlations originally considered by EPR.

Matter-wave amplification and phase conjugation via stimulated dissociation of a molecular Bose-Einstein condensate

We propose a scheme for parametric amplification and phase conjugation of an atomic Bose-Einstein condensate (BEC) via stimulated dissociation of a BEC of molecular dimers consisting of bosonic atoms. This can potentially be realized via coherent Raman transitions or using a magnetic Feshbach resonance. We show that the interaction of a small incoming atomic BEC with a (stationary) molecular BEC can produce two counterpropagating atomic beams - an amplified atomic BEC and its phase-conjugate or "time-reversed" replica. The two beams can possess strong quantum correlation in the relative particle number, with squeezed number-difference fluctuations.

First-principles quantum simulations of dissociation of molecular condensates: Atom correlations in momentum space

We investigate the quantum many-body dynamics of dissociation of a Bose-Einstein condensate of molecular dimers into pairs of constituent bosonic atoms and analyze the resulting atom-atom correlations. The quantum fields of both the molecules and atoms are simulated from first principles in three dimensions using the positive-P representation method. This allows us to provide an exact treatment of the molecular field depletion and s-wave scattering interactions between the particles, as well as to extend the analysis to nonuniform systems. In the simplest uniform case, we find that the major source of atom-atom decorrelation is atom-atom recombination which produces molecules outside the initially occupied condensate mode. The unwanted molecules are formed from dissociated atom pairs with nonopposite momenta. The net effect of this process-which becomes increasingly significant for dissociation durations corresponding to more than about 40% conversion-is to reduce the atom-atom correlations. In addition, for nonuniform systems we find that mode mixing due to inhomogeneity can result in further degradation of the correlation signal. We characterize the correlation strength via the degree of squeezing of particle number-difference fluctuations in a certain momentum-space volume and show that the correlation strength can be increased if the signals are binned into larger counting volumes.