Lattice density-functional theory on graphene

A density-functional approach on the hexagonal graphene lattice is developed using an exact numerical solution to the Hubbard model as the reference system. Both nearest-neighbour and up to third nearest-neighbour hoppings are considered and exchange-correlation potentials within the local density approximation are parameterized for both variants. The method is used to calculate the ground-state energy and density of graphene flakes and infinite graphene sheet. The results are found to agree with exact diagonalization for small systems, also if local impurities are present. In addition, correct ground-state spin is found in the case of large triangular and bowtie flakes out of the scope of exact diagonalization methods.

Nonlinear conduction in charge-ordered Pr0.63Ca0.37MnO3: Effect of magnetic fields

Nonlinear conduction in a single crystal of charge-ordered Pr0.63Ca0.37MnO3 has bren investigated in an applied magnetic field. In zero field, the nonlinear conduction, which starts at T< T-CO, can give rise to a region of negative differential resistance (NDR) which shows up below the Neel temperature. Application of a magnetic field Inhibits the appearance of NDR and makes the nonlinear conduction strongly hysteritic on cycling of the bias current. This is most severe in the temperature range where the charge-ordered state melts in an applied magnetic field. Our experiment strongly suggests that application of a magnetic field in the charge-ordering regime causes a coexistence of two phases.

Spontaneous DC Current Generation in a Resistively Shunted Semiconductor Superlattice Driven by a TeraHertz Field

We study a resistively shunted semiconductor superlattice subject to a high-frequency electric field. Using a balance equation approach that incorporates the influence of the electric circuit, we determine numerically a range of amplitude and frequency of the ac field for which a dc bias and current are generated spontaneously and show that this region is likely accessible to current experiments. Our simulations reveal that the Bloch frequency corresponding to the spontaneous dc bias is approximately an integer multiple of the ac field frequency.

Dynamical Instabilities and Deterministic Chaos in Ballistic Electron Motion in Semiconductor Superlattices

We consider the motion of ballistic electrons within a superlattice miniband under the influence of an alternating electric field. We show that the interaction of electrons with the self-consistent electromagnetic field generated by the electron current may lead to the transition from regular to chaotic dynamics. We estimate the conditions for the experimental observation of this deterministic chaos and discuss the similarities of the superlattice system with the other condensed matter and quantum optical systems.

Symmetry-breaking and chaos in electron transport in semiconductor superlattices

We study the motion of electrons in a single miniband of a semiconductor superlattice driven by THz electric field polarized along the growth direction. We work in the semiclassical balance-equation model, including different elastic and inelastic scattering rates, and incorporating the self-consistent electric field generated by electron motion. We explore regions of complex dynamics, which can include chaotic behaviour and symmetry-breaking. We estimate the magnitudes of dc current and dc voltage that spontaneously appear in regions of broken-symmetry for parameters characteristic of modern semiconductor superlattices. This work complements PRL 80(1998)2669 [ cond-mat/9709026 ].

Dissipative Chaos in Semiconductor Superlattices

We consider the motion of ballistic electrons in a miniband of a semiconductor superlattice (SSL) under the influence of an external, time-periodic electric field. We use the semi-classical balance-equation approach which incorporates elastic and inelastic scattering (as dissipation) and the self-consistent field generated by the electron motion. The coupling of electrons in the miniband to the self-consistent field produces a cooperative nonlinear oscillatory mode which, when interacting with the oscillatory external field and the intrinsic Bloch-type oscillatory mode, can lead to complicated dynamics, including dissipative chaos. For a range of values of the dissipation parameters we determine the regions in the amplitude-frequency plane of the external field in which chaos can occur. Our results suggest that for terahertz external fields of the amplitudes achieved by present-day free electron lasers, chaos may be observable in SSLs. We clarify the nature of this novel nonlinear dynamics in the superlattice-external field system by exploring analogies to the Dicke model of an ensemble of two-level atoms coupled with a resonant cavity field and to Josephson junctions.

Fractional and unquantized dc voltage generation in THz-driven semiconductor superlattices

We consider the spontaneous creation of a dc voltage across a strongly coupled semiconductor superlattice subjected to THz radiation. We show that the dc voltage may be approximately proportional either to an integer or to a half- integer multiple of the frequency of the applied ac field, depending on the ratio of the characteristic scattering rates of conducting electrons. For the case of an ac field frequency less than the characteristic scattering rates, we demonstrate the generation of an unquantized dc voltage.

Possible Exotic Phases in the One-Dimensional Extended Hubbard Model

We investigate numerically the ground state phase diagram of the one-dimensional extended Hubbard model, including an on--site interaction U and a nearest--neighbor interaction V. We focus on the ground state phases of the model in the V >> U region, where previous studies have suggested the possibility of dominant superconducting pairing fluctuations before the system phase separates at a critical value V=V_PS. Using quantum Monte Carlo methods on lattices much larger than in previous Lanczos diagonalization studies, we determine the boundary of phase separation, the Luttinger Liquid correlation exponent K_rho, and other correlation functions in this region. We find that phase separation occurs for V significantly smaller than previously reported. In addition, for negative U, we find that a uniform state re-enters from phase separation as the electron density is increased towards half filling. For V < V_PS, our results show that superconducting fluctuations are not dominant. The system behaves asymptotically as a Luttinger Liquid with K_rho < 1, but we also find strong low-energy (but gapped) charge-density fluctuations at a momentum not expected for a standard Luttinger Liquid.

Quantum Monte Carlo in the Interaction Representation --- Application to a Spin-Peierls Model

A quantum Monte Carlo algorithm is constructed starting from the standard perturbation expansion in the interaction representation. The resulting configuration space is strongly related to that of the Stochastic Series Expansion (SSE) method, which is based on a direct power series expansion of exp(-beta*H). Sampling procedures previously developed for the SSE method can therefore be used also in the interaction representation formulation. The new method is first tested on the S=1/2 Heisenberg chain. Then, as an application to a model of great current interest, a Heisenberg chain including phonon degrees of freedom is studied. Einstein phonons are coupled to the spins via a linear modulation of the nearest-neighbor exchange. The simulation algorithm is implemented in the phonon occupation number basis, without Hilbert space truncations, and is exact. Results are presented for the magnetic properties of the system in a wide temperature regime, including the T-->0 limit where the chain undergoes a spin-Peierls transition. Some aspects of the phonon dynamics are also discussed. The results suggest that the effects of dynamic phonons in spin-Peierls compounds such as GeCuO3 and NaV2O5 must be included in order to obtain a correct quantitative description of their magnetic properties, both above and below the dimerization temperature.

Steady flow dynamics during granular impact.

We study experimentally and computationally the dynamics of granular flow during impacts where intruders strike a collection of disks from above. In the regime where granular force dynamics are much more rapid than the intruder motion, we find that the particle flow near the intruder is proportional to the instantaneous intruder speed; it is essentially constant when normalized by that speed. The granular flow is nearly divergence free and remains in balance with the intruder, despite the latter's rapid deceleration. Simulations indicate that this observation is insensitive to grain properties, which can be explained by the separation of time scales between intergrain force dynamics and intruder dynamics. Assuming there is a comparable separation of time scales, we expect that our results are applicable to a broad class of dynamic or transient granular flows. Our results suggest that descriptions of static-in-time granular flows might be extended or modified to describe these dynamic flows. Additionally, we find that accurate grain-grain interactions are not necessary to correctly capture the granular flow in this regime.

Orientation, Flow, and Clogging in a Two-Dimensional Hopper: Ellipses vs. Disks

Two-dimensional (2D) hopper flow of disks has been extensively studied. Here, we investigate hopper flow of ellipses with aspect ratio $\alpha = 2$, and we contrast that behavior to the flow of disks. We use a quasi-2D hopper containing photoelastic particles to obtain stress/force information. We simultaneously measure the particle motion and stress. We determine several properties, including discharge rates, jamming probabilities, and the number of particles in clogging arches. For both particle types, the size of the opening, $D$, relative to the size of particles, $\ell$ is an important dimensionless measure. The orientation of the ellipses plays an important role in flow rheology and clogging. The alignment of contacting ellipses enhances the probability of forming stable arches. This study offers insight for applications involving the flow of granular materials consisting of ellipsoidal shapes, and possibly other non-spherical shapes.

Modeling non-harmonic behavior of materials from experimental inelastic neutron scattering and thermal expansion measurements

Based on thermodynamic principles, we derive expressions quantifying the non-harmonic vibrational behavior of materials, which are rigorous yet easily evaluated from experimentally available data for the thermal expansion coefficient and the phonon density of states. These experimentally- derived quantities are valuable to benchmark first-principles theoretical predictions of harmonic and non-harmonic thermal behaviors using perturbation theory, ab initio molecular-dynamics, or Monte-Carlo simulations. We illustrate this analysis by computing the harmonic, dilational, and anharmonic contributions to the entropy, internal energy, and free energy of elemental aluminum and the ordered compound FeSi over a wide range of temperature. Results agree well with previous data in the literature and provide an efficient approach to estimate anharmonic effects in materials.

Intrinsic dielectric response in ferroelectric nano-capacitors

Measurements on 'free-standing' single-crystal barium titanate capacitors with thickness down to 75 nm show a dielectric response typical of large single crystals, rather than conventional thin films. There is a notable absence of any broadening or temperature shift of the dielectric peak or loss tangent. Peak dielectric constants of similar to25 000 are observed, and Curie-Weiss analysis demonstrates first order transformation behaviour. This is in dramatic contrast to results on conventionally deposited thin film capacitor heterostructures, which show large dielectric peak broadening and temperature shifts (e.g. Parker et al 2002 Appl. Phys. Lett. 81 340), as well as an apparent change in the nature-of the paraelectric-ferroelectric transition from first to second order. Our data are compatible with a recent model by Bratkovsky and Levanyuk (2004 Preprint cond-mat/0402100), which attributes dielectric peak broadening to gradient terms that will exist in any thin film capacitor heterostructure. The observed recovery of first order transformation behaviour is consistent with the absence of significant substrate clamping in our experiment, as modelled by Pertsev et al (1998,Phys. Rev. Lett. 80 1988), and illustrates that the second order behaviour seen in conventionally deposited thin films cannot be attributed to the effects of reduced dimensionality in the system, nor to the influence of an intrinsic universal interfacial capacitance associated with the electrode- ferroelectric interface.

Element-specific depth profile of magnetism and stoichiometry at the La0.67Sr0.33MnO3/BiFeO3 interface

Depth-sensitive magnetic, structural and chemical characterization is important in the understanding and optimization of novel physical phenomena emerging at interfaces of transition metal oxide heterostructures. In a simultaneous approach we have used polarized neutron and resonant X-ray reflectometry to determine the magnetic profile across atomically sharp interfaces of ferromagnetic La0.67Sr0.33MnO3 / multiferroic BiFeO3 bi-layers with sub-nanometer resolution. In particular, the X-ray resonant magnetic reflectivity measurements at the Fe and Mn resonance edges allowed us to determine the element specific depth profile of the ferromagnetic moments in both the La0.67Sr0.33MnO3 and BiFeO3 layers. Our measurements indicate a magnetically diluted interface layer within the La0.67Sr0.33MnO3 layer, in contrast to previous observations on inversely deposited layers. Additional resonant X-ray reflection measurements indicate a region of an altered Mn- and O-content at the interface, with a thickness matching that of the magnetic diluted layer, as origin of the reduction of the magnetic moment.

Bright vortex solitons in bose condensates

We suggest the possibility of observing and studying bright vortex solitons in attractive Bose-Einstein condensates in three dimensions with a radial trap. Such systems lie on the verge of critical stability and we discuss the conditions of their stability. We study the interaction between two such solitons. Unlike the text-book solitons in one dimension, the interaction between two radially trapped and axially free three-dimensional solitons is inelastic in nature and involves exchange of particles and deformation in shape. The interaction remains repulsive for all phase δ between them except for δ ≈ 0.