Numerical test of the Cardy-Jacobsen conjecture in the site-diluted Potts model in three dimensions

We present a microcanonical Monte Carlo simulation of the site-diluted Potts model in three dimensions with eight internal states, partly carried out on the citizen supercomputer Ibercivis. Upon dilution, the pure model’s first-order transition becomes of the second order at a tricritical point. We compute accurately the critical exponents at the tricritical point. As expected from the Cardy-Jacobsen conjecture, they are compatible with their random field Ising model counterpart. The conclusion is further reinforced by comparison with older data for the Potts model with four states.

On the well-posedness of a mathematical model for lithium-ion batteries.

We discuss the well-posedness of a mathematical model that is used in the literature for the simulation of lithium-ion batteries. First, a mathematical model based on a macrohomogeneous approach is presented, following previous work. Then it is shown, from a physical and a mathematical point of view, that a boundary condition widely used in the literature is not correct. Although the errors could be just sign typos (which can be explained as carelessness in the use of d/dx versus d/dn, with n the outward unit vector) and authors using this model probably use the correct boundary condition when they solve it in order to do simulations, readers should be aware of the right choice. Therefore, the deduction of the correct boundary condition is done here, and a mathematical study of the well-posedness of the corresponding problem is presented.

Search for very high energy gamma radiation from the radio bright region DR4 of the SNR G78.2+2.1

Data from the HEGRA air shower array are used to set an upper limit on the emission of gamma-radiation above 25 (18) TeV from the direction of the radio bright region DR4 within the SNR G78.2 + 2.1 of 2.5 (7.1). 10^-13 cm^-2 sec^-1. The shock front of SNR G78.2 + 2.1 probably recently overtook the molecular cloud Gong 8 which then acts as a target for the cosmic rays produced within the SNR, thus leading to the expectation of enhanced gamma-radiation. Using a model of Drury, Aharonian and Völk which assumes that SNRs are the sources of galactic cosmic rays via first order Fermi acceleration, we calculated a theoretical prediction for the gamma-ray flux from the DR4 region and compared it with our experimental flux limit. Our 'best estimate' value for the predicted flux lies a factor of about 18 above the upper limit for gamma-ray energies above 25 TeV. Possible reasons for this discrepancy are discussed.

Phase diagram of a polydisperse soft-spheres model for liquids and colloids

The phase diagram of soft spheres with size dispersion is studied by means of an optimized Monte Carlo algorithm which allows us to equilibrate below the kinetic glass transition for all size distributions. The system ubiquitously undergoes a first-order freezing transition. While for a small size dispersion the frozen phase has a crystalline structure, large density inhomogeneities appear in the highly disperse systems. Studying the interplay between the equilibrium phase diagram and the kinetic glass transition, we argue that the experimentally found terminal polydispersity of colloids is a purely kinetic phenomenon.

Development of a Conceptual Site Model for Evaluating Seawater Intrusion under Sea Level Rise Scenarios Using Analytical Methods

Island County is located in the Puget Sound of Washington State and includes several islands, the largest of which is Whidbey Island. Central Whidbey Island was chosen as the project site, as residents use groundwater for their water supply and seawater intrusion near the coast is known to contaminate this resource. In 1989, Island County adopted a Saltwater Intrusion Policy and used chloride concentrations in existing wells in order to define and map “risk zones.” In 2005, this method of defining vulnerability was updated with the use of water level elevations in conjunction with chloride concentrations. The result of this work was a revised map of seawater intrusion vulnerability that is currently in use by Island County. This groundwater management strategy is defined as trigger-level management and is largely a reactive tool. In order to evaluate trends in the hydrogeologic processes at the site, including seawater intrusion under sea level rise scenarios, this report presents a workflow where groundwater flow and discharge to the sea are quantified using a revised conceptual site model. The revised conceptual site model used several simplifying assumptions that allow for first-order quantitative predictions of seawater intrusion using analytical methods. Data from water well reports included lithologic and well construction information, static water levels, and aquifer tests for specific capacity. Results from specific capacity tests define the relationship between discharge and drawdown and were input for a modified Theis equation to solve for transmissivity (Arihood, 2009). Components of the conceptual site model were created in ArcGIS and included interpolation of water level elevation, creation of groundwater basins, and the calculation of net recharge and groundwater discharge for each basin. The revised conceptual site model was then used to hypothesize regarding hydrogeologic processes based on observed trends in groundwater flow. Hypotheses used to explain a reduction in aquifer thickness and hydraulic gradient were: (1) A large increase in transmissivity occurring near the coast. (2) The reduced aquifer thickness and hydraulic gradient were the result of seawater intrusion. (3) Data used to create the conceptual site model were insufficient to resolve trends in groundwater flow. For Hypothesis 2, analytical solutions for groundwater flow under Dupuit assumptions were applied in order to evaluate seawater intrusion under projected sea level rise scenarios. Results indicated that a rise in sea level has little impact on the position of a saltwater wedge; however, a reduction in recharge has significant consequences. Future work should evaluate groundwater flow using an expanded monitoring well network and aquifer recharge should be promoted by reducing surface water runoff.

Piggery pond sludge as a nitrogen source for crops. 1. Mineral N supply estimated from laboratory incubations and field application of stockpiled and wet sludge

Piggery pond sludge (PPS) was applied, as-collected (Wet PPS) and following stockpiling for 12 months ( Stockpiled PPS), to a sandy Sodosol and clay Vertosol at sites on the Darling Downs of Queensland. Laboratory measures of N availability were carried out on unamended and PPS-amended soils to investigate their value in estimating supplementary N needs of crops in Australia's northern grains region. Cumulative net N mineralised from the long-term ( 30 weeks) leached aerobic incubation was described by a first-order single exponential model. The mineralisation rate constant (0.057/week) was not significantly different between Control and PPS treatments or across soil types, when the amounts of initial mineral N applied in PPS treatments were excluded. Potentially mineralisable N (N-o) was significantly increased by the application of Wet PPS, and increased with increasing rate of application. Application of Wet PPS significantly increased the total amount of inorganic N leached compared with the Control treatments. Mineral N applied in Wet PPS contributed as much to the total mineral N status of the soil as did that which mineralised over time from organic N. Rates of CO2 evolution during 30 weeks of aerobic leached incubation indicated that the Stockpiled PPS was more stabilised (19-28% of applied organic C mineralised) than the Wet PPS (35-58% of applied organic C mineralised), due to higher lignin content in the former. Net nitrate-N produced following 12 weeks of aerobic non-leached incubation was highly correlated with net nitrate-N leached during 12 weeks of aerobic incubation (R-2 = 0.96), although it was

Kinetic studies of the metathesis of 1-hexene using Re2O7/γ-Al2O3 and the synthesis and metathesis of oxazolines

The kinetics of the metathesis of 1-hexene using Re2O7/-Al_2O_3 as the catalyst were investigated under a variety of conditions. The experiments were carried out under high vacuum conditions. The product solutions were characterised by gas liquid chromatography and mass spectroscopy. The initial kinetics of the metathesis of 1-hexene showed that the reaction was first order in the weight of the catalyst and second order in the concentration of 1-hexene. A kinetic scheme which correlated the experimental data with the metallocarbene chain mechanism postulated by Herisson and Chauvin and the kinetics of the reaction was explained using a model based on the Langmuir-Hinshelwood theory. The low conversion of 1-hexene to its products is due to termination reactions which most likely occur by the decomposition of the metallocyclobutane intermediate to produce a cyclopropane derivative and an inactive centre. The optimum temperature for the metathesis of 1-hexene over Re_2O_7/-Al2O3 is 45oC and above this temperature, the rate of metathesis decreases rapidly. Co-catalysts alter the active sites for metathesis so that the catalyst is more selective to the metathesis of 1-hexene. However, the regeneration of metathesis activity is much worse for promoted catalysts than for the unpromoted. The synthesis and metathesis of 4,4-dimethyl-2-allowbreak (9-decenyl)-1,3-oxazoline and 4,4-dimethyl-2-allowbreak (3-pentenyl)-1,3-oxazoline was attempted and the products were analysed by thin layer chromatography, infra-red, 13C and 1H nmr and mass spectroscopy. Obtaining the oxazolines in a good yield with high purity was difficult and consequently metathesis of the impure products did not occur.

Contrast integration over area is extensive: a three-stage model of spatial summation

Classical studies of area summation measure contrast detection thresholds as a function of grating diameter. Unfortunately, (i) this approach is compromised by retinal inhomogeneity and (ii) it potentially confounds summation of signal with summation of internal noise. The Swiss cheese stimulus of T. S. Meese and R. J. Summers (2007) and the closely related Battenberg stimulus of T. S. Meese (2010) were designed to avoid these problems by keeping target diameter constant and modulating interdigitated checks of first-order carrier contrast within the stimulus region. This approach has revealed a contrast integration process with greater potency than the classical model of spatial probability summation. Here, we used Swiss cheese stimuli to investigate the spatial limits of contrast integration over a range of carrier frequencies (1–16 c/deg) and raised plaid modulator frequencies (0.25–32 cycles/check). Subthreshold summation for interdigitated carrier pairs remained strong (~4 to 6 dB) up to 4 to 8 cycles/check. Our computational analysis of these results implied linear signal combination (following square-law transduction) over either (i) 12 carrier cycles or more or (ii) 1.27 deg or more. Our model has three stages of summation: short-range summation within linear receptive fields, medium-range integration to compute contrast energy for multiple patches of the image, and long-range pooling of the contrast integrators by probability summation. Our analysis legitimizes the inclusion of widespread integration of signal (and noise) within hierarchical image processing models. It also confirms the individual differences in the spatial extent of integration that emerge from our approach.

Transmission comparison of ultra-long Raman fibre laser based amplification with first and dual order Raman amplification using 10×118 Gbit/s DP-QPSK

Experimental investigations of 10×118 Gbit/s DP-QPSK WDM transmission using three types of distributed Raman amplification techniques are presented. Novel ultra-long Raman fibre laser based amplification with second order counter-propagated pumping is compared with conventional first order and dual order counter-pumped Raman amplification. We demonstrate that URFL based amplification can extend the transmission reach up to a distance of 7520 km in comparison with 5010 km and 6180 km using first order and dual order Raman amplification respectively. © 2014 IEEE.

Linear binocular combination of responses to contrast modulation:contrast-weighted summation in first- and second-order vision

Binocular combination for first-order (luminancedefined) stimuli has been widely studied, but we know rather little about this binocular process for spatial modulations of contrast (second-order stimuli). We used phase-matching and amplitude-matching tasks to assess binocular combination of second-order phase and modulation depth simultaneously. With fixed modulation in one eye, we found that binocularly perceived phase was shifted, and perceived amplitude increased almost linearly as modulation depth in the other eye increased. At larger disparities, the phase shift was larger and the amplitude change was smaller. The degree of interocular correlation of the carriers had no influence. These results can be explained by an initial extraction of the contrast envelopes before binocular combination (consistent with the lack of dependence on carrier correlation) followed by a weighted linear summation of second-order modulations in which the weights (gains) for each eye are driven by the first-order carrier contrasts as previously found for first-order binocular combination. Perceived modulation depth fell markedly with increasing phase disparity unlike previous findings that perceived first-order contrast was almost independent of phase disparity. We present a simple revision to a widely used interocular gain-control theory that unifies first- and second-order binocular summation with a single principle-contrast-weighted summation-and we further elaborate the model for first-order combination. Conclusion: Second-order combination is controlled by first-order contrast.

Detecting Precipitation Climate Changes: An Approach Based on a Stochastic Daily Precipitation Model

2002 Mathematics Subject Classification: 62M10.

Langevin dynamics, large deviations and instantons for the quasi-geostrophic model and two-dimensional Euler equations

We investigate a class of simple models for Langevin dynamics of turbulent flows, including the one-layer quasi-geostrophic equation and the two-dimensional Euler equations. Starting from a path integral representation of the transition probability, we compute the most probable fluctuation paths from one attractor to any state within its basin of attraction. We prove that such fluctuation paths are the time reversed trajectories of the relaxation paths for a corresponding dual dynamics, which are also within the framework of quasi-geostrophic Langevin dynamics. Cases with or without detailed balance are studied. We discuss a specific example for which the stationary measure displays either a second order (continuous) or a first order (discontinuous) phase transition and a tricritical point. In situations where a first order phase transition is observed, the dynamics are bistable. Then, the transition paths between two coexisting attractors are instantons (fluctuation paths from an attractor to a saddle), which are related to the relaxation paths of the corresponding dual dynamics. For this example, we show how one can analytically determine the instantons and compute the transition probabilities for rare transitions between two attractors.

Development of a high-frequency surface-wave radar cross section model for icebergs

In this thesis, the first-order radar cross section (RCS) of an iceberg is derived and simulated. This analysis takes place in the context of a monostatic high frequency surface wave radar with a vertical dipole source that is driven by a pulsed waveform. The starting point of this work is a general electric field equation derived previ- ously for an arbitrarily shaped iceberg region surrounded by an ocean surface. The condition of monostatic backscatter is applied to this general field equation and the resulting expression is inverse Fourier transformed. In the time domain the excitation current of the transmit antenna is specified to be a pulsed sinusoid signal. The result- ing electric field equation is simplified and its physical significance is assessed. The field equation is then further simplified by restricting the iceberg's size to fit within a single radar patch width. The power received by the radar is calculated using this electric field equation. Comparing the received power with the radar range equation gives a general expression for the iceberg RCS. The iceberg RCS equation is found to depend on several parameters including the geometry of the iceberg, the radar frequency, and the electrical parameters of both the iceberg and the ocean surface. The RCS is rewritten in a form suitable for simulations and simulations are carried out for rectangularly shaped icebergs. Simulation results are discussed and are found to be consistent with existing research.

Numerical study of the enlarged O(5) symmetry of the 3D antiferromagnetic RP^(2) spin model

We investigate by means of Monte Carlo simulation and finite-size scaling analysis the critical properties of the three dimensional O (5) non-linear σ model and of the antiferromagnetic RP^(2) model, both of them regularized on a lattice. High accuracy estimates are obtained for the critical exponents, universal dimensionless quantities and critical couplings. It is concluded that both models belong to the same universality class, provided that rather non-standard identifications are made for the momentum-space propagator of the RP^(2) model. We have also investigated the phase diagram of the RP^(2) model extended by a second-neighbor interaction. A rich phase diagram is found, where most of the phase transitions are of the first order.