Average ground-state energy of finite Fermi systems

Semiclassical theories such as the Thomas-Fermi and Wigner-Kirkwood methods give a good description of the smooth average part of the total energy of a Fermi gas in some external potential when the chemical potential is varied. However, in systems with a fixed number of particles N, these methods overbind the actual average of the quantum energy as N is varied. We describe a theory that accounts for this effect. Numerical illustrations are discussed for fermions trapped in a harmonic oscillator potential and in a hard-wall cavity, and for self-consistent calculations of atomic nuclei. In the latter case, the influence of deformations on the average behavior of the energy is also considered.

Entropy of a correlated system of nucleons

Realistic nucleon-nucleon interactions induce correlations to the nuclear many-body system, which lead to a fragmentation of the single-particle strength over a wide range of energies and momenta. We address the question of how this fragmentation affects the thermodynamical properties of nuclear matter. In particular, we show that the entropy can be computed with the help of a spectral function, which can be evaluated in terms of the self-energy obtained in the self-consistent Green's function approach. Results for the density and temperature dependences of the entropy per particle for symmetric nuclear matter are presented and compared to the results of lowest order finite-temperature Brueckner-Hartree-Fock calculations. The effects of correlations on the calculated entropy are small, if the appropriate quasiparticle approximation is used. The results demonstrate the thermodynamical consistency of the self-consistent T-matrix approximation for the evaluation of the Green's functions.

Liquid-gas phase transition in nuclear matter from realistic many-body approaches

The existence of a liquid-gas phase transition for hot nuclear systems at subsaturation densities is a well-established prediction of finite-temperature nuclear many-body theory. In this paper, we discuss for the first time the properties of such a phase transition for homogeneous nuclear matter within the self-consistent Green's function approach. We find a substantial decrease of the critical temperature with respect to the Brueckner-Hartree-Fock approximation. Even within the same approximation, the use of two different realistic nucleon-nucleon interactions gives rise to large differences in the properties of the critical point.

Dielectric formulation for Many-Particle Systems with General Interactions

An extension of the self-consistent field approach formulation by Cohen in the preceding paper is proposed in order to include the most general kind of two-body interactions, i.e., interactions depending on position, momenta, spin, isotopic spin, etc. The dielectric function is replaced by a dielectric matrix. The evaluation of the energies involves the computation of a matrix inversion and trace.

(K)over-bar* mesons in dense matter

We study the properties of (K) over bar* mesons in nuclear matter using a unitary approach in coupled channels within the framework of the local hidden gauge formalism and incorporating the (K) over bar pi decay channel in matter. The in-medium (K) over bar *N interaction accounts for Pauli blocking effects and incorporates the (K) over bar* self-energy in a self-consistent manner. We also obtain the (K) over bar* (off-shell) spectral function and analyze its behavior at finite density and momentum. At a normal nuclear matter density, the (K) over bar* meson feels a moderately attractive potential, while the (K) over bar* width becomes five times larger than in free space. We estimate the transparency ratio of the gamma A -> K+K*(-) A` reaction, which we propose as a feasible scenario at the present facilities to detect changes in the properties of the (K) over bar* meson in nuclear medium.

Optical-model potential for electron and positron elastic scattering by atoms

An optical-model potential for systematic calculations of elastic scattering of electrons and positrons by atoms and positive ions is proposed. The electrostatic interaction is determined from the Dirac-Hartree-Fock self-consistent atomic electron density. In the case of electron projectiles, the exchange interaction is described by means of the local-approximation of Furness and McCarthy. The correlation-polarization potential is obtained by combining the correlation potential derived from the local density approximation with a long-range polarization interaction, which is represented by means of a Buckingham potential with an empirical energy-dependent cutoff parameter. The absorption potential is obtained from the local-density approximation, using the Born-Ochkur approximation and the Lindhard dielectric function to describe the binary collisions with a free-electron gas. The strength of the absorption potential is adjusted by means of an empirical parameter, which has been determined by fitting available absolute elastic differential cross-section data for noble gases and mercury. The Dirac partial-wave analysis with this optical-model potential provides a realistic description of elastic scattering of electrons and positrons with energies in the range from ~100 eV up to ~5 keV. At higher energies, correlation-polarization and absorption corrections are small and the usual static-exchange approximation is sufficiently accurate for most practical purposes.

Stability of de Sitter spacetime under isotropic perturbations in semiclassical gravity

A spatially flat Robertson-Walker spacetime driven by a cosmological constant is nonconformally coupled to a massless scalar field. The equations of semiclassical gravity are explicitly solved for this case, and a self-consistent de Sitter solution associated with the Bunch-Davies vacuum state is found (the effect of the quantum field is to shift slightly the effective cosmological constant). Furthermore, it is shown that the corrected de Sitter spacetime is stable under spatially isotropic perturbations of the metric and the quantum state. These results are independent of the free renormalization parameters.

Adsorption kinetics in the presence of external fields

In this article we present a detailed analysis of the kinetics of a class of sequential adsorption models that take into account the effect of externally applied fields (as an electric field, or a shear rate) on the adsorption. The excluded volume interactions related to the finite size of the adsorbing particles are modified by the external fields. As a result, new adsorption mechanisms appear with respect to the ones used to describe the kinetics in a quiescent fluid. In particular, if the adsorbing particles are allowed to roll over preadsorbed ones, adsorption becomes non local even in the simplest geometry. An exact analytic theory cannot be developed, but we introduce a self-consistent theory that turns out to agree with the simulation results over all the range of the parameters.

Adsorption kinetics in the presence of external fields

In this article we present a detailed analysis of the kinetics of a class of sequential adsorption models that take into account the effect of externally applied fields (as an electric field, or a shear rate) on the adsorption. The excluded volume interactions related to the finite size of the adsorbing particles are modified by the external fields. As a result, new adsorption mechanisms appear with respect to the ones used to describe the kinetics in a quiescent fluid. In particular, if the adsorbing particles are allowed to roll over preadsorbed ones, adsorption becomes non local even in the simplest geometry. An exact analytic theory cannot be developed, but we introduce a self-consistent theory that turns out to agree with the simulation results over all the range of the parameters.

Pionic branches in neutron matter and the problem of pion condensation

The pion spectrum for charged and neutral pions is investigated in pure neutron matter, by letting the pions interact with a neutron Fermi sea in a self-consistent scheme that renormalizes simultaneously the mesons, considered the source of the interaction, and the nucleons. The possibility of obtaining different kinds of pion condensates is investigated with the result that they cannot be reached even for values of the spin-spin correlation parameter, g', far below the range commonly accepted.

Neutral atoms in ionic lattices: Excited states of KCl:Ag(0)

The optical-absorption spectrum of a cationic Ag0 atom in a KCl crystal has been studied theoretically by means of a series of cluster models of increasing size. Excitation energies have been determined by means of a multiconfigurational self-consistent field procedure followed by a second-order perturbation correlation treatment. Moreover results obtained within the density-functional framework are also reported. The calculations confirm the assignment of bands I and IV to transitions of the Ag-5s electron into delocalized states with mainly K-4s,4p character. Bands II and III have been assigned to internal transitions on the Ag atom, which correspond to the atomic Ag-4d to Ag-5s transition. We also determine the lowest charge transfer (CT) excitation energy and confirm the assignment of band VI to such a transition. The study of the variation of the CT excitation energy with the Ag-Cl distance R gives additional support to a large displacement of the Cl ions due to the presence of the Ag0 impurity. Moreover, from the present results, it is predicted that on passing to NaCl:Ag0 the CT onset would be out of the optical range while the 5s-5p transition would undergo a redshift of 0.3 eV. These conclusions, which underline the different character of involved orbitals, are consistent with experimental findings. The existence of a CT transition in the optical range for an atom inside an ionic host is explained by a simple model, which also accounts for the differences with the more common 3d systems. The present study sheds also some light on the R dependence of the s2-sp transitions due to s2 ions like Tl+.

First principles LDA+U and GGA+U study of cerium oxides: Dependence on the effective U-parameter

The electronic structure and properties of cerium oxides (CeO2 and Ce2O3) have been studied in the framework of the LDA+U and GGA(PW91)+U implementations of density functional theory. The dependence of selected observables of these materials on the effective U parameter has been investigated in detail. The examined properties include lattice constants, bulk moduli, density of states, and formation energies of CeO2 and Ce2O3. For CeO2, the LDA+U results are in better agreement with experiment than the GGA+U results whereas for the computationally more demanding Ce2O3 both approaches give comparable accuracy. Furthermore, as expected, Ce2O3 is much more sensitive to the choice of the U value. Generally, the PW91 functional provides an optimal agreement with experiment at lower U energies than LDA does. In order to achieve a balanced description of both kinds of materials, and also of nonstoichiometric CeO2¿x phases, an appropriate choice of U is suggested for LDA+U and GGA+U schemes. Nevertheless, an optimum value appears to be property dependent, especially for Ce2O3. Optimum U values are found to be, in general, larger than values determined previously in a self-consistent way.

Thomas-Fermi theory for atomic nuclei revisited

The recently developed semiclassical variational Wigner-Kirkwood (VWK) approach is applied to finite nuclei using external potentials and self-consistent mean fields derived from Skyrme inter-actions and from relativistic mean field theory. VWK consist s of the Thomas-Fermi part plus a pure, perturbative h 2 correction. In external potentials, VWK passes through the average of the quantal values of the accumulated level density and total en energy as a function of the Fermi energy. However, there is a problem of overbinding when the energy per particle is displayed as a function of the particle number. The situation is analyzed comparing spherical and deformed harmonic oscillator potentials. In the self-consistent case, we show for Skyrme forces that VWK binding energies are very close to those obtained from extended Thomas-Fermi functionals of h 4 order, pointing to the rapid convergence of the VWK theory. This satisfying result, however, does not cure the overbinding problem, i.e., the semiclassical energies show more binding than they should. This feature is more pronounced in the case of Skyrme forces than with the relativistic mean field approach. However, even in the latter case the shell correction energy for e.g.208 Pb turns out to be only ∼ −6 MeV what is about a factor two or three off the generally accepted value. As an adhoc remedy, increasing the kinetic energy by 2.5%, leads to shell correction energies well acceptable throughout the periodic table. The general importance of the present studies for other finite Fermi systems, self-bound or in external potentials, is pointed out.

Analysis of bulk and surface contributions in the neutron skin of nuclei

The neutron skin thickness of nuclei is a sensitive probe of the nuclear symmetry energy and has multiple implications for nuclear and astrophysical studies. However, precision measurements of this observable are difficult to obtain. The analysis of the experimental data may imply some assumptions about the bulk or surface nature of the formation of the neutron skin. Here we study the bulk or surface character of neutron skins of nuclei following from calculations with Gogny, Skyrme, and covariant nuclear mean-field interactions. These interactions are successful in describing nuclear charge radii and binding energies but predict different values for neutron skins. We perform the study by fitting two-parameter Fermi distributions to the calculated self-consistent neutron and proton densities. We note that the equivalent sharp radius is a more suitable reference quantity than the half-density radius parameter of the Fermi distributions to discern between the bulk and surface contributions in neutron skins. We present calculations for nuclei in the stability valley and for the isotopic chains of Sn and Pb.

Toward on-the-fly quantum mechanical/molecular mechanical (QM/MM) docking: development and benchmark of a scoring function.

We address the challenges of treating polarization and covalent interactions in docking by developing a hybrid quantum mechanical/molecular mechanical (QM/MM) scoring function based on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. To benchmark this scoring function within the EADock DSS docking algorithm, we created a publicly available dataset of high-quality X-ray structures of zinc metalloproteins ( http://www.molecular-modelling.ch/resources.php ). For zinc-bound ligands (226 complexes), the QM/MM scoring yielded a substantially improved success rate compared to the classical scoring function (77.0% vs 61.5%), while, for allosteric ligands (55 complexes), the success rate remained constant (49.1%). The QM/MM scoring significantly improved the detection of correct zinc-binding geometries and improved the docking success rate by more than 20% for several important drug targets. The performance of both the classical and the QM/MM scoring functions compare favorably to the performance of AutoDock4, AutoDock4Zn, and AutoDock Vina.