Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

Combined Monte Carlo and quantum mechanics study of the solvatochromism of phenol in water. The origin of the blue shift of the lowest pi-pi* transition

A combined and sequential use of Monte Carlo simulations and quantum mechanical calculations is made to analyze the spectral shift of the lowest pi-pi* transition of phenol in water. The solute polarization is included using electrostatic embedded calculations at the MP2/aug-cc-pVDZ level giving a dipole moment of 2.25 D, corresponding to an increase of 76% compared to the calculated gas-phase value. Using statistically uncorrelated configurations sampled from the MC simulation,first-principle size-extensive calculations are performed to obtain the solvatochromic shift. Analysis is then made of the origin of the blue shift. Results both at the optimized geometry and in room-temperature liquid water show that hydrogen bonds of water with phenol promote a red shift when phenol is the proton-donor and a blue shift when phenol is the proton-acceptor. In the case of the optimized clusters the calculated shifts are in very good agreement with results obtained from mass-selected free jet expansion experiments. In the liquid case the contribution of the solute-solvent hydrogen bonds partially cancels and the total shift obtained is dominated by the contribution of the outer solvent water molecules. Our best result, including both inner and outer water molecules, is 570 +/- 35 cm(-1), in very good agreement with the small experimental shift of 460 cm(-1) for the absorption maximum.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Optical second harmonic generation in the centrosymmetric magnetic semiconductors EuTe and EuSe

The magnetic europium chalcogenide semiconductors EuTe and EuSe are investigated by the spectroscopy of second harmonic generation (SHG) in the vicinity of the optical band gap formed by transitions involving the 4f and 5d electronic orbitals of the magnetic Eu(2+) ions. In these materials with centrosymmetric crystal lattice the electric-dipole SHG process is symmetry forbidden so that no signal is observed in zero magnetic field. Signal appears, however, in applied magnetic field with the SHG intensity being proportional to the square of magnetization. The magnetic field and temperature dependencies of the induced SHG allow us to introduce a type of nonlinear optical susceptibility determined by the magnetic-dipole contribution in combination with a spontaneous or induced magnetization. The experimental results can be described qualitatively by a phenomenological model based on a symmetry analysis and are in good quantitative agreement with microscopic model calculations accounting for details of the electronic energy and spin structure.

Edge effects in bilayer graphene nanoribbons: Ab initio total-energy density functional theory calculations

We show that the ground state of zigzag bilayer graphene nanoribbons is nonmagnetic. It also possesses a finite gap, which has a nonmonotonic dependence with the width as a consequence of the competition between bulk and strongly attractive edge interactions. All results were obtained using ab initio total-energy density functional theory calculations with the inclusion of parametrized van der Waals interactions.

Spin filtering and disorder-induced magnetoresistance in carbon nanotubes: Ab initio calculations

Nitrogen-doped carbon nanotubes can provide reactive sites on the porphyrin-like defects. It is well known that many porphyrins have transition-metal atoms, and we have explored transition-metal atoms bonded to those porphyrin-like defects inN-doped carbon nanotubes. The electronic structure and transport are analyzed by means of a combination of density functional theory and recursive Green's function methods. The results determined the heme B-like defect (an iron atom bonded to four nitrogens) is the most stable and has a higher polarization current for a single defect. With randomly positioned heme B defects in nanotubes a few hundred nanometers long, the polarization reaches near 100%, meaning they are effective spin filters. A disorder-induced magnetoresistance effect is also observed in those long nanotubes, and values as high as 20 000% are calculated with nonmagnectic eletrodes.

Barrier-free substitutional doping of graphene sheets with boron atoms: Ab initio calculations

Using ab initio methods, we propose a simple and effective way to substitutionally dope graphene sheets with boron. The method consists of selectively exposing each side of the graphene sheet to different elements. We first expose one side of the membrane to boron while the other side is exposed to nitrogen. Proceeding this way, the B atoms will be spontaneously incorporated into the graphene membrane without any activation barrier. In a second step, the system should be exposed to a H-rich environment, which will remove the CN radical from the layer and form HCN, leading to a perfect substitutional doping.

Amorphous HfO(2) and Hf(1-x)Si(x)O via a melt-and-quench scheme using ab initio molecular dynamics

We have performed ab initio molecular dynamics simulations to generate an atomic structure model of amorphous hafnium oxide (a-HfO(2)) via a melt-and-quench scheme. This structure is analyzed via bond-angle and partial pair distribution functions. These results give a Hf-O average nearest-neighbor distance of 2.2 angstrom, which should be compared to the bulk value, which ranges from 1.96 to 2.54 angstrom. We have also investigated the neutral O vacancy and a substitutional Si impurity for various sites, as well as the amorphous phase of Hf(1-x)Si(x)O(2) for x=0.25, 0375, and 0.5.

Optical absorption and electronic band structure first-principles calculations of alpha-glycine crystals

Light absorption of alpha-glycine crystals grown by slow evaporation at room temperature was measured, indicating a 5.11 +/- 0.02 eV energy band gap. Structural, electronic, and optical absorption properties of alpha-glycine crystals were obtained by first-principles quantum mechanical calculations using density functional theory within the generalized gradient approximation in order to understand this result. To take into account the contribution of core electrons, ultrasoft and norm-conserving pseudopotentials, as well as an all electron approach were considered to compute the electronic density of states and band structure of alpha-glycine crystals. They exhibit three indirect energy band gaps and one direct Gamma-Gamma energy gap around 4.95 eV. The optical absorption related to transitions between the top of the valence band and the bottom of the conduction band involves O 2p valence states and C, O 2p conduction states, with the carboxyl group contributing significantly to the origin of the energy band gap. The calculated optical absorption is highly dependent on the polarization of the incident radiation due to the spatial arrangement of the dipolar glycine molecules; in the case of a polycrystalline sample, the first-principles calculated optical absorption is in good agreement with the measurement when a rigid energy shift is applied.

First-principles study of the adsorption of atomic and molecular hydrogen on BC(2)N nanotubes

The adsorption of atomic and molecular hydrogen on armchair and zigzag boron carbonitride nanotubes is investigated within the ab initio density functional theory. The adsorption of atomic H on the BC(2)N nanotubes presents properties which are promising for nanoelectronic applications. Depending on the adsorption site for the H, the Fermi energy moves toward the bottom of the conduction band or toward the top of the valence band, leading the system to exhibit donor or acceptor characteristics, respectively. The H(2) molecules are physisorbed on the BC(2)N surface for both chiralities. The binding energies for the H(2) molecules are slightly dependent on the adsorption site, and they are near to the range to work as a hydrogen storage medium.

Possible hyperdeformed band in (36)Ar observed in (12)C+(24)Mg elastic scattering

Fifteen strongly oscillating angular distributions of the elastic scattering of (12)C + (24)Mg at energies around the Coulomb barrier (E(c.m). = 10.67-16.00 MeV) are reproduced by adding five Breit-Wigner resonance terms to the l = 2, 4, 6, 7, and 8 elastic S matrix. The nonresonant, background elastic scattering S matrix S(l)(0) is calculated using the Sao Paulo potential. The J = 2, 4, 6, 7, and 8 (h) over bar molecular resonances fit well into a rotational molecular band, together with other higher lying resonances observed in the (16)O + (20)Ne elastic scattering. We propose that the presently observed, largely deformed molecular band corresponds to the hyperdeformed band, which has been found previously in alpha-cluster calculations, as well as in a new Nilsson model calculation. Systematic study of its possible clusterizations predicts the preference of the (12)C + (24)Mg and (16)O + (20)Ne molecular structure, in accordance with our present results.

A quasiclassical trajectory study of the OH plus SO reaction: The role of rotational energy

A full dimensional quasiclassical trajectory study of the OH+SO reaction is presented with the aim of investigating the role of the reactants rotational energy in the reactivity. Different energetic combinations with one and both reactants rotationally excited are studied. A passive method is used to correct zero-point-energy leakage in the classical calculations. The reactive cross sections, for each combination, are calculated and fitted to a capturelike model combined with a factor accounting for recrossing effects. Reactivity decreases as rotational energy is increased in any of both reactants. This fact provides a theoretical support for the experimental dependence of the rate constant on temperature.

Experimental and theoretical study of two-photon absorption in nitrofuran derivatives: Promising compounds for photochemotherapy

We report experimental and theoretical studies of the two-photon absorption spectrum of two nitrofuran derivatives: nitrofurantoine, (1-(5-nitro-2-furfurilideneamine)-hidantoine) and quinifuryl, 2-(5`-nitro-2`-furanyl) ethenyl-4-{N-[4`-(N,N-diethylamino)-1`-methylbutyl]carbamoyl} quinoline. Both molecules are representative of a family of 5-nitrofuran-ethenyl-quinoline drugs that have been demonstrated to display high toxicity to various species of transformed cells in the dark. We determine the two-photon absorption cross-section for both compounds, from 560 to 880 nm, which present peak values of 64 GM for quinifuryl and 20 GM for nitrofurantoine (1 GM = 1 x 10(-50) cm(4).s.photon(-1)). Besides, theoretical calculations employing the linear and quadratic response functions were carried out at the density functional theory level to aid the interpretations of the experimental results. The theoretical results yielded oscillator strengths, two-photon transition probabilities, and transition energies, which are in good agreement with the experimental data. A higher number of allowed electronic transitions was identified for quinifuryl in comparison to nitrofurantoine by the theoretical calculations. Due to the planar structure of both compounds, the differences in the two-photon absorption cross-section values are a consequence of their distinct conjugation lengths. (c) 2011 American Institute of Physics. [doi:10.1063/1.3514911]

Protecting a quantum state from environmental noise by an incompatible finite-time measurement

We show that measurements of finite duration performed on an open two-state system can protect the initial state from a phase-noisy environment, provided the measured observable does not commute with the perturbing interaction. When the measured observable commutes with the environmental interaction, the finite-duration measurement accelerates the rate of decoherence induced by the phase noise. For the description of the measurement of an observable that is incompatible with the interaction between system and environment, we have found an approximate analytical expression, valid at zero temperature and weak coupling with the measuring device. We have tested the validity of the analytical predictions against an exact numerical approach, based on the superoperator-splitting method, that confirms the protection of the initial state of the system. When the coupling between the system and the measuring apparatus increases beyond the range of validity of the analytical approximation, the initial state is still protected by the finite-time measurement, according with the exact numerical calculations.

Theoretical study of one- and two-photon absorption spectra of azoaromatic compounds

In this study, the one- and two-photon absorption spectra of seven azoaromatic compounds (five pseudostilbenes-type and two aminoazobenzenes) were theoretically investigated using the density functional theory combined with the response functions formalism. The equilibrium molecular structure of each compound was obtained at three different levels of theory: Hartree-Fock, density functional theory (DFT), and Moller-Plesset 2. The effect of solvent on the equilibrium structure and the electronic transitions of the compounds were investigated using the polarizable continuum model. For the one-photon absorption, the allowed pi ->pi(*) transition energy showed to be dependent on the molecular structures and the effect of solvent, while the n ->pi(*) and pi ->pi(*)(n) transition energies exhibited only a slight dependence. An inversion between the bands corresponding to the pi ->pi(*) and n ->pi(*) states due to the effect of solvent was observed for the pseudostilbene-type compounds. To characterize the allowed two-photon absorption transitions for azoaromatic compounds, the response functions formalism combined with DFT using the hybrid B3LYP and PBE0 functionals and the long-range corrected CAM-B3LYP functional was employed. The theoretical results support the previous findings based on the three-state model. The model takes into account the ground and two electronic excited states and has already been used to describe and interpret the two-photon absorption spectrum of azoaromatic compounds. The highest energy two-photon allowed transition for the pseudostilbene-type compounds shows to be more effectively affected (similar to 20%) by the torsion of the molecular structure than the lowest allowed transition (similar to 10%). In order to elucidate the effect of the solvent on the two-photon absorption spectra, the lowest allowed two-photon transition (dipolar transition) for each compound was analyzed using a two-state approximation and the polarizable continuum model. The results obtained reveal that the effect of solvent increases drastically the two-photon cross-section of the dipolar transition of the pseudostilbene-type compounds. In general, the features of both one- and two-photon absorption spectra of the azoaromatic compounds are well reproduced by the theoretical calculations.