Fractal and multifractal properties of a family of fractal networks

In this work, we study the fractal and multifractal properties of a family of fractal networks introduced by Gallos et al (2007 Proc. Nat. Acad. Sci. USA 104 7746). In this fractal network model, there is a parameter e which is between 0 and 1, and allows for tuning the level of fractality in the network. Here we examine the multifractal behavior of these networks, the dependence relationship of the fractal dimension and the multifractal parameters on parameter e. First, we find that the empirical fractal dimensions of these networks obtained by our program coincide with the theoretical formula given by Song et al (2006 Nature Phys. 2 275). Then from the shape of the τ(q) and D(q) curves, we find the existence of multifractality in these networks. Last, we find that there exists a linear relationship between the average information dimension 〈D(1)〉 and the parameter e.

Combining transcriptome assemblies from multiple de novo assemblers in the allo-tetraploid plant Nicotiana benthamiana

Background Nicotiana benthamiana is an allo-tetraploid plant, which can be challenging for de novo transcriptome assemblies due to homeologous and duplicated gene copies. Transcripts generated from such genes can be distinct yet highly similar in sequence, with markedly differing expression levels. This can lead to unassembled, partially assembled or mis-assembled contigs. Due to the different properties of de novo assemblers, no one assembler with any one given parameter space can re-assemble all possible transcripts from a transcriptome. Results In an effort to maximise the diversity and completeness of de novo assembled transcripts, we utilised four de novo transcriptome assemblers, TransAbyss, Trinity, SOAPdenovo-Trans, and Oases, using a range of k-mer sizes and different input RNA-seq read counts. We complemented the parameter space biologically by using RNA from 10 plant tissues. We then combined the output of all assemblies into a large super-set of sequences. Using a method from the EvidentialGene pipeline, the combined assembly was reduced from 9.9 million de novo assembled transcripts to about 235,000 of which about 50,000 were classified as primary. Metrics such as average bit-scores, feature response curves and the ability to distinguish paralogous or homeologous transcripts, indicated that the EvidentialGene processed assembly was of high quality. Of 35 RNA silencing gene transcripts, 34 were identified as assembled to full length, whereas in a previous assembly using only one assembler, 9 of these were partially assembled. Conclusions To achieve a high quality transcriptome, it is advantageous to implement and combine the output from as many different de novo assemblers as possible. We have in essence taking the ‘best’ output from each assembler while minimising sequence redundancy. We have also shown that simultaneous assessment of a variety of metrics, not just focused on contig length, is necessary to gauge the quality of assemblies.

Fractal and complex network analyses of protein molecular dynamics

Based on protein molecular dynamics, we investigate the fractal properties of energy, pressure and volume time series using the multifractal detrended fluctuation analysis (MF-DFA) and the topological and fractal properties of their converted horizontal visibility graphs (HVGs). The energy parameters of protein dynamics we considered are bonded potential, angle potential, dihedral potential, improper potential, kinetic energy, Van der Waals potential, electrostatic potential, total energy and potential energy. The shape of the h(q)h(q) curves from MF-DFA indicates that these time series are multifractal. The numerical values of the exponent h(2)h(2) of MF-DFA show that the series of total energy and potential energy are non-stationary and anti-persistent; the other time series are stationary and persistent apart from series of pressure (with H≈0.5H≈0.5 indicating the absence of long-range correlation). The degree distributions of their converted HVGs show that these networks are exponential. The results of fractal analysis show that fractality exists in these converted HVGs. For each energy, pressure or volume parameter, it is found that the values of h(2)h(2) of MF-DFA on the time series, exponent λλ of the exponential degree distribution and fractal dimension dBdB of their converted HVGs do not change much for different proteins (indicating some universality). We also found that after taking average over all proteins, there is a linear relationship between 〈h(2)〉〈h(2)〉 (from MF-DFA on time series) and 〈dB〉〈dB〉 of the converted HVGs for different energy, pressure and volume.

Underlying scaling relationships between solar activity and geomagnetic activity revealed by multifractal analyses

his paper identifies some scaling relationships between solar activity and geomagnetic activity. We examine the scaling properties of hourly data for two geomagnetic indices (ap and AE), two solar indices (solar X-rays Xl and solar flux F10.7), and two inner heliospheric indices (ion density Ni and flow speed Vs) over the period 1995–2001 by the universal multifractal approach and the traditional multifractal analysis. We found that the universal multifractal model (UMM) provides a good fit to the empirical K(q) and τ(q) curves of these time series. The estimated values of the Lévy index α in the UMM indicate that multifractality exists in the time series for ap, AE, Xl, and Ni, while those for F10.7 and Vs are monofractal. The estimated values of the nonconservation parameter H of this model confirm that these time series are conservative which indicate that the mean value of the process is constant for varying resolution. Additionally, the multifractal K(q) and τ(q) curves, and the estimated values of the sparseness parameter C1 of the UMM indicate that there are three pairs of indices displaying similar scaling properties, namely ap and Xl, AE and Ni, and F10.7 and Vs. The similarity in the scaling properties of pairs (ap,Xl) and (AE,Ni) suggests that ap and Xl, AE and Ni are better correlated—in terms of scaling—than previous thought, respectively. But our results still cannot be used to advance forecasting of ap and AE by Xl and Ni, respectively, due to some reasons

Simulating the dissipation of two herbicides using micro paddy lysimeters

A set of packed micro paddy lysimeters, placed in a greenhouse, was used to simulate the dissipation of two herbicides, simetryn and thiobencarb, in a controlled environment. Data from a field monitoring study in 2003, including the soil condition and water balances, were used in the simulation. The herbicides were applied and monitored over a period of 21 d. The water balances under two water management scenarios, intermittent irrigation management (AI) and continuous irrigation management (CI), were simulated. In the AI scenario, the pattern of herbicide dissipation in the surface water of the field were simulated, following the first-order kinetics. In the CI scenario, similarity was observed in most lysimeter and field concentrations, but there were differences in some data points. Dissipation curves of both herbicides in the surface water of the two simulated scenarios were not significantly different (P > 0.05) from the field data except for intercept of the thiobencarb curve in the CI scenario. The distribution of simetryn and thiobencarb in the soil profile after simulation were also similar to the field data. The highest concentrations of both herbicides were found on the topsoil layer at 0-2.5 cm depth. Only a small amount of herbicides moved down to the deeper soil layers. Micro paddy lysimeters are thus a good alternative for the dissipation study of pesticides in the paddy environment.

Finite element analysis investigation of response of pumpkin tissue under uniaxial loading

Finite element analysis (FEA) models of uniaxial loading of pumpkin peel and flesh tissues were developed and validated using experimental results. The tensile model was developed for both linear elastic and plastic material models, the compression model was develop d only with the plastic material model. The outcomes of force versus time curves obtained from FEA models followed similar pattern to the experimental curves however the curve resulted with linear elastic material properties had a higher difference with the experimental curves. The values of predicted forces were determined and compared with the experimental curve. An error indicator was introduced and computed for each case and compared. Additionally Root Mean Square Error (RMSE) values were also calculated for each model and compared. The results of modelling were used to develop material model for peel and flesh tissues in FEA modelling of mechanical peeling of tough skinned vegetables.

Yielding of sensitive clays: micromechanical considerations

Test results reported on several natural sensitive soils show significant anisotropy of the yield curves, which are generally oriented along the coefficient of earth pressure at rest (K-0) axis. An attempt is made in this paper to explain the anisotropy in yielding from microstructural considerations. An elliptic pore, with particle domains aligned along the periphery of the pore, and with the major axis of the pore being oriented along the direction of the in situ major principal stress, is chosen as the unit of microstructure. An analysis of forces at the interdomain contacts around the ellipse is carried out with reference to experimentally determined yield stress conditions of one soil, and a yield criteria is defined. The analysis, with the proposed yield criteria, enables one to define the complete yield curve for any other soil from the results of only two tests (one constant eta compression test with eta close to eta(K?0), where eta is the stress ratio (= q/p) and eta(K?0) is the stress ratio corresponding to anisotropic K-0 compression, and another undrained shear test). Predicted yield curves are compared with experimental yield curves of several soils reported in the literature.

Phase conversions in calcium manganites with changing Ca/Mn ratios and their influence on the electrical transport properties

The phase-interconversions between the spinel-, brownmillerite-, defect rocksalt and perovskite-type structures have been investigated by way of (i) introducing deficiency in A-sites in CaxMn2-xO3 (0.05 <= x <= 1) i.e., by varying Ca/Mn ratio from 0.025 to 1 and (ii) nonstoichiometric CaMnO3-delta (CMO) with 0.02 <= delta <= 1. The temperature dependence of resistivity (rho-T) have been investigated on nonstoichiometric CaMnO3-delta (undoped) as well as the CMO substituted with donor impurities such as La3+, Y3+, Bi3+ or acceptor such as Na1+ ion at the Ca-site. The rho-T characteristics of nonstoichiometric CaMnO3-delta is strongly influenced by oxygen deficiency, which controls the concentration of Mn3+ ions and, in turn, affects the resistivity, rho. The results indicated that the substitution of aliovalent impurities at Ca-site in CaMnO3 has similar effects as of CaMnO3-delta ( undoped) annealed in atmospheres of varying partial pressures whereby electron or hole concentration can be altered, yet the doped samples can be processed in air or atmospheres of higher P-O2. The charge transport mechanisms of nonstoichiometric CaMnO3-delta as against the donor or acceptor doped CaMnO3 (sintered in air, P-O2 similar to 0.2 atm) have been predicted. The rho (T) curves of both donor doped CaMnO3 as well as non-stoichiometric CaMnO3-delta, is predictable by the small polaron hopping (SPH) model, which changes to the variable range hopping (VRH) at low temperatures whereas the acceptor doped CaMnO3 exhibited an activated semiconducting hopping ( ASH) throughout the measured range of temperature (10-500 K).

Species-area and species-individual relationships for tropical trees: A comparison of three 50-ha plots

1 Species-accumulation curves for woody plants were calculated in three tropical forests, based on fully mapped 50-ha plots in wet, old-growth forest in Peninsular Malaysia, in moist, old-growth forest in central Panama, and in dry, previously logged forest in southern India. A total of 610 000 stems were identified to species and mapped to < Im accuracy. Mean species number and stem number were calculated in quadrats as small as 5 m x 5 m to as large as 1000 m x 500 m, for a variety of stem sizes above 10 mm in diameter. Species-area curves were generated by plotting species number as a function of quadrat size; species-individual curves were generated from the same data, but using stem number as the independent variable rather than area. 2 Species-area curves had different forms for stems of different diameters, but species-individual curves were nearly independent of diameter class. With < 10(4) stems, species-individual curves were concave downward on log-log plots, with curves from different forests diverging, but beyond about 104 stems, the log-log curves became nearly linear, with all three sites having a similar slope. This indicates an asymptotic difference in richness between forests: the Malaysian site had 2.7 times as many species as Panama, which in turn was 3.3 times as rich as India. 3 Other details of the species-accumulation relationship were remarkably similar between the three sites. Rectangular quadrats had 5-27% more species than square quadrats of the same area, with longer and narrower quadrats increasingly diverse. Random samples of stems drawn from the entire 50 ha had 10-30% more species than square quadrats with the same number of stems. At both Pasoh and BCI, but not Mudumalai. species richness was slightly higher among intermediate-sized stems (50-100mm in diameter) than in either smaller or larger sizes, These patterns reflect aggregated distributions of individual species, plus weak density-dependent forces that tend to smooth the species abundance distribution and 'loosen' aggregations as stems grow. 4 The results provide support for the view that within each tree community, many species have their abundance and distribution guided more by random drift than deterministic interactions. The drift model predicts that the species-accumulation curve will have a declining slope on a log-log plot, reaching a slope of O.1 in about 50 ha. No other model of community structure can make such a precise prediction. 5 The results demonstrate that diversity studies based on different stem diameters can be compared by sampling identical numbers of stems. Moreover, they indicate that stem counts < 1000 in tropical forests will underestimate the percentage difference in species richness between two diverse sites. Fortunately, standard diversity indices (Fisher's sc, Shannon-Wiener) captured diversity differences in small stem samples more effectively than raw species richness, but both were sample size dependent. Two nonparametric richness estimators (Chao. jackknife) performed poorly, greatly underestimating true species richness.

Fuzzy optimization model for water quality management of a river system - Closure

A fuzzy waste-load allocation model, FWLAM, is developed for water quality management of a river system using fuzzy multiple-objective optimization. An important feature of this model is its capability to incorporate the aspirations and conflicting objectives of the pollution control agency and dischargers. The vagueness associated with specifying the water quality criteria and fraction removal levels is modeled in a fuzzy framework. The goals related to the pollution control agency and dischargers are expressed as fuzzy sets. The membership functions of these fuzzy sets are considered to represent the variation of satisfaction levels of the pollution control agency and dischargers in attaining their respective goals. Two formulations—namely, the MAX-MIN and MAX-BIAS formulations—are proposed for FWLAM. The MAX-MIN formulation maximizes the minimum satisfaction level in the system. The MAX-BIAS formulation maximizes a bias measure, giving a solution that favors the dischargers. Maximization of the bias measure attempts to keep the satisfaction levels of the dischargers away from the minimum satisfaction level and that of the pollution control agency close to the minimum satisfaction level. Most of the conventional water quality management models use waste treatment cost curves that are uncertain and nonlinear. Unlike such models, FWLAM avoids the use of cost curves. Further, the model provides the flexibility for the pollution control agency and dischargers to specify their aspirations independently.

Fatigue Laws Via Functional Equations

In routine industrial design, fatigue life estimation is largely based on S-N curves and ad hoc cycle counting algorithms used with Miner's rule for predicting life under complex loading. However, there are well known deficiencies of the conventional approach. Of the many cumulative damage rules that have been proposed, Manson's Double Linear Damage Rule (DLDR) has been the most successful. Here we follow up, through comparisons with experimental data from many sources, on a new approach to empirical fatigue life estimation (A Constructive Empirical Theory for Metal Fatigue Under Block Cyclic Loading', Proceedings of the Royal Society A, in press). The basic modeling approach is first described: it depends on enforcing mathematical consistency between predictions of simple empirical models that include indeterminate functional forms, and published fatigue data from handbooks. This consistency is enforced through setting up and (with luck) solving a functional equation with three independent variables and six unknown functions. The model, after eliminating or identifying various parameters, retains three fitted parameters; for the experimental data available, one of these may be set to zero. On comparison against data from several different sources, with two fitted parameters, we find that our model works about as well as the DLDR and much better than Miner's rule. We finally discuss some ways in which the model might be used, beyond the scope of the DLDR.

Numerical prediction of load-displacement behaviors of adhesively bonded joints at different extension rates and temperatures

The present work focuses on simulation of nonlinear mechanical behaviors of adhesively bonded DLS (double lap shear) joints for variable extension rates and temperatures using the implicit ABAQUS solver. Load-displacement curves of DLS joints at nine combinations of extension rates and environmental temperatures are initially obtained by conducting tensile tests in a UTM. The joint specimens are made from dual phase (DP) steel coupons bonded with a rubber-toughened adhesive. It is shown that the shell-solid model of a DLS joint, in which substrates are modeled with shell elements and adhesive with solid elements, can effectively predict the mechanical behavior of the joint. Exponent Drucker-Prager or Von Mises yield criterion together with nonlinear isotropic hardening is used for the simulation of DLS joint tests. It has been found that at a low temperature (-20 degrees C), both Von Mises and exponent Drucker-Prager criteria give close prediction of experimental load-extension curves. However. at a high temperature (82 degrees C), Von Mises condition tends to yield a perceptibly softer joint behavior, while the corresponding response obtained using exponent Drucker-Prager criterion is much closer to the experimental load-displacement curve.

Consolidation Behavior of Soils

Based on Terzaghi's consolidation theory, percent of consolidation, U, versus the time factor, T, relationship for constant/linear excess pore water pressure distribution, it is possible to generate theoretical log10(H2/t) versus U curves where H is the length of the drainage path of a consolidating layer, and t is the time for different known values of the coefficient of consolidation, cν. A method has been developed wherein both the theoretical and experimental behavior of soils during consolidation can be simultaneously compared and studied on the same plot. The experimental log10(H2/t) versus U curves have been compared with the theoretical curves. The deviations of the experimental behavior from the theory are explained in terms of initial compression and secondary compression. Analysis of results indicates that the secondary compression essentially starts from about 60% consolidation. A simple procedure is presented for calculating the value of cv from the δ-t data using log10(H2/t) versus U plot.

Melanopsin mediated post-illumination pupil response in early age-related macular degeneration

Purpose To determine whether melanopsin expressing intrinsically photosensitive Retinal Ganglion Cell (ipRGC) inputs to the pupil light reflex (PLR) are affected in early age-related macular degeneration (AMD). Methods The PLR was measured in 40 participants (20 early AMD and 20 age-matched controls) using a custom-built Maxwellian-view pupillometer. Sinusoidal stimuli (0.5 Hz, 11.9 s duration, 35.6° diameter) were presented to the study eye and the consensual pupil response was measured for stimuli with high melanopsin excitation (464nm; blue) and with low melanopsin excitation (638 nm; red) that biased activation to the outer retina. Two melanopsin PLR metrics were quantified: the Phase Amplitude Percentage (PAP) during the sinusoidal stimulus presentation and the Post-Illumination Pupil Response (PIPR). The PLR during stimulus presentation was analyzed using latency to constriction, transient pupil response and maximum pupil constriction metrics. Diagnostic accuracy was evaluated using receiver operating characteristic (ROC) curves. Results The blue PIPR was significantly less sustained in the early AMD group (p<0.001). The red PIPR was not significantly different between groups (p>0.05). The PAP and blue stimulus constriction amplitude were significantly lower in the early AMD group (p < 0.05). There was no significant difference between groups in the latency or transient amplitude for both stimuli (p>0.05). ROC analysis showed excellent diagnostic accuracy for the blue PIPR metrics (AUC>0.9). Conclusions This is the initial report that the melanopsin controlled PIPR is dysfunctional in early AMD. The non-invasive, objective measurement of the ipRGC controlled PIPR has excellent diagnostic accuracy for early AMD.

Ion beam processing of oriented CuO films deposited on (100) YSZ by laser ablation

A systematic study of Ar ion implantation in cupric oxide films has been reported. Oriented CuO films were deposited by pulsed excimer laser ablation technique on (1 0 0) YSZ substrates. X-ray diffraction (XRD) spectra showed the highly oriented nature of the deposited CuO films. The films were subjected to ion bombardment for studies of damage formation, Implantations were carried out using 100 keV Arf over a dose range between 5 x 10(12) and 5 x 10(15) ions/cm(2). The as-deposited and ion beam processed samples were characterized by XRD technique and resistance versus temperature (R-T) measurements. The activation energies for electrical conduction were found from In [R] versus 1/T curves. Defects play an important role in the conduction mechanism in the implanted samples. The conductivity of the film increases, and the corresponding activation energy decreases with respect to the dose value.