基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

Colliding and Moving Bose-Einstein Condensates: Studies of superfluidity and optical tweezers for condensate transport

In this thesis, two different sets of experiments are described. The first is an exploration of the microscopic superfluidity of dilute gaseous Bose- Einstein condensates. The second set of experiments were performed using transported condensates in a new BEC apparatus. Superfluidity was probed by moving impurities through a trapped condensate. The impurities were created using an optical Raman transition, which transferred a small fraction of the atoms into an untrapped hyperfine state. A dramatic reduction in the collisions between the moving impurities and the condensate was observed when the velocity of the impurities was close to the speed of sound of the condensate. This reduction was attributed to the superfluid properties of a BEC. In addition, we observed an increase in the collisional density as the number of impurity atoms increased. This enhancement is an indication of bosonic stimulation by the occupied final states. This stimulation was observed both at small and large velocities relative to the speed of sound. A theoretical calculation of the effect of finite temperature indicated that collision rate should be enhanced at small velocities due to thermal excitations. However, in the current experiments we were insensitive to this effect. Finally, the factor of two between the collisional rate between indistinguishable and distinguishable atoms was confirmed. A new BEC apparatus that can transport condensates using optical tweezers was constructed. Condensates containing 10-15 million sodium atoms were produced in 20 s using conventional BEC production techniques. These condensates were then transferred into an optical trap that was translated from the ‘production chamber’ into a separate vacuum chamber: the ‘science chamber’. Typically, we transferred 2-3 million condensed atoms in less than 2 s. This transport technique avoids optical and mechanical constrainsts of conventional condensate experiments and allows for the possibility of novel experiments. In the first experiments using transported BEC, we loaded condensed atoms from the optical tweezers into both macroscopic and miniaturized magnetic traps. Using microfabricated wires on a silicon chip, we observed excitation-less propagation of a BEC in a magnetic waveguide. The condensates fragmented when brought very close to the wire surface indicating that imperfections in the fabrication process might limit future experiments. Finally, we generated a continuous BEC source by periodically replenishing a condensate held in an optical reservoir trap using fresh condensates delivered using optical tweezers. More than a million condensed atoms were always present in the continuous source, raising the possibility of realizing a truly continuous atom lase.

Finite size effects in the presence of a chemical potential: A study in the classical nonlinear O(2) sigma model

In the presence of a chemical potential, the physics of level crossings leads to singularities at zero temperature, even when the spatial volume is finite. These singularities are smoothed out at a finite temperature but leave behind nontrivial finite size effects which must be understood in order to extract thermodynamic quantities using Monte Carlo methods, particularly close to critical points. We illustrate some of these issues using the classical nonlinear O(2) sigma model with a coupling β and chemical potential μ on a 2+1-dimensional Euclidean lattice. In the conventional formulation this model suffers from a sign problem at nonzero chemical potential and hence cannot be studied with the Wolff cluster algorithm. However, when formulated in terms of the worldline of particles, the sign problem is absent, and the model can be studied efficiently with the "worm algorithm." Using this method we study the finite size effects that arise due to the chemical potential and develop an effective quantum mechanical approach to capture the effects. As a side result we obtain energy levels of up to four particles as a function of the box size and uncover a part of the phase diagram in the (β,μ) plane. © 2010 The American Physical Society.

Reconstructing the dynamics of a movable mirror in a detuned optical cavity

We consider the dynamics of a movable mirror in a Fabry-Perot cavity coupled through radiation pressure to the cavity field and in contact with a thermal bath at finite temperature. In contrast to previous approaches, we consider arbitrary values of the effective detuning between the cavity and an external input field. We analyse the radiation-pressure effect on the Brownian motion of the mirror and its significance in the density noise spectrum of the output cavity field. Important properties of the mirror dynamics can be gathered directly from this noise spectrum. The presented reconstruction provides an experimentally useful tool in the characterization of the energy and rigidity of the mirror as modified by the coupling with light. We also give a quantitative analysis of the recent experimental observation of self-cooling of a micromechanical oscillator.

Panel unit root tests in the presence of cross-sectional dependence: finite sample performance and an application

This paper examines the finite sample properties of three testing regimes for the null hypothesis of a panel unit root against stationary alternatives in the presence of cross-sectional correlation. The regimes of Bai and Ng (2004), Moon and Perron (2004) and Pesaran (2007) are assessed in the presence of multiple factors and also other non-standard situations. The behaviour of some information criteria used to determine the number of factors in a panel is examined and new information criteria with improved properties in small-N panels proposed. An application to the efficient markets hypothesis is also provided. The null hypothesis of a panel random walk is not rejected by any of the tests, supporting the efficient markets hypothesis in the financial services sector of the Australian Stock Exchange.

Multipartite nonlocality in a thermalized Ising spin chain

We study multipartite correlations and nonlocality in an isotropic Ising ring under transverse magnetic field at both zero and finite temperature. We highlight parity-induced differences between the multipartite Bell-like functions used in order to quantify the degree of nonlocality within a ring state and reveal a mechanism for the passive protection of multipartite quantum correlations against thermal spoiling effects that is clearly related to the macroscopic properties of the ring model.

Entangling two distant oscillators with a quantum reservoir

The generation of entanglement between two oscillators that interact via a common reservoir is theoretically studied. The reservoir is modeled by a one-dimensional harmonic crystal initially in thermal equilibrium. Starting from a separable state, the oscillators can become entangled after a transient time, that is of the order of the thermalization time scale. This behaviour is observed at finite temperature even when the oscillators are at a distance significantly larger than the crystal's interparticle spacing. The underlying physical mechanisms can be explained by the dynamical properties of the collective variables of the two oscillators which may decouple from or be squeezed by the reservoir. Our predictions can be tested with an ion chain in a linear Paul trap. Copyright (C) EPLA, 2011

Studies on grain refinement and alloying additions on the microstructure and mechanical properties of Mg-8Zn-4Al alloy

MAGNESIUM ALLOYS have strong potential for weight reduction in a wide range of technical applications because of their low density compared to other structural metallic materials. Therefore, an extensive growth of magnesium alloys usage in the automobile sector is expected in the coming years to enhance the fuel efficiency through mass reduction. The drawback associated with the use of commercially cheaper Mg-Al based alloys, such as AZ91, AM60 and AM50 are their inferior creep properties above 100ºC due to the presence of discontinuous Mg17A112 phases at the grain boundaries. Although rare earth-based magnesium alloys show better mechanical properties, it is not economically viable to use these alloys in auto industries. Recently, many new Mg-Al based alloy systems have been developed for high temperature applications, which do not contain the Mg17Al12 phase. It has been proved that the addition of a high percentage of zinc (which depends upon the percentage of Al) to binary Mg-Al alloys also ensures the complete removal of the Mg17Al12 phase and hence exhibits superior high temperature properties.ZA84 alloy is one such system, which has 8%Zn in it (Mg-8Zn-4Al-0.2Mn, all are in wt %) and shows superior creep resistance compared to AZ and AM series alloys. These alloys are mostly used in die casting industries. However, there are certain large and heavy components, made up of this alloy by sand castings that show lower mechanical properties because of their coarse microstructure. Moreover, further improvement in their high temperature behaviour through microstructural modification is also an essential task to make this alloy suitable for the replacement of high strength aluminium alloys used in automobile industry. Grain refinement is an effective way to improve the tensile behaviour of engineering alloys. In fact, grain refinement of Mg-Al based alloys is well documented in literature. However, there is no grain refiner commercially available in the market for Mg-Al alloys. It is also reported in the literature that the microstructure of AZ91 alloy is modified through the minor elemental additions such as Sb, Si, Sr, Ca, etc., which enhance its high temperature properties because of the formation of new stable intermetallics. The same strategy can be used with the ZA84 alloy system to improve its high temperature properties further without sacrificing the other properties. The primary objective of the present research work, “Studies on grain refinement and alloying additions on the microstructure and mechanical properties of Mg-8Zn-4Al alloy” is twofold: 1. To investigate the role of individual and combined additions of Sb and Ca on the microstructure and mechanical properties of ZA84 alloy. 2. To synthesis a novel Mg-1wt%Al4C3 master alloy for grain refinement of ZA84 alloy and investigate its effects on mechanical properties.

Estimation of temperature effects on fission barriers

An accurate mass formula at finite temperature has been used to obtain a more precise estimation of temperature effects on fission barriers calculated within the liquid drop model.

Correlations in hot asymmetric nuclear matter

Bulk and single-particle properties of hot hyperonic matter are studied within the Brueckner-Hartree-Fock approximation extended to finite temperature. The bare interaction in the nucleon sector is the Argonne V18 potential supplemented with an effective three-body force to reproduce the saturating properties of nuclear matter. The modern Nijmegen NSC97e potential is employed for the hyperon-nucleon and hyperon-hyperon interactions. The effect of temperature on the in-medium effective interaction is found to be, in general, very small and the single-particle potentials differ by at most 25% for temperatures in the range from 0 to 60 MeV. The bulk properties of infinite matter of baryons, either nuclear isospin symmetric or a Beta-stable composition that includes a nonzero fraction of hyperons, are obtained. It is found that the presence of hyperons can modify the thermodynamical properties of the system in a non-negligible way.

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.

Temperature effects on a network of dissipative quantum harmonic oscillators: collective damping and dispersion processes

In this paper we extend the results presented in (de Ponte, Mizrahi and Moussa 2007 Phys. Rev. A 76 032101) to treat quantitatively the effects of reservoirs at finite temperature in a bosonic dissipative network: a chain of coupled harmonic oscillators whatever its topology, i.e., whichever the way the oscillators are coupled together, the strength of their couplings and their natural frequencies. Starting with the case where distinct reservoirs are considered, each one coupled to a corresponding oscillator, we also analyze the case where a common reservoir is assigned to the whole network. Master equations are derived for both situations and both regimes of weak and strong coupling strengths between the network oscillators. Solutions of these master equations are presented through the normal ordered characteristic function. These solutions are shown to be significantly involved when temperature effects are considered, making difficult the analysis of collective decoherence and dispersion in dissipative bosonic networks. To circumvent these difficulties, we turn to the Wigner distribution function which enables us to present a technique to estimate the decoherence time of network states. Our technique proceeds by computing separately the effects of dispersion and the attenuation of the interference terms of the Wigner function. A detailed analysis of the dispersion mechanism is also presented through the evolution of the Wigner function. The interesting collective dispersion effects are discussed and applied to the analysis of decoherence of a class of network states. Finally, the entropy and the entanglement of a pure bipartite system are discussed.

HIGH-TEMPERATURE SUPERFLUIDITY WITH ABNORMAL FERMIONIC OCCUPANCY

We examine a Lipkin based two-level pairing model at finite temperature and in the thermodynamic limit. Whereas at T = 0 the model exhibits a superconducting ground state for sufficiently high values of the coupling constant, a partially superconducting phase in which some of the particles are paired, is found to survive at high temperatures in a special treatment. This phase is a mixture of abnormally-occupied eigenstates, which lie at higher energy, of the interactionless model Hamiltonian.

The temperature dependence of the QCD running coupling

We study the running of the QCD coupling with the momentum squared (Q 2) and the temperature scales in the high temperature limit (T > Tc), using a mass dependent renormalization scheme to build the Renormalization Group Equations. The approach used guaranty gauge invariance, through the use of the Hard Thermal Loop approximation, and independence of the vertex chosen to renormalize the coupling. In general, the dependence of the coupling with the temperature is not logarithmical, although in the region Q2 ∼ T2 the logarithm approximation is reasonable. Finally, as known from Debye screening, color charge is screened in the coupling. The number of flavors, however, is anti-screened.

Critical currents and melting temperature of a two-dimensional vortex lattice with periodic pinning

The critical current and melting temperature of a vortex system are analyzed. Calculations are made for a two-dimensional film at finite temperature with two kinds of periodic pinning: hexagonal and Kagomé. A transport current parallel and perpendicular to the main axis of the pinning arrays is applied and molecular dynamics simulations are used to calculate the vortex velocities to obtain the critical currents. The structure factor and displacements of vortices at zero transport current are used to obtain the melting temperature for both pinning arrays. The critical currents are higher for the hexagonal pinning lattice and anisotropic for both pinning arrays. This anisotropy is stronger with temperature for the hexagonal array. For the Kagomé pinning lattice, our analysis shows a multi stage phase melting; that is, as we increase the temperature, each different dynamic phase melts before reaching the melting temperature. Both the melting temperature and critical currents are larger for the hexagonal lattice, indicating the role for the interstitial vortices in decreasing the pinning strength. © 2012 Springer Science+Business Media New York.