Dynamics of Barrierless and Activated Chemical Reactions in a Dispersive Medium within the Fractional Diffusion Equation Approach

Barrierless chemical reactions have often been modeled as a Brownian motion on a one-dimensional harmonic potential energy surface with a position-dependent reaction sink or window located near the minimum of the surface. This simple (but highly successful) description leads to a nonexponential survival probability only at small to intermediate times but exponential decay in the long-time limit. However, in several reactive events involving proteins and glasses, the reactions are found to exhibit a strongly nonexponential (power law) decay kinetics even in the long time. In order to address such reactions, here, we introduce a model of barrierless chemical reaction where the motion along the reaction coordinate sustains dispersive diffusion. A complete analytical solution of the model can be obtained only in the frequency domain, but an asymptotic solution is obtained in the limit of long time. In this case, the asymptotic long-time decay of the survival probability is a power law of the Mittag−Leffler functional form. When the barrier height is increased, the decay of the survival probability still remains nonexponential, in contrast to the ordinary Brownian motion case where the rate is given by the Smoluchowski limit of the well-known Kramers' expression. Interestingly, the reaction under dispersive diffusion is shown to exhibit strong dependence on the initial state of the system, thus predicting a strong dependence on the excitation wavelength for photoisomerization reactions in a dispersive medium. The theory also predicts a fractional viscosity dependence of the rate, which is often observed in the reactions occurring in complex environments.

Fractional damping: Stochastic origin and finite approximations

Fractional-order derivatives appear in various engineering applications including models for viscoelastic damping. Damping behavior of materials, if modeled using linear, constant coefficient differential equations, cannot include the long memory that fractional-order derivatives require. However, sufficiently great rnicrostructural disorder can lead, statistically, to macroscopic behavior well approximated by fractional order derivatives. The idea has appeared in the physics literature, but may interest an engineering audience. This idea in turn leads to an infinite-dimensional system without memory; a routine Galerkin projection on that infinite-dimensional system leads to a finite dimensional system of ordinary differential equations (ODEs) (integer order) that matches the fractional-order behavior over user-specifiable, but finite, frequency ranges. For extreme frequencies (small or large), the approximation is poor. This is unavoidable, and users interested in such extremes or in the fundamental aspects of true fractional derivatives must take note of it. However, mismatch in extreme frequencies outside the range of interest for a particular model of a real material may have little engineering impact.

Polymer melt dynamics: Microscopic roots of fractional viscoelasticity

The rheological properties of polymer melts and other complex macromolecular fluids are often successfully modeled by phenomenological constitutive equations containing fractional differential operators. We suggest a molecular basis for such fractional equations in terms of the generalized Langevin equation (GLE) that underlies the renormalized Rouse model developed by Schweizer [J. Chem. Phys. 91, 5802 (1989)]. The GLE describes the dynamics of the segments of a tagged chain under the action of random forces originating in the fast fluctuations of the surrounding polymer matrix. By representing these random forces as fractional Gaussian noise, and transforming the GLE into an equivalent diffusion equation for the density of the tagged chain segments, we obtain an analytical expression for the dynamic shear relaxation modulus G(t), which we then show decays as a power law in time. This power-law relaxation is the root of fractional viscoelastic behavior.

A Model of Anomalous Enzyme-Catalyzed Gel Degradation Kinetics

We show that a model of target location involving n noninteracting particles moving subdiffusively along a line segment (a generalization of a model introduced by Sokolov et al. [Biophys. J. 2005, 89, 895.]) provides a basis for understanding recent experiments by Pelta et al. [Phys. Rev. Lett. 2007, 98, 228302.] on the kinetics of diffusion-limited gel degradation. These experiments find that the time t(c) taken by the enzyme thermolysin to completely hydrolyze a gel varies inversely as roughly the 3/2 power of the initial enzyme concentration [E]. In general, however, this time would be expected to vary either as [E](-1) or as [E](-2), depending on whether the Brownian diffusion of the enzyme to the site of cleavage took place along the network chains (1-d diffusion) or through the pore spaces (3-d diffusion). In our model, the unusual dependence of t(c) on [E] is explained in terms of a reaction-diffusion equation that is formulated in terms of fractional rather than ordinary time derivatives.

General model of a foam bed reactor

A general model of a foam bed reactor has been developed which rigorously accounts for the extent of gas absorption with chemical reaction occurring in both the storage and foam sections. Its applicability extends to a wide spectrum of reaction velocities. The possibilities of the predominance of the bulk-liquid reaction in the storage section or the absorption with reaction in the foam section can be handled as merely special cases of the general analysis. The importance of foam for carrying out a particular gas-liquid reaction is characterised by a criterion in terms of the fractional rate of reaction in the foam section. Trends of variations in the concentrations of dissolved free A, solute B, and gas-phase A with time of operation of the reactor are presented. The nature of the variation in the fractional rate of reaction in the foam section with time, at different reaction velocities, and the effect of the liquid flow rate (across the storage section) on the transience are also illustrated. Finally, the predictions of the general model have been validated using the available experimental data on the oxidation of sodium sulphide in a foam bed reactor. The agreement between the experimental and the present theoretical information is fairly good, apart from being more insightful than all the previous models of this reactor.

Sensitivity of wetland methane emissions to model assumptions: application and model testing against site observations

Abstract. Methane emissions from natural wetlands and rice paddies constitute a large proportion of atmospheric methane, but the magnitude and year-to-year variation of these methane sources is still unpredictable. Here we describe and evaluate the integration of a methane biogeochemical model (CLM4Me; Riley et al., 2011) into the Community Land Model 4.0 (CLM4CN) in order to better explain spatial and temporal variations in methane emissions. We test new functions for soil pH and redox potential that impact microbial methane production in soils. We also constrain aerenchyma in plants in always-inundated areas in order to better represent wetland vegetation. Satellite inundated fraction is explicitly prescribed in the model because there are large differences between simulated fractional inundation and satellite observations. A rice paddy module is also incorporated into the model, where the fraction of land used for rice production is explicitly prescribed. The model is evaluated at the site level with vegetation cover and water table prescribed from measurements. Explicit site level evaluations of simulated methane emissions are quite different than evaluating the grid cell averaged emissions against available measurements. Using a baseline set of parameter values, our model-estimated average global wetland emissions for the period 1993–2004 were 256 Tg CH4 yr−1, and rice paddy emissions in the year 2000 were 42 Tg CH4 yr−1. Tropical wetlands contributed 201 Tg CH4 yr−1, or 78 % of the global wetland flux. Northern latitude (>50 N) systems contributed 12 Tg CH4 yr−1. We expect this latter number may be an underestimate due to the low high-latitude inundated area captured by satellites and unrealistically low high-latitude productivity and soil carbon predicted by CLM4. Sensitivity analysis showed a large range (150–346 Tg CH4 yr−1) in predicted global methane emissions. The large range was sensitive to: (1) the amount of methane transported through aerenchyma, (2) soil pH (± 100 Tg CH4 yr−1), and (3) redox inhibition (± 45 Tg CH4 yr−1).

Microscopic study of the 2/5 fractional quantum Hall edge

This paper reports on our study of the edge of the 2/5 fractional quantum Hall state, which is more complicated than the edge of the 1/3 state because of the presence of edge sectors corresponding to different partitions of composite fermions in the lowest two Lambda levels. The addition of an electron at the edge is a nonperturbative process and it is not a priori obvious in what manner the added electron distributes itself over these sectors. We show, from a microscopic calculation, that when an electron is added at the edge of the ground state in the [N(1), N(2)] sector, where N(1) and N(2) are the numbers of composite fermions in the lowest two Lambda levels, the resulting state lies in either [N(1) + 1, N(2)] or [N(1), N(2) + 1] sectors; adding an electron at the edge is thus equivalent to adding a composite fermion at the edge. The coupling to other sectors of the form [N(1) + 1 + k, N(2) - k], k integer, is negligible in the asymptotically low-energy limit. This study also allows a detailed comparison with the two-boson model of the 2/5 edge. We compute the spectral weights and find that while the individual spectral weights are complicated and nonuniversal, their sum is consistent with an effective two-boson description of the 2/5 edge.

Experimental demonstration of a transparent graphene millimetre wave absorber with 28% fractional bandwidth at 140 GHz.

The development of transparent radio-frequency electronics has been limited, until recently, by the lack of suitable materials. Naturally thin and transparent graphene may lead to disruptive innovations in such applications. Here, we realize optically transparent broadband absorbers operating in the millimetre wave regime achieved by stacking graphene bearing quartz substrates on a ground plate. Broadband absorption is a result of mutually coupled Fabry-Perot resonators represented by each graphene-quartz substrate. An analytical model has been developed to predict the absorption performance and the angular dependence of the absorber. Using a repeated transfer-and-etch process, multilayer graphene was processed to control its surface resistivity. Millimetre wave reflectometer measurements of the stacked graphene-quartz absorbers demonstrated excellent broadband absorption of 90% with a 28% fractional bandwidth from 125-165 GHz. Our data suggests that the absorbers' operation can also be extended to microwave and low-terahertz bands with negligible loss in performance.

Cardiac beta ARK1 inhibition prolongs survival and augments beta blocker therapy in a mouse model of severe heart failure.

Chronic human heart failure is characterized by abnormalities in beta-adrenergic receptor (betaAR) signaling, including increased levels of betaAR kinase 1 (betaARK1), which seems critical to the pathogenesis of the disease. To determine whether inhibition of betaARK1 is sufficient to rescue a model of severe heart failure, we mated transgenic mice overexpressing a peptide inhibitor of betaARK1 (betaARKct) with transgenic mice overexpressing the sarcoplasmic reticulum Ca(2+)-binding protein, calsequestrin (CSQ). CSQ mice have a severe cardiomyopathy and markedly shortened survival (9 +/- 1 weeks). In contrast, CSQ/betaARKct mice exhibited a significant increase in mean survival age (15 +/- 1 weeks; P < 0.0001) and showed less cardiac dilation, and cardiac function was significantly improved (CSQ vs. CSQ/betaARKct, left ventricular end diastolic dimension 5.60 +/- 0.17 mm vs. 4.19 +/- 0.09 mm, P < 0.005; % fractional shortening, 15 +/- 2 vs. 36 +/- 2, P < 0.005). The enhancement of the survival rate in CSQ/betaARKct mice was substantially potentiated by chronic treatment with the betaAR antagonist metoprolol (CSQ/betaARKct nontreated vs. CSQ/betaARKct metoprolol treated, 15 +/- 1 weeks vs. 25 +/- 2 weeks, P < 0.0001). Thus, overexpression of the betaARKct resulted in a marked prolongation in survival and improved cardiac function in a mouse model of severe cardiomyopathy that can be potentiated with beta-blocker therapy. These data demonstrate a significant synergy between an established heart-failure treatment and the strategy of betaARK1 inhibition.

Analytical calculation of far infrared spectra of ice in terms of a molecular model

Wideband far infrared (FIR) spectra of complex permittivity e(p) of ice are calculated in terms of a simple analytical theory based on the method of dipolar autocorrelation functions. The molecular model represents a revision of the model recently presented for liquid water in Adv. Chem. Phys. 127 (2003) 65. A composite two-fractional model is proposed. The model is characterised by three phenomenological potential wells corresponding to the three FIR bands observed in ice. The first fraction comprises dipoles reorienting in a rather narrow and deep hat-like well; these dipoles generate the librational band centred at the frequency approximate to 880 cm(-1). The second fraction comprises elastically interacting particles; they generate two nearby bands placed around frequency 200 cm(-1). For description of one of these bands the harmonic oscillator (HO) model is used, in which translational oscillations of two charged molecules along the H-bond are considered. The other band is produced by the H-bond stretch, which governs hindered rotation of a rigid dipole. Such a motion and its dielectric response are described in terms of a new cut parabolic (CP) model applicable for any vibration amplitude. The composite hat-HO-CP model results in a smooth epsilon(nu) ice spectrum, which does not resemble the noise-like spectra of ice met in the known literature. The proposed theory satisfactorily agrees with the experimental ice spectrum measured at - 7 degrees C. The calculated longitudinal optic-transverse optic (LO-TO) splitting occurring at approximate to 250 cm(-1) qualitatively agrees with the measured data. (c) 2004 Elsevier B.V. All rights reserved.

A new look at the pulsating DB white dwarf GD 358: Line-of-sight velocity measurements and constraints on model atmospheres

We report on our findings of the bright, pulsating, helium atmosphere white dwarf GD 358, based on time-resolved optical spectrophotometry. We identify 5 real pulsation modes and at least 6 combination modes at frequencies consistent with those found in previous observations. The measured Doppler shifts from our spectra show variations with amplitudes of up to 5.5 km s-1 at the frequencies inferred from the flux variations. We conclude that these are variations in the line-of-sight velocities associated with the pulsational motion. We use the observed flux and velocity amplitudes and phases to test theoretical predictions within the convective driving framework, and compare these with similar observations of the hydrogen atmosphere white dwarf pulsators (DAVs). The wavelength dependence of the fractional pulsation amplitudes (chromatic amplitudes) allows us to conclude that all five real modes share the same spherical degree, most likely, l=1. This is consistent with previous identifications based solely on photometry. We find that a high signal-to-noise mean spectrum on its own is not enough to determine the atmospheric parameters and that there are small but significant discrepancies between the observations and model atmospheres. The source of these remains to be identified. While we infer Teff =24 kK and log g ~ 8.0 from the mean spectrum, the chromatic amplitudes, which are a measure of the derivative of the flux with respect to the temperature, unambiguously favour a higher effective temperature, 27 kK, which is more in line with independent determinations from ultra-violet spectra.

Artificial surfaces of intertwined square spirals:A CPW model

Intertwining planar spirals arranged in doubly periodic arrays enables a substantially subwavelength response of the unit cell smaller than 1/40 of wavelength with large fractional bandwidths. These properties are important for application at low frequencies, conformal curved surfaces, or with compact radiators. It is shown that interleaving counter-wound spiral arms extended into adjacent unit cells dramatically increase the array equivalent capacitance while reducing the inductance. A coplanar waveguide (CPW) model has been developed to analytically estimate the equivalent capacitance and inductance of intertwined spiral array elements in terms of their geometrical parameters. The proposed CPW model is shown to provide an accurate prediction of the fundamental resonance frequency and can be instrumental in the design of the arrays for a specified frequency response. © 2012 IEEE.

Analytical model of interwoven spiral arrays

Interweaving planar spiral conductors in doubly periodic arrays enable substantially sub-wavelength resonant response along with broadening fractional bandwidth. A self-contained analytical model is proposed to accurately predict the characteristics of the intertwined quadrifilar spiral array near the fundamental resonance. The model, based upon a multiconductor transmission line (MTL) approach, provides physical insight into the unique properties of the distributed interactions between the interleaved counter-wound spiral arms extended beyond a single unit cell and elucidates the mechanisms underlying the array performance at normal and oblique incidence of TE and TM polarised waves. The developed MTL model is instrumental in the design of the artificial surfaces with the specified response.

Forecasting overseas visitors into the United Kingdom using continuous time and autoregressive fractional integrated moving average models with discrete data

This paper applies Gaussian estimation methods to continuous time models for modelling overseas visitors into the UK. The use of continuous time modelling is widely used in economics and finance but not in tourism forecasting. Using monthly data for 1986–2010, various continuous time models are estimated and compared to autoregressive integrated moving average (ARIMA) and autoregressive fractionally integrated moving average (ARFIMA) models. Dynamic forecasts are obtained over different periods. The empirical results show that the ARIMA model performs very well, but that the constant elasticity of variance (CEV) continuous time model has the lowest root mean squared error (RMSE) over a short period.

Fractional order modelling of zero length column desorption response for adsorbents with variable particle sizes

This manuscript analyses the data generated by a Zero Length Column (ZLC) diffusion experimental set-up, for 1,3 Di-isopropyl benzene in a 100% alumina matrix with variable particle size. The time evolution of the phenomena resembles those of fractional order systems, namely those with a fast initial transient followed by long and slow tails. The experimental measurements are best fitted with the Harris model revealing a power law behavior.