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.

Many-Body Effects on the rho(xx) Ringlike Structures in Two-Subband Wells

The longitudinal resistivity rho(xx) of two-dimensional electron gases formed in wells with two subbands displays ringlike structures when plotted in a density-magnetic-field diagram, due to the crossings of spin-split Landau levels (LLs) from distinct subbands. Using spin density functional theory and linear response, we investigate the shape and spin polarization of these structures as a function of temperature and magnetic-field tilt angle. We find that (i) some of the rings ""break'' at sufficiently low temperatures due to a quantum Hall ferromagnetic phase transition, thus exhibiting a high degree of spin polarization (similar to 50%) within, consistent with the NMR data of Zhang et al. [Phys. Rev. Lett. 98, 246802 (2007)], and (ii) for increasing tilting angles the interplay between the anticrossings due to inter-LL couplings and the exchange-correlation effects leads to a collapse of the rings at some critical angle theta(c), in agreement with the data of Guo et al. [Phys. Rev. B 78, 233305 (2008)].

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov-Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics. (C) 2007 Elsevier B.V. All rights reserved.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

QUADRANT-SEARCHING: A NEW TECHNIQUE TO REDUCE THE SEARCH SPACE OF THE INVERSE PROBLEM OF EIT USING BEM TO THE DIRECT PROBLEM

The image reconstruction using the EIT (Electrical Impedance Tomography) technique is a nonlinear and ill-posed inverse problem which demands a powerful direct or iterative method. A typical approach for solving the problem is to minimize an error functional using an iterative method. In this case, an initial solution close enough to the global minimum is mandatory to ensure the convergence to the correct minimum in an appropriate time interval. The aim of this paper is to present a new, simple and low cost technique (quadrant-searching) to reduce the search space and consequently to obtain an initial solution of the inverse problem of EIT. This technique calculates the error functional for four different contrast distributions placing a large prospective inclusion in the four quadrants of the domain. Comparing the four values of the error functional it is possible to get conclusions about the internal electric contrast. For this purpose, initially we performed tests to assess the accuracy of the BEM (Boundary Element Method) when applied to the direct problem of the EIT and to verify the behavior of error functional surface in the search space. Finally, numerical tests have been performed to verify the new technique.

(85)Rb Bose-Einstein condensate with tunable interaction: A quantum many body approach

We present a quantum many body approach with van der Waal type of interaction to achieve (85)Rb Bose-Einstein condensate with tunable interaction which has been produced by magnetic field induced Feshbach resonance in the JILA experiment. (C) 2008 Elsevier B.V. All rights reserved.

A Goldstone theorem in thermal relativistic quantum field theory

We prove a Goldstone theorem in thermal relativistic quantum field theory, which relates spontaneous symmetry breaking to the rate of spacelike decay of the two-point function. The critical rate of fall-off coincides with that of the massless free scalar field theory. Related results and open problems are briefly discussed. (C) 2011 American Institute of Physics. [doi:10.1063/1.3526961]

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.

Evolutionary dynamics on rugged fitness landscapes: Exact dynamics and information theoretical aspects

The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the random energy model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find in addition to these three regimes, a ferromagnetic regime where the overlap and the mean-fitness are frozen. In this case, the system relaxes to equilibrium in a finite time. The relevance of our results to information processing aspects of evolution is discussed.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Opinion dynamics of learning agents: does seeking consensus lead to disagreement?

We study opinion dynamics in a population of interacting adaptive agents voting on a set of issues represented by vectors. We consider agents who can classify issues into one of two categories and can arrive at their opinions using an adaptive algorithm. Adaptation comes from learning and the information for the learning process comes from interacting with other neighboring agents and trying to change the internal state in order to concur with their opinions. The change in the internal state is driven by the information contained in the issue and in the opinion of the other agent. We present results in a simple yet rich context where each agent uses a Boolean perceptron to state their opinion. If the update occurs with information asynchronously exchanged among pairs of agents, then the typical case, if the number of issues is kept small, is the evolution into a society torn by the emergence of factions with extreme opposite beliefs. This occurs even when seeking consensus with agents with opposite opinions. If the number of issues is large, the dynamics becomes trapped, the society does not evolve into factions and a distribution of moderate opinions is observed. The synchronous case is technically simpler and is studied by formulating the problem in terms of differential equations that describe the evolution of order parameters that measure the consensus between pairs of agents. We show that for a large number of issues and unidirectional information flow, global consensus is a fixed point; however, the approach to this consensus is glassy for large societies.

SiO(2) in density functional theory and beyond

We present the first-principle electronic structure calculation on an amorphous material including many-body corrections within the GW approximation. We show that the inclusion of the local field effects in the exchange-correlation potential is crucial to quantitatively describe amorphous systems and defect states. We show that the mobility gap of amorphous silica coincides with the band gap of quartz, contrary to the traditional picture and the densityfunctional theory results. (C) 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

A pattern recognition approach to complex networks

Complex networks exist in many areas of science such as biology, neuroscience, engineering, and sociology. The growing development of this area has led to the introduction of several topological and dynamical measurements, which describe and quantify the structure of networks. Such characterization is essential not only for the modeling of real systems but also for the study of dynamic processes that may take place in them. However, it is not easy to use several measurements for the analysis of complex networks, due to the correlation between them and the difficulty of their visualization. To overcome these limitations, we propose an effective and comprehensive approach for the analysis of complex networks, which allows the visualization of several measurements in a few projections that contain the largest data variance and the classification of networks into three levels of detail, vertices, communities, and the global topology. We also demonstrate the efficiency and the universality of the proposed methods in a series of real-world networks in the three levels.

On exact time averages of a massive Poisson particle

In this work we study, under the Stratonovich definition, the problem of the damped oscillatory massive particle subject to a heterogeneous Poisson noise characterized by a rate of events, lambda(t), and a magnitude, Phi, following an exponential distribution. We tackle the problem by performing exact time averages over the noise in a similar way to previous works analysing the problem of the Brownian particle. From this procedure we obtain the long-term equilibrium distributions of position and velocity as well as analytical asymptotic expressions for the injection and dissipation of energy terms. Considerations on the emergence of stochastic resonance in this type of system are also set forth.

Birkhoffian Systems in Infinite Dimensional Manifolds

We generalize the theory of Kobayashi and Oliva (On the Birkhoff Approach to Classical Mechanics. Resenhas do Instituto de Matematica e Estatistica da Universidade de Sao Paulo, 2003) to infinite dimensional Banach manifolds with a view towards applications in partial differential equations.