Monte Carlo investigation of the critical behavior of Stavskaya's probabilistic cellular automaton

Stavskaya's model is a one-dimensional probabilistic cellular automaton (PCA) introduced in the end of the 1960s as an example of a model displaying a nonequilibrium phase transition. Although its absorbing state phase transition is well understood nowadays, the model never received a full numerical treatment to investigate its critical behavior. In this Brief Report we characterize the critical behavior of Stavskaya's PCA by means of Monte Carlo simulations and finite-size scaling analysis. The critical exponents of the model are calculated and indicate that its phase transition belongs to the directed percolation universality class of critical behavior, as would be expected on the basis of the directed percolation conjecture. We also explicitly establish the relationship of the model with the Domany-Kinzel PCA on its directed site percolation line, a connection that seems to have gone unnoticed in the literature so far.

Renyi entropy and parity oscillations of anisotropic spin-s Heisenberg chains in a magnetic field

Using the density matrix renormalization group, we investigate the Renyi entropy of the anisotropic spin-s Heisenberg chains in a z-magnetic field. We considered the half-odd-integer spin-s chains, with s = 1/2, 3/2, and 5/2, and periodic and open boundary conditions. In the case of the spin-1/2 chain we were able to obtain accurate estimates of the new parity exponents p(alpha)((p)) and p(alpha)((o)) that gives the power-law decay of the oscillations of the alpha-Renyi entropy for periodic and open boundary conditions, respectively. We confirm the relations of these exponents with the Luttinger parameter K, as proposed by Calabrese et al. [Phys. Rev. Lett. 104, 095701 (2010)]. Moreover, the predicted periodicity of the oscillating term was also observed for some nonzero values of the magnetization m. We show that for s > 1/2 the amplitudes of the oscillations are quite small and get accurate estimates of p(alpha)((p)) and p(alpha)((o)) become a challenge. Although our estimates of the new universal exponents p(alpha)((p)) and p(alpha)((o)) for the spin-3/2 chain are not so accurate, they are consistent with the theoretical predictions.

Precise determination of quantum critical points by the violation of the entropic area law

Finite-size scaling analysis turns out to be a powerful tool to calculate the phase diagram as well as the critical properties of two-dimensional classical statistical mechanics models and quantum Hamiltonians in one dimension. The most used method to locate quantum critical points is the so-called crossing method, where the estimates are obtained by comparing the mass gaps of two distinct lattice sizes. The success of this method is due to its simplicity and the ability to provide accurate results even considering relatively small lattice sizes. In this paper, we introduce an estimator that locates quantum critical points by exploring the known distinct behavior of the entanglement entropy in critical and noncritical systems. As a benchmark test, we use this new estimator to locate the critical point of the quantum Ising chain and the critical line of the spin-1 Blume-Capel quantum chain. The tricritical point of this last model is also obtained. Comparison with the standard crossing method is also presented. The method we propose is simple to implement in practice, particularly in density matrix renormalization group calculations, and provides us, like the crossing method, amazingly accurate results for quite small lattice sizes. Our applications show that the proposed method has several advantages, as compared with the standard crossing method, and we believe it will become popular in future numerical studies.

Nonadiabatic Coherent Evolution of Two-Level Systems under Spontaneous Decay

In this Letter we extend current perspectives in engineering reservoirs by producing a time-dependent master equation leading to a nonstationary superposition equilibrium state that can be nonadiabatically controlled by the system-reservoir parameters. Working with an ion trapped inside a nonideal cavity, we first engineer effective interactions, which allow us to achieve two classes of decoherence-free evolution of superpositions of the ground and excited ionic levels: those with a time-dependent azimuthal or polar angle. As an application, we generalize the purpose of an earlier study [Phys. Rev. Lett. 96, 150403 (2006)], showing how to observe the geometric phases acquired by the protected nonstationary states even under nonadiabatic evolution.

Entanglement in spatially inhomogeneous many-fermion systems

We investigate entanglement of strongly interacting fermions in spatially inhomogeneous environments. To quantify entanglement in the presence of spatial inhomogeneity, we propose a local-density approximation (LDA) to the entanglement entropy, and a nested LDA scheme to evaluate the entanglement entropy on inhomogeneous density profiles. These ideas are applied to models of electrons in superlattice structures with different modulation patterns, electrons in a metallic wire in the presence of impurities, and phase-separated states in harmonically confined many-fermion systems, such as electrons in quantum dots and atoms in optical traps. We find that the entanglement entropy of inhomogeneous systems is strikingly different from that of homogeneous systems.

Effect of spatial inhomogeneity on the mapping between strongly interacting fermions and weakly interacting spins

A combined analytical and numerical study is performed of the mapping between strongly interacting fermions and weakly interacting spins, in the framework of the Hubbard, t-J, and Heisenberg models. While for spatially homogeneous models in the thermodynamic limit the mapping is thoroughly understood, we here focus on aspects that become relevant in spatially inhomogeneous situations, such as the effect of boundaries, impurities, superlattices, and interfaces. We consider parameter regimes that are relevant for traditional applications of these models, such as electrons in cuprates and manganites, and for more recent applications to atoms in optical lattices. The rate of the mapping as a function of the interaction strength is determined from the Bethe-Ansatz for infinite systems and from numerical diagonalization for finite systems. We show analytically that if translational symmetry is broken through the presence of impurities, the mapping persists and is, in a certain sense, as local as possible, provided the spin-spin interaction between two sites of the Heisenberg model is calculated from the harmonic mean of the onsite Coulomb interaction on adjacent sites of the Hubbard model. Numerical calculations corroborate these findings also in interfaces and superlattices, where analytical calculations are more complicated.

Nonclassical correlation in NMR quadrupolar systems

The existence of quantum correlation (as revealed by quantum discord), other than entanglement and its role in quantum-information processing (QIP), is a current subject for discussion. In particular, it has been suggested that this nonclassical correlation may provide computational speedup for some quantum algorithms. In this regard, bulk nuclear magnetic resonance (NMR) has been successfully used as a test bench for many QIP implementations, although it has also been continuously criticized for not presenting entanglement in most of the systems used so far. In this paper, we report a theoretical and experimental study on the dynamics of quantum and classical correlations in an NMR quadrupolar system. We present a method for computing the correlations from experimental NMR deviation-density matrices and show that, given the action of the nuclear-spin environment, the relaxation produces a monotonic time decay in the correlations. Although the experimental realizations were performed in a specific quadrupolar system, the main results presented here can be applied to whichever system uses a deviation-density matrix formalism.

Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is known about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they underestimate charge gaps.

Reverse engineering in many-body quantum physics: Correspondence between many-body systems and effective single-particle equations

The mapping, exact or approximate, of a many-body problem onto an effective single-body problem is one of the most widely used conceptual and computational tools of physics. Here, we propose and investigate the inverse map of effective approximate single-particle equations onto the corresponding many-particle system. This approach allows us to understand which interacting system a given single-particle approximation is actually describing, and how far this is from the original physical many-body system. We illustrate the resulting reverse engineering process by means of the Kohn-Sham equations of density-functional theory. In this application, our procedure sheds light on the nonlocality of the density-potential mapping of density-functional theory, and on the self-interaction error inherent in approximate density functionals.

POLLING SYSTEMS WITH PARAMETER REGENERATION, THE GENERAL CASE

We consider a polling model with multiple stations, each with Poisson arrivals and a queue of infinite capacity. The service regime is exhaustive and there is Jacksonian feedback of served customers. What is new here is that when the server comes to a station it chooses the service rate and the feedback parameters at random; these remain valid during the whole stay of the server at that station. We give criteria for recurrence, transience and existence of the sth moment of the return time to the empty state for this model. This paper generalizes the model, when only two stations accept arriving jobs, which was considered in [Ann. Appl. Probab. 17 (2007) 1447-1473]. Our results are stated in terms of Lyapunov exponents for random matrices. From the recurrence criteria it can be seen that the polling model with parameter regeneration can exhibit the unusual phenomenon of null recurrence over a thick region of parameter space.

Chromic materials for responsive surface-enhanced resonance Raman scattering systems: a nanometric pH sensor

The use of chromic materials for responsive surface-enhanced resonance Raman scattering (SERRS) based nanosensors is reported. The potential of nano-chromic SERRS is demonstrated with the use of the halochrome methyl yellow to fabricate an ultrasensitive pH optical sensor. Some of the challenges of the incorporation of chromic materials with metal nanostructures are addressed through the use of computational calculations and a comparison to measured SERRS and surface-enhanced Raman scattering (SERS) spectra is presented. A strong correlation between the measured SERRS and the medium's proton concentration is demonstrated for the pH range 2-6. The high sensitivity achieved by the use of resonance Raman conditions is shown through responsive SERRS measurements from only femtolitres of volume and with the concentration of the reporting molecules approaching the single molecule regime.

Effect of sugarcane harvesting systems on soil carbon stocks in Brazil: an examination of existing data

Agricultural management practices that promote net carbon (C) accumulation in the soil have been considered as an important potential mitigation option to combat global warming. The change in the sugarcane harvesting system, to one which incorporates C into the soil from crop residues, is the focus of this work. The main objective was to assess and discuss the changes in soil organic C stocks caused by the conversion of burnt to unburnt sugarcane harvesting systems in Brazil, when considering the main soils and climates associated with this crop. For this purpose, a dataset was obtained from a literature review of soils under sugarcane in Brazil. Although not necessarily from experimental studies, only paired comparisons were examined, and for each site the dominant soil type, topography and climate were similar. The results show a mean annual C accumulation rate of 1.5 Mg ha-1 year-1 for the surface to 30-cm depth (0.73 and 2.04 Mg ha-1 year-1 for sandy and clay soils, respectively) caused by the conversion from a burnt to an unburnt sugarcane harvesting system. The findings suggest that soil should be included in future studies related to life cycle assessment and C footprint of Brazilian sugarcane ethanol.

Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets

Outgassing of carbon dioxide (CO(2)) from rivers and streams to the atmosphere is a major loss term in the coupled terrestrial-aquatic carbon cycle of major low-gradient river systems (the term ""river system"" encompasses the rivers and streams of all sizes that compose the drainage network in a river basin). However, the magnitude and controls on this important carbon flux are not well quantified. We measured carbon dioxide flux rates (F(CO2)), gas transfer velocity (k), and partial pressures (p(CO2)) in rivers and streams of the Amazon and Mekong river systems in South America and Southeast Asia, respectively. F(CO2) and k values were significantly higher in small rivers and streams (channels <100 m wide) than in large rivers (channels >100 m wide). Small rivers and streams also had substantially higher variability in k values than large rivers. Observed F(CO2) and k values suggest that previous estimates of basinwide CO(2) evasion from tropical rivers and wetlands have been conservative and are likely to be revised upward substantially in the future. Data from the present study combined with data compiled from the literature collectively suggest that the physical control of gas exchange velocities and fluxes in low-gradient river systems makes a transition from the dominance of wind control at the largest spatial scales (in estuaries and river mainstems) toward increasing importance of water current velocity and depth at progressively smaller channel dimensions upstream. These results highlight the importance of incorporating scale-appropriate k values into basinwide models of whole ecosystem carbon balance.

Soil carbon stocks under no-tillage mulch-based cropping systems in the Brazilian Cerrado: An on-farm synchronic assessment

No-tillage mulch-based (NTM) cropping systems have been widely adopted by farmers in the Brazilian savanna region (Cerrado biome). We hypothesized that this new type of management should have a profound impact on soil organic carbon (SOC) at regional scale and consequently on climate change mitigation. The objective of this study was thus to quantify the SOC storage potential of NTM in the oxisols of the Cerrado using a synchronic approach that is based on a chronosequence of fields of different years under NTM. The study consisted of three phases: (1) a farm/cropping system survey to identify the main types of NTM systems to be chosen for the chronosequence; (2) a field survey to identify a homogeneous set of situations for the chronosequence and (3) the characterization of the chronosequence to assess the SOC storage potential. The main NTM system practiced by farmers is an annual succession of soybean (Glycine max)or maize (Zea mays) with another cereal crop. This cropping system covers 54% of the total cultivated area in the region. At the regional level, soil organic C concentrations from NTM fields were closely correlated with clay + silt content of the soil (r(2) = 0.64). No significant correlation was observed (r(2) = 0.07), however, between these two variables when we only considered the fields with a clay + silt content in the 500-700 g kg(-1) range. The final chronosequence of NTM fields was therefore based on a subsample of eight fields, within this textural range. The SOC stocks in the 0-30 cm topsoil layer of these selected fields varied between 4.2 and 6.7 kg C m(-2) and increased on average (r(2) = 0.97) with 0.19 kg C m(-2) year(-1). After 12 years of NTM management, SOC stocks were no longer significantly different from the stocks under natural Cerrado vegetation (p < 0.05), whereas a 23-year-old conventionally tilled and cropped field showed SOC stocks that were about 30% below this level. Confirming our hypotheses, this study clearly illustrated the high potential of NTM systems in increasing SOC storage under tropical conditions, and how a synchronic approach may be used to assess efficiently such modification on farmers` fields, identifying and excluding non desirable sources of heterogeneity (management, soils and climate). (C) 2010 Elsevier B.V. All rights reserved.

Net greenhouse gas fluxes in Brazilian ethanol production systems

Biofuels are both a promising solution to global warming mitigation and a potential contributor to the problem. Several life cycle assessments of bioethanol have been conducted to address these questions. We performed a synthesis of the available data on Brazilian ethanol production focusing on greenhouse gas (GHG) emissions and carbon (C) sinks in the agricultural and industrial phases. Emissions of carbon dioxide (CO(2)) from fossil fuels, methane (CH(4)) and nitrous oxide (N(2)O) from sources commonly included in C footprints, such as fossil fuel usage, biomass burning, nitrogen fertilizer application, liming and litter decomposition were accounted for. In addition, black carbon (BC) emissions from burning biomass and soil C sequestration were included in the balance. Most of the annual emissions per hectare are in the agricultural phase, both in the burned system (2209 out of a total of 2398 kg C(eq)), and in the unburned system (559 out of 748 kg C(eq)). Although nitrogen fertilizer emissions are large, 111 kg C(eq) ha-1 yr-1, the largest single source of emissions is biomass burning in the manual harvest system, with a large amount of both GHG (196 kg C(eq) ha-1 yr-1). and BC (1536 kg C(eq) ha-1 yr-1). Besides avoiding emissions from biomass burning, harvesting sugarcane mechanically without burning tends to increase soil C stocks, providing a C sink of 1500 kg C ha-1 yr-1 in the 30 cm layer. The data show a C output: input ratio of 1.4 for ethanol produced under the conventionally burned and manual harvest compared with 6.5 for the mechanized harvest without burning, signifying the importance of conservation agricultural systems in bioethanol feedstock production.