Analytical potential energy function and spectroscopic parameters for the ground and excited states of NaH

set of energies at different internuclear distances for the ground electronic state and two excited electronic states of NaH molecule have been calculated using valence internally contracted multireference configuration interaction(MRCI) including Davidson correction and three basis sets. Then, a potential energy curve (PEC) for each state was determined by extrapolating MRCI energies to the complete basis sets limit. Based on the PECs, accurate vibrational energy levels and rotational constants were determined. The computational PECs are were fitted to analytical potential energy functions using the Murrell-Sorbie potential function. Then, accurate spectroscopic parameters were calculated. Compared with experimental results, values obtained with the basis set extrapolation yield a potential energy curve that gives accurate vibrational energy levels, rotational constants and spectroscopic parameters for the NaH molecule. (C) 2004 Elsevier B.V. All rights reserved.

Potential energy landscape and robustness of a gene regulatory network: Toggle switch

Finding a multidimensional potential landscape is the key for addressing important global issues, such as the robustness of cellular networks. We have uncovered the underlying potential energy landscape of a simple gene regulatory network: a toggle switch. This was realized by explicitly constructing the steady state probability of the gene switch in the protein concentration space in the presence of the intrinsic statistical fluctuations due to the small number of proteins in the cell. We explored the global phase space for the system. We found that the protein synthesis rate and the unbinding rate of proteins to the gene were small relative to the protein degradation rate; the gene switch is monostable with only one stable basin of attraction. When both the protein synthesis rate and the unbinding rate of proteins to the gene are large compared with the protein degradation rate, two global basins of attraction emerge for a toggle switch. These basins correspond to the biologically stable functional states. The potential energy barrier between the two basins determines the time scale of conversion from one to the other. We found as the protein synthesis rate and protein unbinding rate to the gene relative to the protein degradation rate became larger, the potential energy barrier became larger. This also corresponded to systems with less noise or the fluctuations on the protein numbers.

Effect of quantization of vibrations on the structural properties of crystals

We study the structural effects produced by the quantization of vibrational degrees of freedom in periodic crystals at zero temperature. To this end we introduce a methodology based on mapping a suitable subspace of the vibrational manifold and solving the Schrödinger equation in it. A number of increasingly accurate approximations ranging from the quasiharmonic approximation (QHA) to the vibrational self-consistent field (VSCF) method and the exact solution are described. A thorough analysis of the approximations is presented for model monatomic and hydrogen-bonded chains, and results are presented for a linear H-F chain where the potential-energy surface is obtained via first-principles electronic structure calculations. We focus on quantum nuclear effects on the lattice constant and show that the VSCF is an excellent approximation, meaning that correlation between modes is not extremely important. The QHA is excellent for covalently bonded mildly anharmonic systems, but it fails for hydrogen-bonded ones. In the latter, the zero-point energy exhibits a nonanalytic behavior at the lattice constant where the H atoms center, which leads to a spurious secondary minimum in the quantum-corrected energy curve. An inexpensive anharmonic approximation of noninteracting modes appears to produce rather good results for hydrogen-bonded chains for small system sizes. However, it converges to the incorrect QHA results for increasing size. Isotope effects are studied for the first-principles H-F chain. We show how the lattice constant and the H-F distance increase with decreasing mass and how the QHA proves to be insufficient to reproduce this behavior.

Quasi-classical model of non-destructive wavepacket manipulation by intense ultrashort nonresonant laser pulses

A quasi-classical model (QCM) of nuclear wavepacket generation, modification and imaging by three intense ultrafast near-infrared laser pulses has been developed. Intensities in excess of 10(13) W cm(-2) are studied, the laser radiation is non-resonant and pulse durations are in the few-cycle regime, hence significantly removed from the conditions typical of coherent control and femtochemistry. The 1s sigma ground state of the D-2 precursor is projected onto the available electronic states in D-2(+) (1s sigma(g) ground and 2p sigma(u) dissociative) and D+ + D+ (Coulomb explosion) by tunnel ionization by an ultrashort 'pump' pulse, and relative populations are found numerically. A generalized non-adiabatic treatment allows the dependence of the initial vibrational population distribution on laser intensity to be calculated. The wavepacket is approximated as a classical ensemble of particles moving on the 1s sigma(g) potential energy surface (PES), and hence follow trajectories of different amplitudes and frequencies depending on the initial vibrational state. The 'control' pulse introduces a time-dependent polarization of the molecular orbital, causing the PES to be modified according to the dynamic Stark effect and the transition dipole. The trajectories adjust in amplitude, frequency and phase-offset as work is done on or by the resulting force; comparing the perturbed and unperturbed trajectories allows the final vibrational state populations and phases to be determined. The action of the 'probe' pulse is represented by a discrete internuclear boundary, such that elements of the ensemble at a larger internuclear separation are assumed to be photodissociated. The vibrational populations predicted by the QCM are compared to recent quantum simulations (Niederhausen and Thumm 2008 Phys. Rev. A 77 013404), and a remarkable agreement has been found. The applicability of this model to femtosecond and attosecond time-scale experiments is discussed and the relation to established femtochemistry and coherent control techniques are explored.

Hydrolysis characteristics and kinetics of waste hay biomass as a potential energy crop for fermentable sugars production using autoclave parr reactor system

The use of the organic fraction of municipal solid waste crops has received considerable attention as a sustainable feedstock that can replace fossil fuels for the production of renewable energy. Therefore, municipal bin-waste in the form of hay was investigated as a potential energy crop for fermentable sugars production. Hydrolysis of hay by dilute phosphoric acid was carried out in autoclave parr reactor, where reactor temperature (135-200 degrees c) and acid concentration (2.5-10% (w/w)) were examined. Analysis of the decomposition rate of hemicellulosic biomass was undertaken using HPLC of the reaction products. Xylose production reached a maximum value of 13.5 g/100 g dry mass corresponding to a yield of 67% at the best identified conditions of 2.5 wt% H3PO4, 175 degrees C, 10 min reaction time, and at 5 wt% H3PO4, 150 degrees C, and 5 min reaction time. For glucose, an average yield of 25% was obtained at 5 wt% H3PO4, 175 degrees C and 30 min. Glucose degradation to HMF was achieved at 10 wt% H3PO4 and 200 degrees C. The maximum yield for produced arabinose was an average of 3 g/100 g dry. mass corresponding to 100% of the total possible arabinose. The kinetic study of the acid hydrolysis was also carried out using the Saeman and the Two-fraction models. It was found for both models that the kinetic constants (k) depend on the acid concentration and temperature. For xylose and arabinose it was found that the rate of formation was more favoured than the rate of degradation. By contrast, for glucose it was found that glucose degradation was occurring faster than glucose formation. It can be concluded that dilute phosphoric acid hydrolysis of hay crop is feasible for the production of fermentable sugars which are essential for bioethanol synthesis. 

Exploring the free energy surfaces of clusters using reconnaissance metadynamics

A new approach is proposed for exploring the low-energy structures of small to medium-sized aggregates of atoms and molecules. This approach uses the recently proposed reconnaissance metadynamics method [G. A. Tribello, M. Ceriotti, and M. Parrinello. Proc. Natl. Acad. Sci. U.S.A. 107(41), 17509 (2010)] in tandem with collective variables that describe the average structure of the coordination sphere around the atoms/molecules. We demonstrate this method on both Lennard-Jones and water clusters and show how it is able to quickly find the global minimum in the potential energy surface, while exploring the finite temperature free energy surface. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3628676]

Effets de la concentration des défauts sur la surface d'énergie potentielle du silicium amorphe

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Caratérisation de la surface d’énergie potentielle des matériaux complexes et son application sur la cinétique du SiO2/Si

Dans ce rapport de mémoire, nous avons utilisé les méthodes numériques telles que la dynamique moléculaire (code de Lammps) et ART-cinétique. Ce dernier est un algorithme de Monte Carlo cinétique hors réseau avec construction du catalogue d'événements à la volée qui incorpore exactement tous les effets élastiques. Dans la première partie, nous avons comparé et évalué des divers algorithmes de la recherche du minimum global sur une surface d'énergie potentielle des matériaux complexes. Ces divers algorithmes choisis sont essentiellement ceux qui utilisent le principe Bell-Evans-Polanyi pour explorer la surface d'énergie potentielle. Cette étude nous a permis de comprendre d'une part, les étapes nécessaires pour un matériau complexe d'échapper d'un minimum local vers un autre et d'autre part de contrôler les recherches pour vite trouver le minimum global. En plus, ces travaux nous ont amené à comprendre la force de ces méthodes sur la cinétique de l'évolution structurale de ces matériaux complexes. Dans la deuxième partie, nous avons mis en place un outil de simulation (le potentiel ReaxFF couplé avec ART-cinétique) capable d'étudier les étapes et les processus d'oxydation du silicium pendant des temps long comparable expérimentalement. Pour valider le système mis en place, nous avons effectué des tests sur les premières étapes d'oxydation du silicium. Les résultats obtenus sont en accord avec la littérature. Cet outil va être utilisé pour comprendre les vrais processus de l'oxydation et les transitions possibles des atomes d'oxygène à la surface du silicium associée avec les énergies de barrière, des questions qui sont des défis pour l'industrie micro-électronique.

Nuclear curvature energy in relativistic models

The difficulties arising in the calculation of the nuclear curvature energy are analyzed in detail, especially with reference to relativistic models. It is underlined that the implicit dependence on curvature of the quantal wave functions is directly accessible only in a semiclassical framework. It is shown that also in the relativistic models quantal and semiclassical calculations of the curvature energy are in good agreement.