Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.

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.

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.

Unraveling effects of disorder on the electronic structure of SiO(2) from first principles

We present a first-principles systematic study of the electronic structure of SiO(2) including the crystalline polymorphs alpha quartz and beta cristobalite, and different types of disorder leading to the amorphous phase. We start from calculations within density functional theory and proceed to more sophisticated quasiparticle calculations according to the GW scheme. Our results show that different origins of disorder have also different impact on atomic and electronic-density fluctuations, which affect the electronic structure and, in particular, the size of the mobility gap in each case.

Anomalous lattice parameter of magnetic semiconductor alloys

The addition of transition metals to III-V semiconductors radically changes their electronic, magnetic, and structural properties. We show by ab initio calculations that in contrast to the conventional semiconductor alloys, the lattice parameter in magnetic semiconductor alloys, including those with diluted concentration, strongly deviates from Vegard's law. We find a direct correlation between the magnetic moment and the anion-transition metal bond lengths and derive a simple and general formula that determines the lattice parameter of a particular magnetic semiconductor by considering both the composition and magnetic moment. This dependence can explain some experimentally observed anomalies and stimulate other kind of investigations.

Gradient-dependent density functionals of the Perdew-Burke-Ernzerhof type for atoms, molecules, and solids

One of the standard generalized-gradient approximations (GGAs) in use in modern electronic-structure theory [Perdew-Burke-Ernzerhof (PBE) GGA] and a recently proposed modification designed specifically for solids (PBEsol) are identified as particular members of a family of functionals taking their parameters from different properties of homogeneous or inhomogeneous electron liquids. Three further members of this family are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules, and solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, is not always the best performing member of the family and that PBEsol, in spite of having been constructed specifically for solids, is not the best for solids. The performance of GGAs for finite systems is found to sensitively depend on the choice of constraints stemming from infinite systems. Guidelines both for users and for developers of density functionals emerge from this work.

Electronic, structural, and transport properties of Ni-doped graphene nanoribbons

We have investigated the electronic and transport properties of zigzag Ni-adsorbed graphene nanoribbons (Ni/GNRs) using ab initio calculations. We find that the Ni adatoms lying along the edge of zigzag GNRs represent the energetically most stable configuration, with an energy difference of approximately 0.3 eV when compared to the adsorption in the middle of the ribbon. The carbon atoms at the ribbon edges still present nonzero magnetic moments as in the pristine GNR even though there is a quenching by a factor of almost five in the value of the local magnetic moments at the C atoms bonded to the Ni. This quenching decays relatively fast and at approximately 9 A from the Ni adsorption site the magnetic moments have already values close to the pristine ribbon. At the opposite edge and at the central carbon atoms the changes in the magnetic moments are negligible. The energetic preference for the antiparallel alignment between the magnetization at the opposite edges of the ribbon is still maintained upon Ni adsorption. We find many Ni d-related states within an energy window of 1 eV above and below the Fermi energy, which gives rise to a spin-dependent charge transport. These results suggest the possibility of manufacturing spin devices based on GNRs doped with Ni atoms.

Designing real nanotube-based gas sensors

Using a combination of density functional theory and recursive Green's functions techniques, we present a full description of a large scale sensor, accounting for disorder and different coverages. Here, we use this method to demonstrate the functionality of nitrogen-rich carbon nanotubes as ammonia sensors as an example. We show how the molecules one wishes to detect bind to the most relevant defects on the nanotube, describe how these interactions lead to changes in the electronic transport properties of each isolated defect, and demonstrate that there are significative resistance changes even in the presence of disorder, elucidating how a realistic nanosensor works.

Comparative study of the adsorption and dissociation of vinylacetic acid and acrylic acid on silicon (001)

In this work, we employ the state of the art pseudopotential method, within a generalized gradient approximation to the density functional theory, to investigate the adsorption process of acrylic acid (AAc) and vinylacetic acid (VAA) on the silicon surface. Our total energy calculations support the proposed experimental process, as it indicates that the chemisorption of the molecule is as follows: The gas phase VAA (AAc) adsorbs molecularly to the electrophilic surface Si atom and then dissociates into H(2)C = CH - COO and H, bonded to the electrophilic and nucleophilic surface silicon dimer atoms, respectively. The activation energy for both processes correspond to thermal activations that are smaller than the usual growth temperature. In addition, the electronic structure, calculated vibrational modes, and theoretical scanning tunneling microscopy images are discussed, with a view to contribute to further experimental investigations.

Elementary Excitations at Magnetic Surfaces and Their Spin Dependence

The elementary surface excitations are studied by spin-polarized electron energy loss spectroscopy on a prototype oxide surface [an oxygen passivated Fe(001)-p(1 x 1) surface], where the various excitations coexist. For the first time, the surface phonons and magnons are measured simultaneously and are distinguished based on their different spin nature. The dispersion relation of all excitations is probed over the entire Brillouin zone. The different phonon modes observed in our experiment are described by means of ab initio calculations.

Two-photon absorption spectra of carotenoids compounds

Carotenoids are biosynthetic organic pigments that constitute an important class of one-dimensional pi-conjugated organic molecules with enormous potential for application in biophotonic devices. In this context, we studied the degenerate two-photon absorption (2PA) cross-section spectra of two carotenoid compounds (beta-carotene and beta-apo-8'-carotenal) employing the conventional and white-light-continuum Z-scan techniques and quantum chemistry calculations. Because carotenoids coexist at room temperature as a mixture of isomers, the 2PA spectra reported here are due to samples containing a distribution of isomers, presenting distinct conjugation length and conformation. We show that these compounds present a defined structure on the 2PA spectra, that peaks at 650 nm with an absorption cross-section of approximately 5000 GM, for both compounds. In addition, we observed a 2PA band at 990 nm for beta-apo-8'-carotenal, which was attributed to a overlapping of I(I)B(u) +-like and 2(I)Ag(-)-like states, which are strongly one- and two-photon allowed, respectively. Spectroscopic parameters of the electronic transitions to singlet-excited states, which are directly related to photophysical properties of these compounds, were obtained by fitting the 2PA spectra using the sum-over-states approach. The analysis and interpretations of the 2PA spectra of the investigated carotenoids were supported by theoretical predictions of one- and two-photon transitions carried out using the response functions formalism within the density functional theory framework, using the long-range corrected CAM-B3LYP functional. (C) 2011 American Institute of Physics. [doi:10.1063/1.3590157]

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]

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.

Ground-state functional for quantum spin chains with bond defects

We have developed a nonlocal functional of the exchange interaction for the ground-state energy of quantum spin chains described by the Heisenberg Hamiltonian. An alternating chain is used to obtain the correlation energy and a local unit-cell approximation is defined in the context of the density-functional theory. The agreement with our exact numerical data, for small chains, is significantly better than a previous formulation, even for chains with several ferromagnetic or antiferromagnetic bond defects. The results can be particularly relevant in the study of finite spin-1/2 Heisenberg chains, with exchange couplings changing, magnitude, or even sign, from bond-to-bond.