Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.

Libros de texto de matemáticas en España en los siglos XVIII y XIX

A finales del siglo XVIII, en Europa el conocimiento científico se había desarrollado extraordinariamente. Surgen los nombres de Lavoisier, Ritcher, Coulomb y Celsius entre otros muchos. Se enuncian leyes en química y física; junto a ellas también florece la matemática de la mano de Euler, Lagrange, D«Alambert, Monge, por citar sólo unos cuantos. Mientras tanto, el atraso de las matemáticas españolas se debía, entre otras causas, al pobre estado en que se encontraban las universidades: aún de tipo medieval y de carácter eclesiástico. Esto lo evidencia Fray Benito Jerónimo Feijoo en la carta titulada Causas del atraso que se padece en España en orden a las ciencias naturales, y el Marqués de la Ensenada quien, en 1748, se lo expresa al rey Fernando VI. Las deficiencias de las universidades tenían que ver con la enseñanza memorística, textos anticuados e interés primordial por disciplinas como derecho, teología y filosofía en detrimento de las matemáticas y las ciencias.

Singular stresses at the tip of a crack terminating at frictional interface of fibre reinforced composites

The stress singularities at the tip of a crack that terminates at a frictional interface between two layers in anisotropic composites are investigated. The order of stress singularities is determined by solving the characteristic equations obtained from the boundary conditions and the frictional interface conditions for the cases concerned. The interface is assumed to be governed by Coulomb's law of friction. Numerical results are presented for the cases with a crack terminating at a frictional interface of a fibre reinforced composite, and it is shown that there is a big difference of stress singularities between cases with and without considering friction along the interface.

Single-photon tunneling via localized surface plasmons

Strong evidence of a single-photon tunneling effect, a direct analog of single-electron tunneling, has been obtained in the measurements of light tunneling through individual subwavelength pinholes in a gold film covered with a layer of polydiacetylene. The transmission of some pinholes reached saturation because of the optical nonlinearity of polydiacetylene at a very low light intensity of a few thousand photons per second. This result is explained theoretically in terms of a "photon blockade," similar to the Coulomb blockade phenomenon observed in single-electron tunneling experiments. Single-photon tunneling may find applications in the fields of quantum communication and information processing.

High harmonic generation by halogen anions and noble gas atoms in a laser field

A comparative study of high harmonic generation (HHG) by atoms and ions with active p-electrons is carried out in the theoretical framework of the rescattering mechanism. The substate with m(l) = 0, i.e. zero orbital momentum projection along the electric vector of a linearly polarized laser wave, is found to give the major contribution to the HHG rate. Our calculations for HHG by an H atom in an excited 2p-state demonstrate that the rate for recombination into a final state with a different value of m(l) (= +/- 1), is higher for lower harmonic orders N, while for higher N (beyond the plateau domain) the difference vanishes. For species with closed electron shells, the m(l)-changing transitions are forbidden by the Pauli exclusion principle. We report absolute HHG rates for halogen ions and noble gas atoms at various intensities. These results demonstrate that the Coulomb binding potential of the atoms considerably enhances both the ionization and recombination steps in the rescattering process. However, the weak binding energy of the anions allows lower orders of HHG to be efficiently produced at relatively low intensities, from which we conclude that observation of HHG by an anion is experimentally feasible.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

Angular distribution for the elastic scattering of electrons from Ar+(3s^2 3p^5 ^2P) above the first inelastic threshold

The measured angular differential cross section (DCS) for the elastic scattering of electrons from Ar+(3s2 3p5 2P) at the collision energy of 16 eV is presented. By solving the Hartree-Fock equations, we calculate the corresponding theoretical DCS including the coupling between the orbital angular momenta and spin of the incident electron and those of the target ion and also relaxation effects. Since the collision energy is above one inelastic threshold for the transition 3s2 3p5 2P–3s 3p6 2S, we consider the effects on the DCS of inelastic absorption processes and elastic resonances. The measurements deviate significantly from the Rutherford cross section over the full angular range observed, especially in the region of a deep minimum centered at approximately 75°. Our theory and an uncoupled, unrelaxed method using a local, spherically symmetric potential by Manson [Phys. Rev. 182, 97 (1969)] both reproduce the overall shape of the measured DCS, although the coupled Hartree-Fock approach describes the depth of the minimum more accurately. The minimum is shallower in the present theory owing to our lower average value for the d-wave non-Coulomb phase shift s2, which is due to the high sensitivity of s2 to the different scattering potentials used in the two models. The present measurements and calculations therefore show the importance of including coupling and relaxation effects when accurately modeling electron-ion collisions. The phase shifts obtained by fitting to the measurements are compared with the values of Manson and the present method.

Ion Association in [bmim][PF6]/Naphthalene Mixtures: An Experimental and Computational Study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity a of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction (x), whereas the conductivity per ion displays a clear peak at x approximate to 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x <0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.

Latin American migration in the central area of Iquique: New fronts of residential location and unequal forms of access to housing

Historically the central area of the city of Iquique has been established as residential space migrants choosing from different backgrounds , however since the late 2000s migration flows are diversified being mostly Latin American immigrants who live in precarious conditions , accessing tugurizados properties , deteriorated in an increasingly growing informal market. The results presented here are derived from quantitative residential location of migrants , as well as the implementation of 13 in-depth interviews . From these results emerge that Latin American migrants access to the same places where once lived internal migrants, however they inhabit a restrictive market , uneven and inadequate living conditions lease, but allows them to articulate residence and proximity to industrial networks , social and popular trade.

Structure of warm dense matter via angularly resolved x-ray scatter

In this paper we describe experimental results on angularly resolved x-ray scatter from a sample of warm dense aluminium that has been created by double sided laser-driven shock compression. The experiment was carried out on the Central Laser Facility of the Rutherford Appleton Laboratory, using the VULCAN laser. The form of the angularly resolved scatter cross-section was compared with predictions based on a series of molecular dynamics simulations with an embedded atom potential, a Yukakwa potential and a bare Coulomb potential. The importance of screening is evident from the comparison and the embedded atom model seems to match experiment better than the Yukawa potential.

The Perils of Moviegoing in America, 1896-1950

The Perils of Moviegoing in America is a film history that examines the various physical and (perceived) moral dangers facing audiences during the first fifty years of film exhibition.

Chapter 1: “Conflagration”
As early as 1897, a major fire broke out at a film exhibition in San Francisco, with flames burning the projectionist and nearby audience members. From that point until the widespread adoption of safety stock in 1950, fires were a very common movie-going experience. Hundreds of audience members lost their lives in literally thousands of theatre fires, ranging from early nickelodeons to the movie palaces of the thirties and forties.

Chapter 2: “Thieves Among Us”
Bandits robbed movie theatres on hundreds of occasions from the early days of film exhibition through the end of the Great Depression. They held up ticket booths, and they dynamited theatre safes. They also shot theatre managers, ushers, and audience members, as a great many of the robberies occurred while movies were playing on the screens inside.

Chapter 3: “Bombs Away”
Bombings at movie theatres became common in small towns and large cities on literally hundreds of occasions from 1914 to the start of World War II. Some were incendiary bombs, and some were stench bombs; both could be fatal, whether due to explosions or to the trampling of panicked moviegoers

Chapter 4: “It’s Catching”
Widespread movie-going in the early 20th century provoked an outcry from numerous doctors and optometrists who believed that viewing films could do irreparable harm to the vision of audience members. Medical publications (including the Journal of the American Medical Association) published major studies on this perceived problem, which then filtered into popular-audience magazines and newspapers.

Chapter 5: “The Devil’s Apothecary Shops”
Sitting in the dark with complete strangers proved worrisome for many early filmgoers, who had good reason to be concerned. Darkness meant that prostitutes could easily work in the balconies of some movie theatres, as could “mashers” who molested female patrons (and sometimes children) after the lights were dimmed. That was all in addition to the various murderers who used the cover of darkness to commit their crimes at movie theatres.

Chapter 6: “Blue Sundays”
Blue laws were those regulations that prohibited businesses from operating on Sundays. Most communities across the US had such legislation on their books, which by the nickelodeon era were at odds with the thousands of filmgoers who went to the movies every Sunday. Theatre managers were often arrested, making newspaper headlines over and over again. Police sometimes even arrested entire film audiences as accomplices in the Blue Law violations.

Chapter 7: “Something for Nothing”
In an effort to bolster ticket sales, many movie theatres in the 1910s began to hold lotteries in which lucky audience members won cash prizes; by the time of the Great Depression, lotteries like “Bank Night” became a common aspect of the theatre-going enterprise. However, reception studies have generally overlooked the intense (and sometimes coordinated) efforts by police, politicians, and preachers to end this practice, which they viewed as illegal and immoral gambling.

Quasi-classical model of non-destructive wavepacket manipulation by intense ultrashort nonresonant laser pulses

A quasi-classical model (QCM) of nuclear wavepacket generation, modification and imaging by three intense ultrafast near-infrared laser pulses has been developed. Intensities in excess of 10(13) W cm(-2) are studied, the laser radiation is non-resonant and pulse durations are in the few-cycle regime, hence significantly removed from the conditions typical of coherent control and femtochemistry. The 1s sigma ground state of the D-2 precursor is projected onto the available electronic states in D-2(+) (1s sigma(g) ground and 2p sigma(u) dissociative) and D+ + D+ (Coulomb explosion) by tunnel ionization by an ultrashort 'pump' pulse, and relative populations are found numerically. A generalized non-adiabatic treatment allows the dependence of the initial vibrational population distribution on laser intensity to be calculated. The wavepacket is approximated as a classical ensemble of particles moving on the 1s sigma(g) potential energy surface (PES), and hence follow trajectories of different amplitudes and frequencies depending on the initial vibrational state. The 'control' pulse introduces a time-dependent polarization of the molecular orbital, causing the PES to be modified according to the dynamic Stark effect and the transition dipole. The trajectories adjust in amplitude, frequency and phase-offset as work is done on or by the resulting force; comparing the perturbed and unperturbed trajectories allows the final vibrational state populations and phases to be determined. The action of the 'probe' pulse is represented by a discrete internuclear boundary, such that elements of the ensemble at a larger internuclear separation are assumed to be photodissociated. The vibrational populations predicted by the QCM are compared to recent quantum simulations (Niederhausen and Thumm 2008 Phys. Rev. A 77 013404), and a remarkable agreement has been found. The applicability of this model to femtosecond and attosecond time-scale experiments is discussed and the relation to established femtochemistry and coherent control techniques are explored.

The death of massive stars - I. Observational constraints on the progenitors of Type II-P supernovae

We present the results of a 10.5-yr, volume-limited (28-Mpc) search for supernova (SN) progenitor stars. In doing so we compile all SNe discovered within this volume (132, of which 27 per cent are Type Ia) and determine the relative rates of each subtype from literature studies. The core-collapse SNe break down into 59 per cent II-P and 29 per cent Ib/c, with the remainder being IIb (5 per cent), IIn (4 per cent) and II-L (3 per cent). There have been 20 II-P SNe with high-quality optical or near-infrared pre-explosion images that allow a meaningful search for the progenitor stars. In five cases they are clearly red supergiants, one case is unconstrained, two fall on compact coeval star clusters and the other twelve have no progenitor detected. We review and update all the available data for the host galaxies and SN environments (distance, metallicity and extinction) and determine masses and upper mass estimates for these 20 progenitor stars using the STARS stellar evolutionary code and a single consistent homogeneous method. A maximum likelihood calculation suggests that the minimum stellar mass for a Type II-P to form is m(min) = 8.5(-1.5)(+1) M-circle dot and the maximum mass for II-P progenitors is m(max) = 16.5 +/- 1.5 M-circle dot, assuming a Salpeter initial mass function holds for the progenitor population (in the range Gamma = -1.35(-0.7)(+0.3)). The minimum mass is consistent with current estimates for the upper limit to white dwarf progenitor masses, but the maximum mass does not appear consistent with massive star populations in Local Group galaxies. Red supergiants in the Local Group have masses up to 25 M-circle dot and the minimum mass to produce a Wolf-Rayet star in single star evolution (between solar and LMC metallicity) is similarly 25-30 M-circle dot. The reason we have not detected any high-mass red supergiant progenitors above 17 M-circle dot is unclear, but we estimate that it is statistically significant at 2.4 sigma confidence. Two simple reasons for this could be that we have systematically underestimated the progenitor masses due to dust extinction or that stars between 17-25 M-circle dot produce other kinds of SNe which are not II-P. We discuss these possibilities and find that neither provides a satisfactory solution. We term this discrepancy the 'red supergiant problem' and speculate that these stars could have core masses high enough to form black holes and SNe which are too faint to have been detected. We compare the Ni-56 masses ejected in the SNe to the progenitor mass estimates and find that low-luminosity SNe with low Ni-56 production are most likely to arise from explosions of low-mass progenitors near the mass threshold that can produce a core-collapse.

Mesophases in Nearly 2D Room-Temperature Ionic Liquids

Computer simulations of (i) a [C(12)mim][Tf2N] film of nanometric thickness squeezed at kbar pressure by a piecewise parabolic confining potential reveal a mesoscopic in-plane density and composition modulation reminiscent of mesophases seen in 3D samples of the same room-temperature ionic liquid (RTIL). Near 2D confinement, enforced by a high normal load, as well as relatively long aliphatic chains are strictly required for the mesophase formation, as confirmed by computations for two related systems made of (ii) the same [C(12)mim][Tf2N] adsorbed at a neutral solid surface and (iii) a shorter-chain RTIL ([C(4)mim][Tf2N]) trapped in the potential well of part i. No in-plane modulation is seen for ii and iii. In case ii, the optimal arrangement of charge and neutral tails is achieved by layering parallel to the surface, while, in case iii, weaker dispersion and packing interactions are unable to bring aliphatic tails together into mesoscopic islands, against overwhelming entropy and Coulomb forces. The onset of in-plane mesophases could greatly affect the properties of long-chain RTILs used as lubricants.