Regular Tunnelling Sequences in Mixed Systems

We show that the pattern of tunnelling rates can display a vivid and regular pattern when the classical dynamics is of mixed chaotic/regular type. We consider the situation in which the dominant tunnelling route connects to a stable periodic orbit and this orbit is surrounded by a regular island which supports a number of quantum states. We derive an explicit semiclassical expression for the positions and tunnelling rates of these states by use of a complexified trace formula.Comment: submitted to Physica E as a contribution to the workshop proceedings of "Dynamics of Complex Systems" held at the Max Planck Institute for the Physics of Complex Systems in Dresden from March 30 to June 15, 199

Scarring and the statistics of tunnelling

We show that the statistics of tunnelling can be dramatically affected by scarring and derive distributions quantifying this effect. Strong deviations from the prediction of random matrix theory can be explained quantitatively by modifying the Gaussian distribution which describes wavefunction statistics. The modified distribution depends on classical parameters which are determined completely by linearised dynamics around a periodic orbit. This distribution generalises the scarring theory of Kaplan [Phys. Rev. Lett. {\bf 80}, 2582 (1998)] to describe the statistics of the components of the wavefunction in a complete basis, rather than overlaps with single Gaussian wavepackets. In particular it is shown that correlations in the components of the wavefunction are present, which can strongly influence tunnelling-rate statistics. The resulting distribution for tunnelling rates is tested successfully on a two-dimensional double-well potential.Comment: 20 pages, 4 figures, submitted to Ann. Phy

A Matrix Element for Chaotic Tunnelling Rates and Scarring Intensities

It is shown that tunnelling splittings in ergodic double wells and resonant widths in ergodic metastable wells can be approximated as easily-calculated matrix elements involving the wavefunction in the neighbourhood of a certain real orbit. This orbit is a continuation of the complex orbit which crosses the barrier with minimum imaginary action. The matrix element is computed by integrating across the orbit in a surface of section representation, and uses only the wavefunction in the allowed region and the stability properties of the orbit. When the real orbit is periodic, the matrix element is a natural measure of the degree of scarring of the wavefunction. This scarring measure is canonically invariant and independent of the choice of surface of section, within semiclassical error. The result can alternatively be interpretated as the autocorrelation function of the state with respect to a transfer operator which quantises a certain complex surface of section mapping. The formula provides an efficient numerical method to compute tunnelling rates while avoiding the need for the exceedingly precise diagonalisation endemic to numerical tunnelling calculations.Comment: Submitted to Annals of Physics. This work has been submitted to Academic Press for possible publicatio

A simple hourly wind power simulation for the South-West region of Western Australia using MERRA data

A simple simulator capable of generating synthetic hourly values of wind power was developed for the South West region of Western Australia. The global Modern Era Retrospective Analysis for Research and Applications (MERRA) atmospheric database was used to calibrate the simulation with wind speeds 50m above ground level. Analysis of the MERRA data indicated that the normalised residual of hourly wind speed had a double exponential distribution. A translated square-root transformation function yn=(√(1.96+ ye )−1.4)/0.302 was used to convert this to a normal-like distribution so that autoregressive (AR) time series analysis could be used. There was a significant dependency in this time series on the last three hours, so a third order AR model was used to generate hourly 50m wind speed residuals. The MERRA daily average 50m wind speed was found to have a Weibull-like distribution, so a square root conversion was used on the data to obtain a normal distribution. The time series for this distribution was found to have a significant dependency on the values for the last two days, so a second order AR model was also used in the simulation to generate synthetic time series values for the square root of the daily average wind speed. Seasonal, daily, diurnal, and hourly components were added to generate synthetic time series values of total 50m wind speed. To scale this wind speed to turbine hub height, a time varying wind shear factor model was created and calibrated using measured data at a coastal and an inland site. Standard wind turbine power curves were modified to produce an estimate of wind farm power output from the hub-height wind speed. Comparison with measured grid supervisory control and data acquisition (SCADA) data indicated that the simulation generated conservative power output values. The simulation was compared to two other models: a Weibull distribution model, and an AR model with normally distributed residuals. The statistical fit with the SCADA data was found to be closer than these two models. Spatial correlation using only the MERRA data was found to be higher than the SCADA data, indicating that there is still a further source of variability to be accounted for. Hence the simulation spatial correlation was calibrated to previously reported findings, which were similar to the SCADA data

Synchrotron radiation study of the relation between structure and strain in polyurethane elastomers

This paper describes a system for the study of the relation between structure and applied strain in thermoplastic polyurethane elastomers using the Australian National Beamline Facility at the Photon Factory, KEK, Tsukuba, Japan. The system uses the sagittal focusing monochromator at beamline 20B to provide a high-intensity focused beam which then falls on the specimen mounted in a miniature tensometer mounted in the unique vacuum diffractometer (BIGDIFF). Imaging plates were used to record simultaneously SAXS and WAXS patterns from the specimen at a particular strain. The change in SAXS and WAXS patterns with loading and unloading was recorded using a ten-plate imaging-plate changer

A large-scale renewable electricity supply system by 2030: Solar, wind, energy efficiency, storage and inertia for the South West Interconnected System (SWIS) in Western Australia

An interactive web tool was created to simulate 100% renewable electricity supply scenarios for the South-West Interconnected System (SWIS) in the south-west of Western Australia. The SWIS is isolated from other grids and currently has no available hydropower. Hence it makes a good case study of how supply and demand might be balanced on an hour-by-hour basis and grid stability maintained without the benefit of energy import/export or pumped hydroelectric storage. The tool included regional models for wind and solar power, so that hypothetical power stations were not confined to sites with existing wind farms or solar power stations, or sites with measurements of wind speed and solar radiation. A generic model for solar thermal storage and simple models for energy efficiency, distributed battery storage and power to gas storage were also developed. Due to the urgency of climate change mitigation a rapid construction schedule of completion by 2030, rather than the more common target of 2050, was set. A scenario with high wind generation, and scenarios with varying levels of solar power, wind power, distributed battery storage, energy efficiency improvements and power to gas systems were considered. The battery storage system and PV arrays were configured to provide synthetic inertia to maintain grid stability (with a small loss in capacity for each), and existing synchronous generators were kept spinning with no fuel input, adding a small increase to the electrical load demand. The level of synthetic inertia provided by battery storage was estimated for each scenario. The results indicated that a balanced mix of solar PV, solar thermal, efficiency, and storage were the most feasible to be built on a rapid time scale. The required capacity and build rate of the generation and storage systems would be reduced if energy efficiency improvements were implemented on a more rapid schedule compared to the current global improvement rate. The scenario with very high levels of wind power (∼80% generation) were found to be capable of meeting SWIS reliability criteria if very large amounts of distributed storage or some high capacity seasonal reserve generation system such as power to gas were present. High levels of battery storage capacity and efficiency improvement could be as effective as a power to gas system. It was confirmed that all scenarios provided the same or greater levels of inertia than presently provided by conventional generators. This tool showed that it is possible to examine renewable energy scenarios for regional electricity networks without high computing power

A realistic example of chaotic tunneling: The hydrogen atom in parallel static electric and magnetic fields

Statistics of tunneling rates in the presence of chaotic classical dynamics is discussed on a realistic example: a hydrogen atom placed in parallel uniform static electric and magnetic fields, where tunneling is followed by ionization along the fields direction. Depending on the magnetic quantum number, one may observe either a standard Porter-Thomas distribution of tunneling rates or, for strong scarring by a periodic orbit parallel to the external fields, strong deviations from it. For the latter case, a simple model based on random matrix theory gives the correct distribution.Comment: Submitted to Phys. Rev.

Gallavotti-Cohen theorem, Chaotic Hypothesis and the zero-noise limit

The Fluctuation Relation for a stationary state, kept at constant energy by a deterministic thermostat - the Gallavotti-Cohen Theorem -- relies on the ergodic properties of the system considered. We show that when perturbed by an energy-conserving random noise, the relation follows trivially for any system at finite noise amplitude. The time needed to achieve stationarity may stay finite as the noise tends to zero, or it may diverge. In the former case the Gallavotti-Cohen result is recovered, while in the latter case, the crossover time may be computed from the action of `instanton' orbits that bridge attractors and repellors. We suggest that the `Chaotic Hypothesis' of Gallavotti can thus be reformulated as a matter of stochastic stability of the measure in trajectory space. In this form this hypothesis may be directly tested

Signatures of unstable semiclassical trajectories in tunneling

It was found recently that processes of multidimensional tunneling are generally described at high energies by unstable semiclassical trajectories. We study two observational signatures related to the instability of trajectories. First, we find an additional power-law dependence of the tunneling probability on the semiclassical parameter as compared to the standard case of potential tunneling. The second signature is substantial widening of the probability distribution over final-state quantum numbers. These effects are studied using modified semiclassical technique which incorporates stabilization of the tunneling trajectories. The technique is derived from first principles. We obtain expressions for the inclusive and exclusive tunneling probabilities in the case of unstable semiclassical trajectories. We also investigate the "phase transition" between the cases of stable and unstable trajectories across certain "critical" value of energy. Finally, we derive the relation between the semiclassical probabilities of tunneling from the low-lying and highly excited initial states. This puts on firm ground a conjecture made previously in the semiclassical description of collision-induced tunneling in field theory.Comment: Journal version; 48 pages, 16 figure