Second quantized reduced Bloch equations and the exact solutions for pairing Hamiltonian

In this article, we present a set of hierarchy Bloch equations for the reduced density operators in either canonical or grand canonical ensembles in the occupation number representation. They provide a convenient tool for studying the equilibrium quantum statistical mechanics for some model systems. As an example of their applications, we solve the equations for the model system with a pairing Hamiltonian. With the aid of its symplectic group symmetry, we obtain the statistical reduced density matrices with different orders. As a special instance for the solutions, we also get the reduced density matrices of the ground state for a superconductor.     

Quantum information transmission between two qubits through an intermediary photon gas

The theory of the quantum information transmission between two semiconductor two-level quantum dots as two qubits through an intermediary photon gas in a cavity is presented. The reduced density matrix of each two-level quantum dot is the quantum information encoded into this qubit. The quantum information exchange between two distant qubits imbedded in the photon gas is performed in the form of the mutual dependence of their reduced density matrices due to the interaction between the electrons in the qubits and the photon gas. The system of rate equations for the reduced density matrix of the two-qubit system is derived. From the solution of this system of equations it follows the mutual dependence of the reduced density matrices of two distant qubits.     

A core-softened fluid model in disordered porous media. Grand canonical Monte Carlo simulation and integral equations

We have studied the microscopic structure and thermodynamic properties of a core-softened fluid model in disordered matrices of Lennard-Jones particles by using grand canonical Monte Carlo simulation. The dependence of density on the applied chemical potential (adsorption isotherms), pair distribution functions, as well as the heat capacity in different matrices are discussed. The microscopic structure of the model in matrices changes with density similar to the bulk model. Thus one should expect that the structural anomaly persists at least in dilute matrices. The region of densities for the heat capacity anomaly shrinks with increasing matrix density. This behavior is also observed for the diffusion coefficient on density from independent molecular dynamics simulation. Theoretical results for the model have been obtained by using replica Ornstein-Zernike integral equations with hypernetted chain closure. Predictions of the theory generally are in good agreement with simulation data, except for the heat capacity on fluid density. However, possible anomalies of thermodynamic properties for the model in disordered matrices are not captured adequately by the present theory. It seems necessary to develop and apply more elaborated, thermodynamically self-consistent closures to capture these features.     

Reduced density matrices of the anisotropic Heisenberg model

The reduced density matrices of the anisotropic Heisenberg model are studied by means of a functional integral representation based on a generalized Poisson process. Integral equations, which are analogous to the classical Kirkwood-Salzburg equations, are then used to prove the existence of the infinite volume limit of the reduced density matrices, analyticity properties with respect to the fugacity (or magnetic field) and the potentials, and a cluster property, in the low fugacity (high magnetic field) region.The research reported in this paper was supported by the National Science Foundation.     

Hierarchy of equations for reduced density matrices in the case of thermodynamic equilibrium

A hierarchy of equations for s-particle density matrices at thermodynamic equilibrium is obtained, with the equation for the nonequilibrium density matrix as the starting point. When deducing the hierarchy the hypothesis of maximum statistical independence for the density matrices is used. The hierarchy obtained is an analogue of the classical equilibrium BBGKY hierarchy and goes over into it when . It is shown that thermodynamic quantities can be expressed in terms of functions that enter only into the first hierarchy equations. The hierarchy is analysed in detail in the case of a uniform fluid. As an example in which the equations can be solved easily enough, a hard-sphere system wherein triplet correlations are neglected is considered. Different approximations that can be used when solving the equations derived are discussed. Comparisons are made with the results of other theoretical treatments.     

Dissipative evolution of quantum statistical ensembles and nonlinear response to a time-periodic perturbation

We present a detailed discussion of the evolution of a statistical ensemble of quantum mechanical systems coupled weakly to a bath. The Hilbert space of the full system is given by the tensor product between the Hilbert spaces associated with the bath and the bathed system. The statistical states of the ensemble are described in terms of density matrices. Supposing the bath to be held at some - not necessarily thermal - statistical equilibrium and tracing over the bath degrees of freedom, we obtain reduced density matrices defining the statistical states of the bathed system. The master equations describing the evolution of these reduced density matrices are derived under the most general conditions. On time scales that are large with respect to the bath correlation time and with respect to the reciprocal transition frequencies of the bathed system, the resulting evolution of the reduced density matrix of the bathed system is of Markovian type. The detailed balance relations valid for a thermal equilibrium of the bath are derived and the conditions for the validity of the fluctuation-dissipation theorem are given. Based on the general approach, we investigate the non-linear response of the bathed subsystem to a time-periodic perturbation. Summing the perturbation series we obtain the coherences and the populations for arbitrary strengths of the perturbation.Received: 26 November 2003, Published online: 30 January 2004PACS: 05.30.-d Quantum statistical mechanics - 33.35. + r Electron resonance and relaxation - 33.25. + k Nuclear resonance and relaxation     

Functional equations for path integrals

We consider the density matrices that arise in the statistical mechanics of the electron-phonon systems. In the path integral representation the phonon coordinates can be eliminated. This leads to an action that depends on pairs of points on a path, that depends explicitly on time differences, and that contains the phonon occupation numbers. The integral is reduced to a standard form by scaling to the thermal length. We use the technique of integration by parts and add specially chosen generating functionals to the action. We set down functional derivative equations for the source-dependent density matrix and for the mass operator. This allows us to develop a series of approximations for the operator in terms of exact propagators. The crudest approximation is a coherent potential approximation applicable at a general temperature.     

量纲分析法求解辐射热传导自相似解

通过量纲分析,将辐射流体方程简化为一阶常微分方程组的形式,进而结合数值计算求得热波烧蚀自相似解。流体方程涉及的内能、状态方程表达式采用依赖于温度、密度两个参量的形式,且参数可调,适用性较广。源边界条件可为恒温、恒流、恒压、常密等情况。以恒温源为例,给出了X射线加热Au壁时烧蚀质量与烧蚀压的定标关系,与文献结果及MULTI辐射流体程序计算结果吻合较好。     

A dynamics-controlled truncation scheme for the hierarchy of density matrices in semiconductor optics

We discuss the interaction of coherent electromagnetic fields with the semiconductor band edge in a dynamic density matrix model. Due to the influence of the Coulomb-interaction then-point density matrices are coupled in an infinite hierarchy of equations of motion. We show how this hierarchy is related to an expansion of the density matrices in terms of powers of the exciting field. We make use of the above results to set up a closed set of equations of motion involving two-, four-and six-point correlation functions, from which all third order contributions to the polarization can be calculated exactly. Comparison of our treatment of the hierarchy with the widely used RPA decoupling on the two-point level, gives interesting insight into the validity of the RPA. In particular we find, that a RPA-like factorization for two of the relevant density-matrices yields a solution of their respective equations of motion to lowest order in the electric field.     

Will a bose-condensed exciton gas be superfluid?

For a system of Bose-condensed Wannier excitons, the second reduced density matrix for electrons and holes is shown to have the property of off-diagonal long-range order (ODLRO). From the equation of motion for this density matrix, we derive the conservation laws which are relevant for a condensed exciton system. These equations are obtained by projecting the electron-hole density matrices into the exciton space. From the conservation laws, a two-fluid model follows, which describes the superfluid flow of the excitation energy.     

Maximum-entropy closures for kinetic theories of neuronal network dynamics

We analyze (1 + 1)D kinetic equations for neuronal network dynamics, which are derived via an intuitive closure from a Boltzmann-like equation governing the evolution of a one-particle (i.e., one-neuron) probability density function. We demonstrate that this intuitive closure is a generalization of moment closures based on the maximum-entropy principle. By invoking maximum-entropy closures, we show how to systematically extend this kinetic theory to obtain higher-order, kinetic equations and to include coupled networks of both excitatory and inhibitory neurons.     

Masterequation,H-theorem and stationary solution for coupled quantum systems

We extend the Zwanzig projector formalism to coupled systems taking into account the mutual interactions of the reduced density matrices of both systems. In the Born- and Markoff-approximation we end up with a bilinear masterequation for occupation probabilities, in contrast to the usually studied linear equations. We derive theH-theorem for this equation and show that the stationary solution is the canonical or more generally a grand canonical density matrix.     

Decay of correlations in spin systems

We investigate the decay of initial correlations in a spin system where each spin relaxes independently by an intramolecular mechanism. The equation of motion for the spin density matrix is assumed to be the Redfield equation, which is of the form of a quantum mechanical master equation. Our analysis of this problem is based on the techniques of Shuler, Oppenheim, and coworkers, who have studied the decay of correlations in systems which can be described by classical stochastic master equations. We find that the off-diagonal elements of the reduced spin density matrices approach their equilibrium values faster than the diagonal elements. The Ursell functions, which are a measure of the correlations in the system, decay to their zero equilibrium values faster than the spin density matrix except for the furthest off-diagonal elements. Far off-diagonal matrix elements of the spin density matrix approach equilibrium at the same rate as the Ursell functions, which is the important difference between the quantum mechanical model studied here and the classical models studied earlier.Supported in part by the National Science Foundation.     

The BBGKY hierarchy in quantum statistical mechanics

It is shown that the infinite volume limit of the equilibrium reduced density matrices, shown by Ginibre to exist at low densities, satisfy the quantum time independent BBGKY hierarchy of equations. This extends analogous results for classical systems due to Gallavotti.Supported by AFOSR Contract Number F44620-71-C-0013.     

Solving Master Equation for Two-Mode Density Matrices by Virtue of Thermal Entangled State Representation

We extend the approach of solving master equations for density matrices by projecting it onto the thermal entangled state representation (Hong-Yi Fan and Jun-Hua Chen, J. Phys. A35 (2002) 6873) to two-mode case. In this approach the two-photon master equations can be directly and conveniently converted into c-number partial differential equations. As an example, we solve the typical master equation for two-photon process in some limiting cases.     

A direct solver for the least-squares spectral collocation system on rectangular elements for the incompressible Navier–Stokes equations

A least-squares spectral collocation scheme for the Stokes and incompressible Navier–Stokes equations is proposed. The original domain is decomposed into quadrilateral subelements and on the element interfaces continuity of the functions is enforced in the least-squares sense. The collocation conditions and the interface conditions lead to overdetermined systems. These systems are directly solved by QR decomposition of the underlying matrices. By numerical simulations it is shown that the direct method leads to better results than the approach with normal equations. Furthermore, it is shown that the condition numbers can be reduced by introducing the Clenshaw–Curtis quadrature rule for imposing the average pressure to be zero. Finally, our scheme is successfully applied to the regularized and lid-driven cavity flow problems.     

Methodology for the solutions of some reduced Fokker‐Planck equations in high dimensions

In this paper, a new methodology is formulated for solving the reduced Fokker‐Planck (FP) equations in high dimensions based on the idea that the state space of large‐scale nonlinear stochastic dynamic system is split into two subspaces. The FP equation relevant to the nonlinear stochastic dynamic system is then integrated over one of the subspaces. The FP equation for the joint probability density function of the state variables in another subspace is formulated with some techniques. Therefore, the FP equation in high‐dimensional state space is reduced to some FP equations in low‐dimensional state spaces, which are solvable with exponential polynomial closure method. Numerical results are presented and compared with the results from Monte Carlo simulation and those from equivalent linearization to show the effectiveness of the presented solution procedure. It attempts to provide an analytical tool for the probabilistic solutions of the nonlinear stochastic dynamics systems arising from statistical mechanics and other areas of science and engineering.     

Equations for generating functionals in classical and quantum statistical mechanics

This paper develops a formalism for the generating functionals for partial distribution functions in classical statistical mechanics and partial density matrices in quantum statistical mechanics. For the case of a large canonical ensemble, functional equations are written with respect to the functionals introduced. Each functional system creates a system of integral equations for distribution functions or density matrices.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii Fizika, No. 9, pp. 98–102, September, 1971.