Wavelet correlation between subjects: A time-scale data driven analysis for brain mapping using fMRI

Functional magnetic resonance imaging (fMRI) based on BOLD signal has been used to indirectly measure the local neural activity induced by cognitive tasks or stimulation. Most fMRI data analysis is carried out using the general linear model (GLM), a statistical approach which predicts the changes in the observed BOLD response based on an expected hemodynamic response function (HRF). In cases when the task is cognitively complex or in cases of diseases, variations in shape and/or delay may reduce the reliability of results. A novel exploratory method using fMRI data, which attempts to discriminate between neurophysiological signals induced by the stimulation protocol from artifacts or other confounding factors, is introduced in this paper. This new method is based on the fusion between correlation analysis and the discrete wavelet transform, to identify similarities in the time course of the BOLD signal in a group of volunteers. We illustrate the usefulness of this approach by analyzing fMRI data from normal subjects presented with standardized human face pictures expressing different degrees of sadness. The results show that the proposed wavelet correlation analysis has greater statistical power than conventional GLM or time domain intersubject correlation analysis. (C) 2010 Elsevier B.V. All rights reserved.

Optimal quantum trajectories for discrete measurements

Using the method of quantum trajectories we show that a known pure state can be optimally monitored through time when subject to a sequence of discrete measurements. By modifying the way that we extract information from the measurement apparatus we can minimize the average algorithmic information of the measurement record, without changing the unconditional evolution of the measured system. We define an optimal measurement scheme as one which has the lowest average algorithmic information allowed. We also show how it is possible to extract information about system operator averages from the measurement records and their probabilities. The optimal measurement scheme, in the limit of weak coupling, determines the statistics of the variance of the measured variable directly. We discuss the relevance of such measurements for recent experiments in quantum optics.

Dynamic processes in the contingent valuation of an endangered mammal species

Reports experimental results involving 204 members of the public who were asked their willingness to pay for the conservation of the mahogany glider Petaurus gracilis on three occasions: prior to information being provided to them about the glider and other wildlife species; after such information was provided, and after participants had an opportunity to see live specimens of this endangered species. Variations in the mean willingness to pay are analysed. Concerns arise about whether information provision and experience reveal ‘true’ contingent valuations of public goods and about the choice of the relevant contingent valuation measure.

A dynamic model with on-line identification for rotary sugar drying processes

A dynamic modelling methodology, which combines on-line variable estimation and parameter identification with physical laws to form an adaptive model for rotary sugar drying processes, is developed in this paper. In contrast to the conventional rate-based models using empirical transfer coefficients, the heat and mass transfer rates are estimated by using on-line measurements in the new model. Furthermore, a set of improved sectional solid transport equations with localized parameters is developed in this work to reidentified on-line using measurement data, the model is able to closely track the dynamic behaviour of rotary drying processes within a broad range of operational conditions. This adaptive model is validated against experimental data obtained from a pilot-scale rotary sugar dryer. The proposed modelling methodology can be easily incorporated into nonlinear model based control schemes to form a unified modelling and control framework.place the global correlation for the computation of solid retention time. Since a number of key model variables and parameters are identified on-line using measurement data, the model is able to closely track the dynamic behaviour of rotary drying processes within a broad range of operational conditions. This adaptive model is validated against experimental data obtained from a pilot-scale rotary sugar dryer. The proposed modelling methodology can be easily incorporated into nonlinear model based control schemes to form a unified modelling and control framework.

Effect of probe stimulus intensity on the dissociation between autonomic orienting and secondary probe reaction time

Information processing accounts propose that autonomic orienting reflects the amount of resources allocated to process a stimulus. However, secondary task reaction time (RT), a supposed measure of processing resources, has shown a dissociation from autonomic orienting. The present study tested the hypothesis that secondary task RT reflects a serial processing mechanism. Participants (N = 24) were presented with circle and ellipse shapes and asked to count the number of longer-than-usual presentations of one shape (task-relevant) and to ignore presentations of a second shape (task-irrelevant). Concurrent with the counting task, participants performed a secondary RT task to an auditory probe presented at either a high or low intensity and at two different probe positions following shape onset (50 and 300 ms). Electrodermal orienting was larger during task-relevant shapes than during task-irrelevant shapes, but secondary task RT to the high-intensity probe was slower during the latter. In addition, an underadditive interaction between probe stimulus intensity and probe position was found in secondary RT. The findings are consistent with a serial processing model of secondary RT and suggest that the notion of processing stages should be incorporated into current information-processing models of autonomic orienting.

Step size control for efficient discrete element simulation

The step size determines the accuracy of a discrete element simulation. The position and velocity updating calculation uses a pre-calculated table and hence the control of step size can not use the integration formulas for step size control. A step size control scheme for use with the table driven velocity and position calculation uses the difference between the calculation result from one big step and that from two small steps. This variable time step size method chooses the suitable time step size for each particle at each step automatically according to the conditions. Simulation using fixed time step method is compared with that of using variable time step method. The difference in computation time for the same accuracy using a variable step size (compared to the fixed step) depends on the particular problem. For a simple test case the times are roughly similar. However, the variable step size gives the required accuracy on the first run. A fixed step size may require several runs to check the simulation accuracy or a conservative step size that results in longer run times. (C) 2001 Elsevier Science Ltd. All rights reserved.

Theory of pseudomodes in quantum optical processes

This paper deals with non-Markovian behavior in atomic systems coupled to a structured reservoir of quantum electromagnetic field modes, with particular relevance to atoms interacting with the field in high-Q cavities or photonic band-gap materials. In cases such as the former, we show that the pseudomode theory for single-quantum reservoir excitations can be obtained by applying the Fano diagonalization method to a system in which the atomic transitions are coupled to a discrete set of (cavity) quasimodes, which in turn are coupled to a continuum set of (external) quasimodes with slowly varying coupling constants and continuum mode density. Each pseudomode can be identified with a discrete quasimode, which gives structure to the actual reservoir of true modes via the expressions for the equivalent atom-true mode coupling constants. The quasimode theory enables cases of multiple excitation of the reservoir to now be treated via Markovian master equations for the atom-discrete quasimode system. Applications of the theory to one, two, and many discrete quasimodes are made. For a simple photonic band-gap model, where the reservoir structure is associated with the true mode density rather than the coupling constants, the single quantum excitation case appears to be equivalent to a case with two discrete quasimodes.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Analytical validation for a series of marker compounds used to assess renal drug elimination processes

Renal drug elimination is determined by glomerular filtration, tubular secretion, and tubular reabsorption. Changes in the integrity of these processes influence renal drug clearance, and these changes may not be detected by conventional measures of renal function such as creatinine clearance. The aim of the current study was to examine the analytic issues needed to develop a cocktail of marker drugs (fluconazole, rac-pindolol, para-aminohippuric acid, sinistrin) to measure simultaneously the mechanisms contributing to renal clearance. High-performance liquid chromatographic methods of analysis for fluconazole, pindolol, para-aminohippuric acid, and creatinine and an enzymatic assay for sinistrin were developed or modified and then validated to allow determination of each of the compounds in both plasma and urine in the presence of all other marker drugs. A pilot clinical study in one volunteer was conducted to ensure that the assays were suitable for quantitating all the marker drugs to the sensitivity and specificity needed to allow accurate determination of individual renal clearances. The performance of all assays (plasma and urine) complied with published validation criteria. All standard curves displayed linearity over the concentration ranges required, with coefficients of correlation greater than 0.99. The precision of the interday and intraday variabilities of quality controls for each marker in plasma and urine were all less than 11.9% for each marker. Recoveries of markers (and internal standards) in plasma and urine were all at least 90%. All markers investigated were shown to be stable when plasma or urine was frozen and thawed. For all the assays developed, there were no interferences from other markers or endogenous substances. In a pilot clinical study, concentrations of all markers could be accurately and reproducibly determined for a sufficient duration of time after administration to calculate accurate renal clearance for each marker. This article presents details of the analytic techniques developed for measuring concentrations of marker drugs for different renal elimination processes administered as a single dose to define the processes contributing to renal drug elimination.

Non-Markovian quantum optical processes and pseudomodes

Non-Markovian behaviour in atomic systems coupled to a structured reservoir of quantum EM field modes, such as in high Q cavities, is treated using a quasimode description, and the pseudo mode theory for single quantum reservoir excitations is obtained via Fano diagonalisation. The atomic transitions are coupled to a discrete set of (cavity) quasimodes, which are also coupled to a continuum set of (external) quasimodes with slowly varying coupling constants. Each pseudomode corresponds to a cavity quasimode, and the original reservoir structure is obtained in expressions for the equivalent atom-true mode coupling constants. Cases of multiple excitation of the reservoir are now treatable via Markovian master equations for the atom-discrete quasimode system.

Path integrals for continuous-time Markov chains

This note presents a method of evaluating the distribution of a path integral for Markov chains on a countable state space.

Integrals for continuous-time Markov chains

This paper presents a method of evaluating the expected value of a path integral for a general Markov chain on a countable state space. We illustrate the method with reference to several models, including birth-death processes and the birth, death and catastrophe process. (C) 2002 Elsevier Science Inc. All rights reserved.

Lanczos subspace time-independent wave packet calculations of S(D-1)+H-2 reactive scattering

In this paper. we present the results of quantum dynamical simulations of the S (D-1) + H-2 insertion reaction on a newly developed potential energy surface (J. Chem. Phys. 2001, 114, 320). State-to-state reaction probabilities. product state distributions, and initial-state resolved cumulative reaction probabilities from a given incoming reactant channel are obtained from a time-independent wave packet analysis, performed within a single Lanczos subspace. Integral reaction cross sections are then estimated by J-shifting method and compared with the results from molecular beam experiment and QCT calculations.

Kinetics of baculovirus replication and release using real-time quantitative polymerase chain reaction

The study of viral-based processes is hampered by (a) their complex, transient nature, (b) the instability of products, and (c) the lack of accurate diagnostic assays. Here, we describe the use of real-time quantitative polymerase chain reaction to characterize baculoviral infection. Baculovirus DNA content doubles every 1.7 h from 6 h post-infection until replication is halted at the onset of budding. No dynamic equilibrium exists between replication and release, and the kinetics are independent of the cell density at the time of infection. No more than 16% of the intracellular virus copies bud from the cell. (C) 2002 John Wiley & Sons, Inc. Biotechnol Bioeng 77: 476-480, 2002; DOI 10.1002/bit.10126.

Transformation processes, pathways, and possible sources of distinctive polychlorinated dibenzo-p-dioxin signatures in sink environments

In recent years, studies on environmental samples with unusual dibenzo-p-dioxin (PCDD) congener profiles were reported from a range of countries. These profiles, characterized by a dominance of octachlorinated dibenzodioxin (OCDD) and relatively low in dibenzofuran (PCDF) concentrations, could not be attributed to known sources or formation processes. In the present study, the processes that result in these unusual profiles were assessed using the concentrations and isomer signatures of PCDDs from dated estuarine sediment cores in Queensland, Australia. Increases in relative concentrations of lower chlorinated PODS and a relative decrease of OCDD were correlated with time of sediment deposition. Preferred lateral, anaerobic dechlorination of OCDD represents a likely pathway for these changes. In Queensland sediments, these transformations result in a distinct dominance of isomers fully chlorinated in the 1,4,6,9-positions (1,4-patterns), and similar 1,4-patterns were observed in sediments from elsewhere. Consequently, these environmental samples may not reflect the signatures of the original source, and a reevaluation of source inputs was undertaken. Natural formation of PCDDs, which has previously been suggested, is discussed; however, based on the present results and literature comparisons, we propose an alternative scenario. This scenario hypothesizes that an anthropogenic PCDD precursor input (e.g. pentachlorophenol) results in the contamination. These results and hypothesis imply further investigations are warrented into possible anthropogenic sources in areas where natural PCDD formation has been suggested.