Elliptical slice sampling

Many probabilistic models introduce strong dependencies between variables using a latent multivariate Gaussian distribution or a Gaussian process. We present a new Markov chain Monte Carlo algorithm for performing inference in models with multivariate Gaussian priors. Its key properties are: 1) it has simple, generic code applicable to many models, 2) it has no free parameters, 3) it works well for a variety of Gaussian process based models. These properties make our method ideal for use while model building, removing the need to spend time deriving and tuning updates for more complex algorithms.

Size effects on the melting of nickel nanowires: a molecular dynamics study

The melting process of nickel nanowires are simulated by using molecular dynamics with the quantum Sutten-Chen many-body force field. The wires studied were approximately cylindrical in cross-section and periodic boundary conditions were applied along their length; the atoms were arranged initially in a face-centred cubic structure with the [0 0 1] direction parallel to the long axis of the wire. The size effects of the nanowires on the melting temperatures are investigated. We find that for the nanoscale regime, the melting temperatures of Ni nanowires are much lower than that of the bulk and are linear with the reciprocal of the diameter of the nanowire. When a nanowire is heated up above the melting temperature, the neck of the nanowire begins to arise and the diameter of neck decreases rapidly with the equilibrated running time. Finally, the breaking of nanowire arises, which leads to the formation of the spherical clusters. (C) 2004 Elsevier B.V. All rights reserved.

Molecular dynamics studies on the dislocation gliding near a tilt boundary

The gliding behavior of edge dislocation near a grain boundary(QB) in copper under pure shear stresses is simulated by using molecular dynamics(MD) method. Many-body potential incorporating the embedded atom method (EAM) is used. The critical shear stresses for a single disocation to pass across GB surface are obtained at values of sigma(c)=23MPa similar to 68 MPa and 137 MPa similar to 274 MPa for Sigma=165 small angle tilt GB at 300 K and 20 K, respectively. The first result agrees with the experimental yield stress sigma(y)(=42 MPa) quite well. It suggests that there might be one of the reasons of initial plastic yielding caused by single dislocation gliding across GB. In addition, there might be possibility to obtain yield strength from microscopic analysis. Moreover, the experimental value of sigma(y) at low temperature is generally higher than that at room temperature. So, these results are in conformity qualitatively with experimental fact. On the other hand, the Sigma=25 GB is too strong an obstacle to the dislocation. In this case, a dislocation is able to pass across GB under relatively low stress only when it is driven by other dislocations. This is taken to mean that dislocation pile-up must be built up in front of this kind of GB, if this GB may take effect on the process of plastic deformation.

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.

Fermionic quantum systems. Part I: Phase transitions in quantum dots. Part II: Nuclear matter on a lattice

In the first part I perform Hartree-Fock calculations to show that quantum dots (i.e., two-dimensional systems of up to twenty interacting electrons in an external parabolic potential) undergo a gradual transition to a spin-polarized Wigner crystal with increasing magnetic field strength. The phase diagram and ground state energies have been determined. I tried to improve the ground state of the Wigner crystal by introducing a Jastrow ansatz for the wave function and performing a variational Monte Carlo calculation. The existence of so called magic numbers was also investigated. Finally, I also calculated the heat capacity associated with the rotational degree of freedom of deformed many-body states and suggest an experimental method to detect Wigner crystals.

The second part of the thesis investigates infinite nuclear matter on a cubic lattice. The exact thermal formalism describes nucleons with a Hamiltonian that accommodates on-site and next-neighbor parts of the central, spin-exchange and isospin-exchange interaction. Using auxiliary field Monte Carlo methods, I show that energy and basic saturation properties of nuclear matter can be reproduced. A first order phase transition from an uncorrelated Fermi gas to a clustered system is observed by computing mechanical and thermodynamical quantities such as compressibility, heat capacity, entropy and grand potential. The structure of the clusters is investigated with the help two-body correlations. I compare symmetry energy and first sound velocities with literature and find reasonable agreement. I also calculate the energy of pure neutron matter and search for a similar phase transition, but the survey is restricted by the infamous Monte Carlo sign problem. Also, a regularization scheme to extract potential parameters from scattering lengths and effective ranges is investigated.

Theory of s-d exchange scattering in dilute magnetic alloys

The problem of s-d exchange scattering of conduction electrons off localized magnetic moments in dilute magnetic alloys is considered employing formal methods of quantum field theoretical scattering. It is shown that such a treatment not only allows for the first time, the inclusion of multiparticle intermediate states in single particle scattering equations but also results in extremely simple and straight forward mathematical analysis. These equations are proved to be exact in the thermodynamic limit. A self-consistent integral equation for electron self energy is derived and approximately solved. The ground state and physical parameters of dilute magnetic alloys are discussed in terms of the theoretical results. Within the approximation of single particle intermediate states our results reduce to earlier versions. The following additional features are found as a consequence of the inclusion of multiparticle intermediate states;

(i) A non analytic binding energy is pre sent for both, antiferromagnetic (J < o) and ferromagnetic (J > o) couplings of the electron plus impurity system.

(ii) The correct behavior of the energy difference of the conduction electron plus impurity system and the free electron system is found which is free of unphysical singularities present in earlier versions of the theories.

(iii) The ground state of the conduction electron plus impurity system is shown to be a many-body condensate state for J < o and J > o, both. However, a distinction is made between the usual terminology of "Singlet" and "Triplet" ground states and nature of our ground state.

(iv) It is shown that a long range ordering, leading to an ordering of the magnetic moments can result from a contact interaction such as the s-d exchange interaction.

(v) The explicit dependence of the excess specific heat of the Kondo systems is obtained and found to be linear in temperatures as T→ o and T ℓnT for 0.3 T_K ≤ T ≤ 0.6 T_K. A rise in (ΔC/T) for temperatures in the region 0 < T ≤ 0.1 T_K is predicted. These results are found to be in excellent agreement with experiments.

(vi) The existence of a critical temperature for Ferromagnetic coupling (J > o) is shown. On the basis of this the apparent contradiction of the simultaneous existence of giant moments and Kondo effect is resolved.

Cavity Optomechanics at Millikelvin Temperatures

The field of cavity optomechanics, which concerns the coupling of a mechanical object's motion to the electromagnetic field of a high finesse cavity, allows for exquisitely sensitive measurements of mechanical motion, from large-scale gravitational wave detection to microscale accelerometers. Moreover, it provides a potential means to control and engineer the state of a macroscopic mechanical object at the quantum level, provided one can realize sufficiently strong interaction strengths relative to the ambient thermal noise. Recent experiments utilizing the optomechanical interaction to cool mechanical resonators to their motional quantum ground state allow for a variety of quantum engineering applications, including preparation of non-classical mechanical states and coherent optical to microwave conversion. Optomechanical crystals (OMCs), in which bandgaps for both optical and mechanical waves can be introduced through patterning of a material, provide one particularly attractive means for realizing strong interactions between high-frequency mechanical resonators and near-infrared light. Beyond the usual paradigm of cavity optomechanics involving isolated single mechanical elements, OMCs can also be fashioned into planar circuits for photons and phonons, and arrays of optomechanical elements can be interconnected via optical and acoustic waveguides. Such coupled OMC arrays have been proposed as a way to realize quantum optomechanical memories, nanomechanical circuits for continuous variable quantum information processing and phononic quantum networks, and as a platform for engineering and studying quantum many-body physics of optomechanical meta-materials.

However, while ground state occupancies (that is, average phonon occupancies less than one) have been achieved in OMC cavities utilizing laser cooling techniques, parasitic absorption and the concomitant degradation of the mechanical quality factor fundamentally limit this approach. On the other hand, the high mechanical frequency of these systems allows for the possibility of using a dilution refrigerator to simultaneously achieve low thermal occupancy and long mechanical coherence time by passively cooling the device to the millikelvin regime. This thesis describes efforts to realize the measurement of OMC cavities inside a dilution refrigerator, including the development of fridge-compatible optical coupling schemes and the characterization of the heating dynamics of the mechanical resonator at sub-kelvin temperatures.

We will begin by summarizing the theoretical framework used to describe cavity optomechanical systems, as well as a handful of the quantum applications envisioned for such devices. Then, we will present background on the design of the nanobeam OMC cavities used for this work, along with details of the design and characterization of tapered fiber couplers for optical coupling inside the fridge. Finally, we will present measurements of the devices at fridge base temperatures of T_f = 10 mK, using both heterodyne spectroscopy and time-resolved sideband photon counting, as well as detailed analysis of the prospects for future quantum applications based on the observed optically-induced heating.

Strong atom-light interactions along nanostructures: Transition from free-space to nanophotonic interfaces

An exciting frontier in quantum information science is the integration of otherwise "simple'' quantum elements into complex quantum networks. The laboratory realization of even small quantum networks enables the exploration of physical systems that have not heretofore existed in the natural world. Within this context, there is active research to achieve nanoscale quantum optical circuits, for which atoms are trapped near nano-scopic dielectric structures and "wired'' together by photons propagating through the circuit elements. Single atoms and atomic ensembles endow quantum functionality for otherwise linear optical circuits and thereby enable the capability of building quantum networks component by component. Toward these goals, we have experimentally investigated three different systems, from conventional to rather exotic systems : free-space atomic ensembles, optical nano fibers, and photonics crystal waveguides. First, we demonstrate measurement-induced quadripartite entanglement among four quantum memories. Next, following the landmark realization of a nanofiber trap, we demonstrate the implementation of a state-insensitive, compensated nanofiber trap. Finally, we reach more exotic systems based on photonics crystal devices. Beyond conventional topologies of resonators and waveguides, new opportunities emerge from the powerful capabilities of dispersion and modal engineering in photonic crystal waveguides. We have implemented an integrated optical circuit with a photonics crystal waveguide capable of both trapping and interfacing atoms with guided photons, and have observed the collective effect, superradiance, mediated by the guided photons. These advances provide an important capability for engineered light-matter interactions, enabling explorations of novel quantum transport and quantum many-body phenomena.

Entanglement negativity in the harmonic chain out of equilibrium

We study the entanglement in a chain of harmonic oscillators driven out of equilibrium by preparing the two sides of the system at different temperatures, and subsequently joining them together. The steady state is constructed explicitly and the logarithmic negativity is calculated between two adjacent segments of the chain. We find that, for low temperatures, the steady-state entanglement is a sum of contributions pertaining to left-and right-moving excitations emitted from the two reservoirs. In turn, the steady-state entanglement is a simple average of the Gibbs-state values and thus its scaling can be obtained from conformal field theory. A similar averaging behaviour is observed during the entire time evolution. As a particular case, we also discuss a local quench where both sides of the chain are initialized in their respective ground states.

Symmetry restoration and quantumness reestablishment

A realistic quantum many-body system, characterized by a generic microscopic Hamiltonian, is accessible only through approximation methods. The mean field theories, as the simplest practices of approximation methods, commonly serve as a powerful tool, but unfortunately often violate the symmetry of the Hamiltonian. The conventional BCS theory, as an excellent mean field approach, violates the particle number conservation and completely erases quantumness characterized by concurrence and quantum discord between different modes. We restore the symmetry by using the projected BCS theory and the exact numerical solution and find that the lost quantumness is synchronously reestablished. We show that while entanglement remains unchanged with the particle numbers, quantum discord behaves as an extensive quantity with respect to the system size. Surprisingly, discord is hardly dependent on the interaction strengths. The new feature of discord offers promising applications in modern quantum technologies.

Density-potential mapping in the standard and quantum electrodynamical time-dependent density functional theory

131 p.

Distribuição da caspase-3 clivada no cerebelo de ratos durante o desenvolvimento em um modelo de hipóxia-isquemia sistêmica pré-natal

A Hipóxia-isquemia (HI) perinatal é um problema de saúde pública, e ocorrem aproximadamente 1,5 casos de encefalopatias por HI por 1000 nascidos vivos. Dos que sobrevivem 25-60% sofrem de deficiências permanentes do desenvolvimento neurológico, incluindo paralisia cerebral, convulsões, retardo mental, e dificuldade de aprender. Neurônios e oligodendrócitos, especialmente os progenitores, são os mais afetados pela HI. Existem vários modelos de HI, no entanto, poucos levam em consideração as intercorrências maternas, a importância da atividade placentária, e as trocas entre mãe-filho, que são clinicamente observadas em humanos. Robinson estabeleceu um modelo de HI sistêmica pré-natal transitório, onde o fluxo das artérias uterinas da rata grávida era obstruído por 45 minutos no décimo oitavo dia (E18) de gestação. Neste modelo foram observadas alterações que são similares às observadas em cérebros humanos que passaram por hipóxia perinatal, dentre as quais foram relatados aumento no nível de apoptose. Caspase-3 é descrita como uma enzima que atua na apoptose, e é amplamente utilizada como marcador para células apoptóticas. Vários autores vêm mostrando, entretanto, que a enzima caspase-3 pode estar ativada para fins não apoptóticos. No modelo de HI sistêmica pré-natal, foram observados astrogliose na substância branca, morte de oligodendrócitos, lesão em axônios tanto na substância branca como no córtex cerebral, e danos motores. Pouco se sabe da influencia do insulto HI no desenvolvimento do cerebelo, considerando que o cerebelo junto com o córtex motor, contribui para o controle motor. O objetivo desse trabalho foi avaliar a distribuição da caspase-3 clivada durante o desenvolvimento do cerebelo em um modelo de HI pré-natal. Os resultados deste trabalho demonstram que as células caspase-3 clivadas apresentaram duas morfologias distintas em ambos os grupos. Uma onde a caspase-3 foi observada apenas no núcleo, oscilando entre células com imunorreatividade fraca a intensa, e de células com a presença da caspase-3 no corpo celular, nos prolongamentos condensados e presença de fragmentos ao redor do soma, morfologia típica de célula em apoptose. A HI pré-natal, assim como nos hemisférios cerebrais, levou ao aumento de células caspase-3 clivadas com morfologia de progenitores de oligodendrócitos no cerebelo do grupo HI em P2, mas não em P9 e P23. Também foi demonstrado que a HI pré-natal não levou a uma ativação da apoptose em oligodendrócitos, neurônios e microglia (identificados por seus respectivos marcadores, CNPase, NeuN e ED1) apresentando marcação no núcleos de células GFAP+, na substância branca, camada granular e nas células da glia de Bergmann, em P9 e P23 no cerebelo. Podemos concluir que a HI pré-natal aumentou o número de células imunorreativas para a caspase-3 em um período crítico do desenvolvimento da oligodendroglia no cerebelo, e que a diminuição de progenitores de oligodendrócitos no cerebelo decorrente do insulto pré-natal visto em trabalhos anteriores, pode estar relacionada a morte celular por apoptose, embora não se possa descartar a hipótese da participação dessas células que apresentam caspase-3 clivada em outros eventos não apoptóticos desencadeados pela hipóxia-isquemia.

Single-particle probing of edge-state formation in a graphene nanoribbon

We investigate the effect of a perpendicular magnetic field on the single-particle charging spectrum of a graphene quantum dot embedded inline with a nanoribbon. We observe uniform shifts in the single-particle spectrum which coincide with peaks in the magnetoconductance, implicating Landau level condensation and edge state formation as the mechanism underlying magnetic field-enhanced transmission through graphene nanostructures. The experimentally determined ratio of bulk to edge states is supported by single-particle band-structure simulations, while a fourfold beating of the Coulomb blockade transmission amplitude points to many-body interaction effects during Landau level condensation of the ν=0 state. © 2012 American Physical Society.

Accurate screened exchange band structures for the transition metal monoxides MnO, FeO, CoO and NiO.

We report calculations of the band structures and density of states of the four transition metal monoxides MnO, FeO, CoO and NiO using the hybrid density functional sX-LDA ('screened exchange local density approximation'). Late transition metal oxides are prototypical examples of strongly correlated materials, which pose challenges for electronic structure methods. We compare our results with available experimental data and show that our calculations generally yield accurate predictions for the fundamental band gaps and valence bands, in favourable agreement with previously reported theoretical studies. For MnO, the band gaps are still underestimated, suggesting additional many-body effects that are not captured by our screened hybrid functional approach.

Ultrafast collinear scattering and carrier multiplication in graphene.

Graphene is emerging as a viable alternative to conventional optoelectronic, plasmonic and nanophotonic materials. The interaction of light with charge carriers creates an out-of-equilibrium distribution, which relaxes on an ultrafast timescale to a hot Fermi-Dirac distribution, that subsequently cools emitting phonons. Although the slower relaxation mechanisms have been extensively investigated, the initial stages still pose a challenge. Experimentally, they defy the resolution of most pump-probe setups, due to the extremely fast sub-100 fs carrier dynamics. Theoretically, massless Dirac fermions represent a novel many-body problem, fundamentally different from Schrödinger fermions. Here we combine pump-probe spectroscopy with a microscopic theory to investigate electron-electron interactions during the early stages of relaxation. We identify the mechanisms controlling the ultrafast dynamics, in particular the role of collinear scattering. This gives rise to Auger processes, including charge multiplication, which is key in photovoltage generation and photodetectors.