Independent particle descriptions of tunneling using the many-body quantum transport approach

Currents across thin insulators are commonly taken as single electrons moving across classically forbidden regions; this independent particle picture is well-known to describe most tunneling phenomena. Examining quantum transport from a different perspective, i.e., by explicit treatment of electron-electron interactions, we evaluate different single particle approximations with specific application to tunneling in metal-molecule-metal junctions. We find maximizing the overlap of a Slater determinant composed of single-particle states to the many-body current-carrying state is more important than energy minimization for defining single-particle approximations in a system with open boundary conditions. Thus the most suitable single particle effective potential is not one commonly in use by electronic structure methods, such as the Hartree-Fock or Kohn-Sham approximations.

Shell Model Study of Local and Global Energy Barriers in KDP

We perform a study of the energetics of KH2PO4 (KDP) by using a shell model (SM) which was constructed by adjusting the interaction parameters to ab initio calculations, and was fitted to reproduce phonons, polarization-inversion energies and structural properties. We calculate the energy profiles by performing global displacements and local distortions following the ferroelectric (FE) mode pattern in clusters of different sizes embedded in a paraelectric (PE) phase matrix. These properties are expected to be relevant to the PE-FE phase transition. The obtained SM results are compared to corresponding ab initio (AI) data. The global instabilities are found in good agreement for both KDP and DKDP. We also find qualitative good agreement in the KDP structure and even quantitative agreement in the expanded DKDP structure for the local distortions. The SM results reproduce well different trends like increasing instabilities as the cluster sizes grows, as the heavier atoms are included, and as the volume is increased, in accordance with the corresponding data from AI calculations.

Single ionization of helium by antiproton impact

We study the ionization of helium (fie) in collision with antiprotons (p) in the energy range from 10 keV to 1000 keV. We adopt a semiclassical single center close coupling approach in which the wave function for the electron is expanded in a B-spline basis centered on the nucleus of the He atom, The calculations are performed using two different models: the independent particle (IP) model and the one-electron (OE) approximation. The interaction between the active electron and the rest of the atom, i.e. passive electron and nucleus, is represented by a model potential. The results obtained are compared with experimental data as well as with existing theoretical calculations. (c) 2005 Published by Elsevier B.V.

Statistical theory of finite Fermi systems with chaotic excited eigenstates

A theory of strongly interacting Fermi systems of a few particles is developed. At high excit at ion energies (a few times the single-parti cle level spacing) these systems are characterized by an extreme degree of complexity due to strong mixing of the shell-model-based many-part icle basis st at es by the residual two- body interaction. This regime can be described as many-body quantum chaos. Practically, it occurs when the excitation energy of the system is greater than a few single-particle level spacings near the Fermi energy. Physical examples of such systems are compound nuclei, heavy open shell atoms (e.g. rare earths) and multicharged ions, molecules, clusters and quantum dots in solids. The main quantity of the theory is the strength function which describes spreading of the eigenstates over many-part icle basis states (determinants) constructed using the shell-model orbital basis. A nonlinear equation for the strength function is derived, which enables one to describe the eigenstates without diagonalization of the Hamiltonian matrix. We show how to use this approach to calculate mean orbital occupation numbers and matrix elements between chaotic eigenstates and introduce typically statistical variable s such as t emperature in an isolated microscopic Fermi system of a few particles.

Macroscopic limit of time-dependent density-functional theory for adiabatic local approximations of the exchange-correlation kernel

Time-dependent density-functional theory is a rather accurate and efficient way to compute electronic excitations for finite systems. However, in the macroscopic limit (systems of increasing size), for the usual adiabatic random-phase, local-density, or generalized-gradient approximations, one recovers the Kohn-Sham independent-particle picture, and thus the incorrect band gap. To clarify this trend, we investigate the macroscopic limit of the exchange-correlation kernel in such approximations by means of an algebraical analysis complemented with numerical studies of a one-dimensional tight-binding model. We link the failure to shift the Kohn-Sham spectrum of these approximate kernels to the fact that the corresponding operators in the transition space act only on a finite subspace.

Automated mixed dimensional modelling from 2D and 3D CAD models

The motivation for this paper is to present procedures for automatically creating idealised finite element models from the 3D CAD solid geometry of a component. The procedures produce an accurate and efficient analysis model with little effort on the part of the user. The technique is applicable to thin walled components with local complex features and automatically creates analysis models where 3D elements representing the complex regions in the component are embedded in an efficient shell mesh representing the mid-faces of the thin sheet regions. As the resulting models contain elements of more than one dimension, they are referred to as mixed dimensional models. Although these models are computationally more expensive than some of the idealisation techniques currently employed in industry, they do allow the structural behaviour of the model to be analysed more accurately, which is essential if appropriate design decisions are to be made. Also, using these procedures, analysis models can be created automatically whereas the current idealisation techniques are mostly manual, have long preparation times, and are based on engineering judgement. In the paper the idealisation approach is first applied to 2D models that are used to approximate axisymmetric components for analysis. For these models 2D elements representing the complex regions are embedded in a 1D mesh representing the midline of the cross section of the thin sheet regions. Also discussed is the coupling, which is necessary to link the elements of different dimensionality together. Analysis results from a 3D mixed dimensional model created using the techniques in this paper are compared to those from a stiffened shell model and a 3D solid model to demonstrate the improved accuracy of the new approach. At the end of the paper a quantitative analysis of the reduction in computational cost due to shell meshing thin sheet regions demonstrates that the reduction in degrees of freedom is proportional to the square of the aspect ratio of the region, and for long slender solids, the reduction can be proportional to the aspect ratio of the region if appropriate meshing algorithms are used.

Local-field and excitonic effects in the optical response of alpha-alumina

This paper introduces key ingredients of the dielectric response of a-alumina that go beyond an independent-particle (IP) treatment of the valence-electron excitations. The optical-response functions were calculated from first-principles both at the Bethe-Salpeter and the random-phase approximation (RPA) levels. Excitonic effects obtained within the Bethe-Salpeter framework were found essential for reproducing the low-energy part of the experimental spectra (below 15 eV) and the bound exciton in particular. For higher energies, local-field effects introduced through the RPA modified considerably the IP results and provided a satisfactory account of the reflectivity spectra and of the position and shape of the dominant bulk plasmon resonance in the electron energy-loss spectra.

Nonlinear optics from an ab initio approach by means of the dynamical Berry phase: Application to second- and third-harmonic generation in semiconductors

We present an ab initio real-time-based computational approach to study nonlinear optical properties in condensed matter systems that is especially suitable for crystalline solids and periodic nanostructures. The equations of motion and the coupling of the electrons with the external electric field are derived from the Berry-phase formulation of the dynamical polarization [Souza et al., Phys. Rev. B 69, 085106 (2004)]. Many-body effects are introduced by adding single-particle operators to the independent-particle Hamiltonian. We add a Hartree operator to account for crystal local effects and a scissor operator to correct the independent particle band structure for quasiparticle effects. We also discuss the possibility of accurately treating excitonic effects by adding a screened Hartree-Fock self-energy operator. The approach is validated by calculating the second-harmonic generation of SiC and AlAs bulk semiconductors: an excellent agreement is obtained with existing ab initio calculations from response theory in frequency domain [Luppi et al., Phys. Rev. B 82, 235201 (2010)]. We finally show applications to the second-harmonic generation of CdTe and the third-harmonic generation of Si. 

Understanding macroalgal dispersal in a complex hydrodynamic environment: a combined population genetic and physical modelling approach

Gene flow in macroalgal populations can be strongly influenced by spore or gamete dispersal. This, in turn, is influenced by a convolution of the effects of current flow and specific plant reproductive strategies. Although several studies have demonstrated genetic variability in macroalgal populations over a wide range of spatial scales, the associated current data have generally been poorly resolved spatially and temporally. In this study, we used a combination of population genetic analyses and high-resolution hydrodynamic modelling to investigate potential connectivity between populations of the kelp Laminaria digitata in the Strangford Narrows, a narrow channel characterized by strong currents linking the large semi-enclosed sea lough, Strangford Lough, to the Irish Sea. Levels of genetic structuring based on six microsatellite markers were very low, indicating high levels of gene flow and a pattern of isolation-by-distance, where populations are more likely to exchange migrants with geographically proximal populations, but with occasional long-distance dispersal. This was confirmed by the particle tracking model, which showed that, while the majority of spores settle near the release site, there is potential for dispersal over several kilometres. This combined population genetic and modelling approach suggests that the complex hydrodynamic environment at the entrance to Strangford Lough can facilitate dispersal on a scale exceeding that proposed for L. digitata in particular, and the majority of macroalgae in general. The study demonstrates the potential of integrated physical–biological approaches for the prediction of ecological changes resulting from factors such as anthropogenically induced coastal zone changes.

Positron scattering and annihilation on noble-gas atoms

Positron scattering and annihilation on noble-gas atoms is studied ab initio using many-body theory methods for positron energies below the positronium formation threshold. We show that in this energy range, the many-body theory yields accurate numerical results and provides a near-complete understanding of the positron–noble-gas atom system. It accounts for positron-atom and electron-positron correlations, including the polarization of the atom by the positron and the nonperturbative effect of virtual positronium formation. These correlations have a large influence on the scattering dynamics and result in a strong enhancement of the annihilation rates compared to the independent-particle mean-field description. Computed elastic scattering cross sections are found to be in good agreement with recent experimental results and Kohn variational and convergent close-coupling calculations. The calculated values of the annihilation rate parameter Zeff (effective number of electrons participating in annihilation) rise steeply along the sequence of noble-gas atoms due to the increasing strength of the correlation effects, and agree well with experimental data.

Heating Mechanisms for Intermittent Loops in Active Region Cores from AIA/SDO EUV Observations

We investigate intensity variations and energy deposition in five coronal loops in active region cores. These were selected for their strong variability in the AIA/SDO 94 Å intensity channel. We isolate the hot Fe XVIII and Fe XXI components of the 94 Å and 131 Å by modeling and subtracting the "warm" contributions to the emission. HMI/SDO data allow us to focus on "inter-moss" regions in the loops. The detailed evolution of the inter-moss intensity time series reveals loops that are impulsively heated in a mode compatible with a nanoflare storm, with a spike in the hot 131 Å signals leading and the other five EUV emission channels following in progressive cooling order. A sharp increase in electron temperature tends to follow closely after the hot 131 Å signal confirming the impulsive nature of the process. A cooler process of growing emission measure follows more slowly. The Fourier power spectra of the hot 131 Å signals, when averaged over the five loops, present three scaling regimes with break frequencies near 0.1 min–1 and 0.7 min–1. The low frequency regime corresponds to 1/f noise; the intermediate indicates a persistent scaling process and the high frequencies show white noise. Very similar results are found for the energy dissipation in a 2D "hybrid" shell model of loop magneto-turbulence, based on reduced magnetohydrodynamics, that is compatible with nanoflare statistics. We suggest that such turbulent dissipation is the energy source for our loops

Experimental and discrete element modelling of cohesive iron ore fines

In this study, the behaviour of iron ore fines with varying levels of adhesion was investigated using a confined compression test and a uniaxial test. The uniaxial test was conducted using the semi-automated uniaxial EPT tester in which the cohesive strength of a bulk solid is evaluated from an unconfined compression test following a period of consolidation to a pre-defined vertical stress. The iron ore fines were also tested by measuring both the vertical and circumferential strains on the cylindrical container walls under vertical loading in a separate confined compression tester - the K0 tester, to determine the lateral pressure ratio. Discrete Element Method simulations of both experiments were carried out and the predictions were compared with the experimental observations. A recently developed DEM contact model for cohesive solids, an Elasto-Plastic Adhesive model, was used. This particle contact model uses hysteretic non-linear loading and unloading paths and an adhesion parameter which is a function of the maximum contact overlap. The model parameters for the simulations are phenomenologically based to reproduce the key bulk characteristics exhibited by the solid. The simulation results show a good agreement in capturing the stress history dependent behaviour depicted by the flow function of the cohesive iron ore fines while also providing a reasonably good match for the lateral pressure ratio observed during the confined compression K0 tests. This demonstrates the potential for the DEM model to be used in the simulation of bulk handling applications.

Strong second harmonic generation in SiC, ZnO, GaN two-dimensional hexagonal crystals from first-principles many-body calculations

The second harmonic generation (SHG) intensity spectrum of SiC, ZnO, GaN two-dimensional hexagonal crystals is calculated by using a real-time first-principles approach based on Green's function theory [Attaccalite et al., Phys. Rev. B: Condens. Matter Mater. Phys. 2013 88, 235113]. This approach allows one to go beyond the independent particle description used in standard first-principles nonlinear optics calculations by including quasiparticle corrections (by means of the GW approximation), crystal local field effects and excitonic effects. Our results show that the SHG spectra obtained using the latter approach differ significantly from their independent particle counterparts. In particular they show strong excitonic resonances at which the SHG intensity is about two times stronger than within the independent particle approximation. All the systems studied (whose stabilities have been predicted theoretically) are transparent and at the same time exhibit a remarkable SHG intensity in the range of frequencies at which Ti:sapphire and Nd:YAG lasers operate; thus they can be of interest for nanoscale nonlinear frequency conversion devices. Specifically the SHG intensity at 800 nm (1.55 eV) ranges from about 40-80 pm V(-1) in ZnO and GaN to 0.6 nm V(-1) in SiC. The latter value in particular is 1 order of magnitude larger than values in standard nonlinear crystals.