Delocalized water and fluoride contributions to Dyson orbitals for electron detachment from the hydrated fluoride anion

The experimental vertical electron detachment energy (VEDE) of aqueous fluoride, [F(-)(H(2)O)], is approximately 9.8 eV, but spectral assignment is complicated by interference between F(-) 2p and H(2)O 1b(1) orbitals. The electronic structure of [F(-)(H(2)O)] is analyzed with Monte Carlo and ab initio quantum-mechanical calculations. Electron-propagator calculations in the partial third-order approximation yield a VEDE of 9.4 eV. None of the Dyson orbitals corresponding to valence VEDEs consists primarily of F 2p functions. These results and ground-state atomic charges indicate that the final, neutral state is more appropriately described as [F(-)(H(2)O)(+)] than as [F(H(2)O)]. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3431081]

Hyperpolarizabilities of the methanol molecule: A CCSD calculation including vibrational corrections

In this work we present the results for hyperpolarizabilities of the methanol molecule including vibrational corrections and electron correlation effects at the CCSD level. Comparisons to random phase approximation results previously reported show that the electron correlation is in general important for both electronic contribution and vibrational corrections. The role played by the anharmonicities on the calculations of the vibrational corrections has also been analyzed and the obtained results indicate that the anharmonic terms are important for the dc-Pockels and dc-Kerr effects. For the other nonlinear optical properties studied the double-harmonic approximation is found to be suitable. Comparison to available experimental result in gas phase for the dc-second harmonic generation second hyperpolarizability shows a very good agreement with the electronic contribution calculated here while our total value is 14% larger than the experimental value.

The Importance of Protein-Protein Interactions on the pH-Induced Conformational Changes of Bovine Serum Albumin: A Small-Angle X-Ray Scattering Study

The combined effects of concentration and pH on the conformational states of bovine serum albumin (BSA) are investigated by small-angle x-ray scattering. Serum albumins, at physiological conditions, are found at concentrations of similar to 35-45 mg/mL (42 mg/mL in the case of humans). In this work, BSA at three different concentrations (10, 25, and 50 mg/mL) and pH values (2.0-9.0) have been studied. Data were analyzed by means of the Global Fitting procedure, with the protein form factor calculated from human serum albumin (HSA) crystallographic structure and the interference function described, considering repulsive and attractive interaction potentials within a random phase approximation. Small-angle x-ray scattering data show that BSA maintains its native state from pH 4.0 up to 9.0 at all investigated concentrations. A pH-dependence of the absolute net protein charge is shown and the charge number per BSA is quantified to 10(2), 8(l), 13(2), 20(2), and 26(2) for pH values 4.0, 5.4, 7.0, 8.0, and 9.0, respectively. The attractive potential diminishes as BSA concentration increases. The coexistence of monomers and dimers is observed at 50 mg/mL and pH 5.4, near the BSA isoelectric point. Samples at pH 2.0 show a different behavior, because BSA overall shape changes as a function of concentration. At 10 mg/mL, BSA is partially unfolded and a strong repulsive protein-protein interaction occurs due to the high amount of exposed charge. At 25 and 50 mg/mL, BSA undergoes some refolding, which likely results in a molten-globule state. This work concludes by confirming that the protein concentration plays an important role on the pH-unfolded BSA state, due to a delicate compromise between interaction forces and crowding effects.

Liquid polymorphism, order-disorder transitions and anomalous behavior: A Monte Carlo study of the Bell-Lavis model for water

The Bell-Lavis model for liquid water is investigated through numerical simulations. The lattice-gas model on a triangular lattice presents orientational states and is known to present a highly bonded low density phase and a loosely bonded high density phase. We show that the model liquid-liquid transition is continuous, in contradiction with mean-field results on the Husimi cactus and from the cluster variational method. We define an order parameter which allows interpretation of the transition as an order-disorder transition of the bond network. Our results indicate that the order-disorder transition is in the Ising universality class. Previous proposal of an Ehrenfest second order transition is discarded. A detailed investigation of anomalous properties has also been undertaken. The line of density maxima in the HDL phase is stabilized by fluctuations, absent in the mean-field solution. (C) 2009 American Institute of Physics. [doi:10.1063/1.3253297]

Noncommutativity due to spin

Using the Berezin-Marinov pseudoclassical formulation of the spin particle we propose a classical model of spin noncommutativity. In the nonrelativistic case, the Poisson brackets between the coordinates are proportional to the spin angular momentum. The quantization of the model leads to the noncommutativity with mixed spatial and spin degrees of freedom. A modified Pauli equation, describing a spin half particle in an external electromagnetic field is obtained. We show that nonlocality caused by the spin noncommutativity depends on the spin of the particle; for spin zero, nonlocality does not appear, for spin half, Delta x Delta y >= theta(2)/2, etc. In the relativistic case the noncommutative Dirac equation was derived. For that we introduce a new star product. The advantage of our model is that in spite of the presence of noncommutativity and nonlocality, it is Lorentz invariant. Also, in the quasiclassical approximation it gives noncommutativity with a nilpotent parameter.

Physical approximations for the nonlinear evolution of perturbations in inhomogeneous dark energy scenarios

The abundance and distribution of collapsed objects such as galaxy clusters will become an important tool to investigate the nature of dark energy and dark matter. Number counts of very massive objects are sensitive not only to the equation of state of dark energy, which parametrizes the smooth component of its pressure, but also to the sound speed of dark energy, which determines the amount of pressure in inhomogeneous and collapsed structures. Since the evolution of these structures must be followed well into the nonlinear regime, and a fully relativistic framework for this regime does not exist yet, we compare two approximate schemes: the widely used spherical collapse model and the pseudo-Newtonian approach. We show that both approximation schemes convey identical equations for the density contrast, when the pressure perturbation of dark energy is parametrized in terms of an effective sound speed. We also make a comparison of these approximate approaches to general relativity in the linearized regime, which lends some support to the approximations.

Kinetic energy sum spectra in nonmesonic weak decay of hypernuclei

We evaluate the coincidence spectra in the nonmesonic weak decay (NMWD) Lambda N -> nN of Lambda hypernuclei (4)(Lambda)He, (5)(Lambda)He, (12)(Lambda)C, (16)(Lambda)O, and (28)(Lambda)Si, as a function of the sum of kinetic energies E(nN)=E(n)+E(N) for N=n,p. The strangeness-changing transition potential is described by the one-meson-exchange model, with commonly used parametrization. Two versions of the independent-particle shell model (IPSM) are employed to account for the nuclear structure of the final residual nuclei. They are as follows: (a) IPSM-a, where no correlation, except for the Pauli principle, is taken into account and (b) IPSM-b, where the highly excited hole states are considered to be quasistationary and are described by Breit-Wigner distributions, whose widths are estimated from the experimental data. All np and nn spectra exhibit a series of peaks in the energy interval 110 MeV < E(nN)< 170 MeV, one for each occupied shell-model state. Within the IPSM-a, and because of the recoil effect, each peak covers an energy interval proportional to A(-1) , going from congruent to 4 MeV for (28)(Lambda)Si to congruent to 40 MeV for (4)(Lambda)He. Such a description could be pretty fair for the light (4)(Lambda)He and (5)(Lambda)He hypernuclei. For the remaining, heavier, hypernuclei it is very important, however, to consider as well the spreading in strength of the deep-hole states and bring into play the IPSM-b approach. Notwithstanding the nuclear model that is employed the results depend only very weakly on the details of the dynamics involved in the decay process proper. We propose that the IPSM is the appropriate lowest-order approximation for the theoretical calculations of the of kinetic energy sum spectra in the NMWD. It is in comparison to this picture that one should appraise the effects of the final-state interactions and of the two-nucleon-induced decay mode.

Hydro-kinetic approach to relativistic heavy ion collisions

We develop a combined hydro-kinetic approach which incorporates a hydrodynamical expansion of the systems formed in A + A collisions and their dynamical decoupling described by escape probabilities. The method corresponds to a generalized relaxation time (tau(rel)) approximation for the Boltzmann equation applied to inhomogeneous expanding systems; at small tau(rel) it also allows one to catch the viscous effects in hadronic component-hadron-resonance gas. We demonstrate how the approximation of sudden freeze-out can be obtained within this dynamical picture of continuous emission and find that hypersurfaces, corresponding to a sharp freeze-out limit, are momentum dependent. The pion m(T) spectra are computed in the developed hydro-kinetic model, and compared with those obtained from ideal hydrodynamics with the Cooper-Frye isothermal prescription. Our results indicate that there does not exist a universal freeze-out temperature for pions with different momenta, and support an earlier decoupling of higher p(T) particles. By performing numerical simulations for various initial conditions and equations of state we identify several characteristic features of the bulk QCD matter evolution preferred in view of the current analysis of heavy ion collisions at RHIC energies.

Optical third-harmonic spectroscopy of the magnetic semiconductor EuTe

EuTe possesses the centrosymmetric crystal structure m3m of rocksalt type in which the second-harmonic generation is forbidden in electric dipole approximation but the third-harmonic generation (THG) is allowed. We studied the THG spectra of this material and observed several resonances in the vicinity of the band gap at 2.2-2.5 eV and at higher energies up to 4 eV, which are related to four-photon THG processes. The observed resonances are assigned to specific combinations of electronic transitions between the ground 4f(7) state at the top of the valence band and excited 4f(6)5d(1) states of Eu(2+) ions, which form the lowest energy conduction band. Temperature, magnetic field, and rotational anisotropy studies allowed us to distinguish crystallographic and magnetic-field-induced contributions to the THG. A strong modification of THG intensity for the 2.4 eV band and suppression of the THG for the 3.15 eV band was observed in applied magnetic field. Two main features of the THG spectra were assigned to 5d(t(2g)) and 5d(e(g)) subbands at 2.4 eV and 3.15 eV, respectively. A microscopic quantum-mechanical model of the THG response was developed and its conclusions are in qualitative agreement with the experimental results.

Systematic investigation of a family of gradient-dependent functionals for solids

Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. On a test set of 60 solids, we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.

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.

Theory of near-gap second harmonic generation in centrosymmetric magnetic semiconductors: Europium chalcogenides

Second harmonic generation is strictly forbidden in centrosymmetric materials, within the electric dipole approximation. Recently, it was found that the centrosymmetric magnetic semiconductors EuTe and EuSe can generate near-gap second harmonics, if the system is submitted to an external magnetic field. Here, a theoretical model is presented, which well describes the observed phenomena. The model shows that second harmonic generation becomes efficient when the magnetic dipole oscillations between the band-edge excited states of the system, induced by the excitation light, enter the in-phase regime, which can be achieved by applying a magnetic field to the material.

Accurate prediction of the Si/SiO(2) interface band offset using the self-consistent ab initio DFT/LDA-1/2 method

We use the density functional theory/local-density approximation (DFT/LDA)-1/2 method [L. G. Ferreira , Phys. Rev. B 78, 125116 (2008)], which attempts to fix the electron self-energy deficiency of DFT/LDA by half-ionizing the whole Bloch band of the crystal, to calculate the band offsets of two Si/SiO(2) interface models. Our results are similar to those obtained with a ""state-of-the-art"" GW approach [R. Shaltaf , Phys. Rev. Lett. 100, 186401 (2008)], with the advantage of being as computationally inexpensive as the usual DFT/LDA. Our band gap and band offset predictions are in excellent agreement with experiments.

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.

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.