A simple topological model with continuous phase transition

In the area of topological and geometric treatment of phase transitions and symmetry breaking in Hamiltonian systems, in a recent paper some general sufficient conditions for these phenomena in $\mathbb{Z}_2$-symmetric systems (i.e. invariant under reflection of coordinates) have been found out. In this paper we present a simple topological model satisfying the above conditions hoping to enlighten the mechanism which causes this phenomenon in more general physical models. The symmetry breaking is testified by a continuous magnetization with a nonanalytic point in correspondence of a critical temperature which divides the broken symmetry phase from the unbroken one. A particularity with respect to the common pictures of a phase transition is that the nonanalyticity of the magnetization is not accompanied by a nonanalytic behavior of the free energy.Comment: 17 pages, 7 figure

Topological conditions for discrete symmetry breaking and phase transitions

In the framework of a recently proposed topological approach to phase transitions, some sufficient conditions ensuring the presence of the spontaneous breaking of a Z_2 symmetry and of a symmetry-breaking phase transition are introduced and discussed. A very simple model, which we refer to as the hypercubic model, is introduced and solved. The main purpose of this model is that of illustrating the content of the sufficient conditions, but it is interesting also in itself due to its simplicity. Then some mean-field models already known in the literature are discussed in the light of the sufficient conditions introduced here

Uncertainty in the determination of soil hydraulic parameters and its influence on the performance of two hydrological models of different complexity

Data of soil hydraulic properties forms often a limiting factor in unsaturated zone modelling, especially at the larger scales. Investigations for the hydraulic characterization of soils are time-consuming and costly, and the accuracy of the results obtained by the different methodologies is still debated. However, we may wonder how the uncertainty in soil hydraulic parameters relates to the uncertainty of the selected modelling approach. We performed an intensive monitoring study during the cropping season of a 10 ha maize field in Northern Italy. The data were used to: i) compare different methods for determining soil hydraulic parameters and ii) evaluate the effect of the uncertainty in these parameters on different variables (i.e. evapotranspiration, average water content in the root zone, flux at the bottom boundary of the root zone) simulated by two hydrological models of different complexity: SWAP, a widely used model of soil moisture dynamics in unsaturated soils based on Richards equation, and ALHyMUS, a conceptual model of the same dynamics based on a reservoir cascade scheme. We employed five direct and indirect methods to determine soil hydraulic parameters for each horizon of the experimental profile. Two methods were based on a parameter optimization of: a) laboratory measured retention and hydraulic conductivity data and b) field measured retention and hydraulic conductivity data. The remaining three methods were based on the application of widely used Pedo-Transfer Functions: c) Rawls and Brakensiek, d) HYPRES, and e) ROSETTA. Simulations were performed using meteorological, irrigation and crop data measured at the experimental site during the period June – October 2006. Results showed a wide range of soil hydraulic parameter values generated with the different methods, especially for the saturated hydraulic conductivity Ksat and the shape parameter a of the van Genuchten curve. This is reflected in a variability of the modeling results which is, as expected, different for each model and each variable analysed. The variability of the simulated water content in the root zone and of the bottom flux for different soil hydraulic parameter sets is found to be often larger than the difference between modeling results of the two models using the same soil hydraulic parameter set. Also we found that a good agreement in simulated soil moisture patterns may occur even if evapotranspiration and percolation fluxes are significantly different. Therefore multiple output variables should be considered to test the performances of methods and model

Accurate quadratic-response approximation for the self-consistent pseudopotential of semiconductor nanostructures

Quadratic-response theory is shown to provide a conceptually simple but accurate approximation for the self-consistent one-electron potential of semiconductor nanostructures. Numerical examples are presented for GaAs/AlAs and InGaAs/InP (001) superlattices using the local-density approximation to density-functional theory and norm-conserving pseudopotentials without spin-orbit coupling. When the reference crystal is chosen to be the virtual-crystal average of the two bulk constituents, the absolute error in the quadratic-response potential for Gamma(15) valence electrons is about 2 meV for GaAs/AlAs and 5 meV for InGaAs/InP. Low-order multipole expansions of the electron density and potential response are shown to be accurate throughout a small neighborhood of each reciprocal lattice vector, thus providing a further simplification that is confirmed to be valid for slowly varying envelope functions. Although the linear response is about an order of magnitude larger than the quadratic response, the quadratic terms are important both quantitatively (if an accuracy of better than a few tens of meV is desired) and qualitatively (due to their different symmetry and long-range dipole effects).Comment: 16 pages, 20 figures; v2: new section on limitations of theor

Dynamical-charge neutrality at a crystal surface

For both molecules and periodic solids, the ionic dynamical charge tensors which govern the infrared activity are known to obey a dynamical neutrality condition. This condition enforces their sum to vanish (over the whole finite system, or over the crystal cell, respectively). We extend this sum rule to the non trivial case of the surface of a semiinfinite solid and show that, in the case of a polar surface of an insulator, the surface ions cannot have the same dynamical charges as in the bulk. The sum rule is demonstrated through calculations for the Si-terminated SiC(001) surface.Comment: 4 pages, latex file, 1 postscript figure automatically include

Antagonism between salicylate and the cAMP signal controls yeast cell survival and growth recovery from quiescence

Aspirin and its main metabolite salicylate are promising molecules in preventing cancer and metabolic diseases. S. cerevisiae cells have been used to study some of their effects: (i) salicylate induces the reversible inhibition of both glucose transport and the biosyntheses of glucose-derived sugar phosphates, (ii) Aspirin/salicylate causes apoptosis associated with superoxide radical accumulation or early cell necrosis in MnSOD-deficient cells growing in ethanol or in glucose, respectively. So, treatment with (acetyl)-salicylic acid can alter the yeast metabolism and is associated with cell death. We describe here the dramatic effects of salicylate on cellular control of the exit from a quiescence state. The growth recovery of long-term stationary phase cells was strongly inhibited in the presence of salicylate, to a degree proportional to the drug concentration. At high salicylate concentration, growth reactivation was completely repressed and associated with a dramatic loss of cell viability. Strikingly, both of these phenotypes were fully suppressed by increasing the cAMP signal without any variation of the exponential growth rate. Upon nutrient exhaustion, salicylate induced a premature lethal cell cycle arrest in the budded-G2/M phase that cannot be suppressed by PKA activation. We discuss how the dramatic antagonism between cAMP and salicylate could be conserved and impinge common targets in yeast and humans. Targeting quiescence of cancer cells with stem-like properties and their growth recovery from dormancy are major challenges in cancer therapy. If mechanisms underlying cAMP-salicylate antagonism will be defined in our model, this might have significant therapeutic implications

Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion

We describe how to apply the recently developed pole expansion and selected inversion (PEXSI) technique to Kohn-Sham density function theory (DFT) electronic structure calculations that are based on atomic orbital discretization. We give analytic expressions for evaluating the charge density, the total energy, the Helmholtz free energy and the atomic forces (including both the Hellman-Feynman force and the Pulay force) without using the eigenvalues and eigenvectors of the Kohn-Sham Hamiltonian. We also show how to update the chemical potential without using Kohn-Sham eigenvalues. The advantage of using PEXSI is that it has a much lower computational complexity than that associated with the matrix diagonalization procedure. We demonstrate the performance gain by comparing the timing of PEXSI with that of diagonalization on insulating and metallic nanotubes. For these quasi-1D systems, the complexity of PEXSI is linear with respect to the number of atoms. This linear scaling can be observed in our computational experiments when the number of atoms in a nanotube is larger than a few hundreds. Both the wall clock time and the memory requirement of PEXSI is modest. This makes it even possible to perform Kohn-Sham DFT calculations for 10,000-atom nanotubes with a sequential implementation of the selected inversion algorithm. We also perform an accurate geometry optimization calculation on a truncated (8,0) boron-nitride nanotube system containing 1024 atoms. Numerical results indicate that the use of PEXSI does not lead to loss of accuracy required in a practical DFT calculation

Medium polarization isotopic effects on nuclear binding energies

There exist several effective interactions whose parameters are fitted to force mean field predictions to reproduce experimental findings of finite nuclei and calculated properties of infinite nuclear matter. Exploiting this tecnique one can give a good description of nuclear binding energies. We present evidence that further progress can be made by taking into account medium polarization effects associated with surface and pairing vibrations.Comment: 7 pages, 5 figure

Sampling Molecular Conformers in Solution with Quantum Mechanical Accuracy at a Nearly Molecular-Mechanics Cost

We introduce a method to evaluate the relative populations of different conformers of molecular species in solution, aiming at quantum mechanical accuracy, while keeping the computational cost at a nearly molecular-mechanics level. This goal is achieved by combining long classical molecular-dynamics simulations to sample the free-energy landscape of the system, advanced clustering techniques to identify the most relevant conformers, and thermodynamic perturbation theory to correct the resulting populations, using quantum-mechanical energies from density functional theory. A quantitative criterion for assessing the accuracy thus achieved is proposed. The resulting methodology is demonstrated in the specific case of cyanin (cyanidin-3-glucoside) in water solution

Lattice instabilities of PbZrO3/PbTiO3 [1:1] superlattices from first principles

Ab initio phonon calculations for the nonpolar reference structures of the (001), (110), and (111) PbZrO_3/PbTiO_3 [1:1] superlattices are presented. The unstable polar modes in the tetragonal (001) and (110) structures are confined in either the Ti- or the Zr-centered layers and display two-mode behavior, while in the cubic (111) case one-mode behavior is observed. Instabilities with pure oxygen character are observed in all three structures. The implications for the ferroelectric behavior and related properties are discussed.Comment: 12 pages, 2 figures, 7 tables, submitted to PR