Structural and optical properties of GaAs(100) with a thin surface layer doped with chromium

DOE Office of Scientific and Technical Information (OSTI.GOV)

Seredin, P. V., E-mail: paul@phys.vsu.ru; Fedyukin, A. V.; Arsentyev, I. N.

The aim of this study is to explore the structural and optical properties of single-crystal GaAs(100) doped with Cr atoms by burning them into the substrate at high temperatures. The diffusion of chromium into single-crystal GaAs(100) substrates brings about the formation of a thin (~20–40 μm) GaAs:Cr transition layer. In this case, chromium atoms are incorporated into the gallium-arsenide crystal lattice and occupy the regular atomic sites of the metal sublattice. As the chromium diffusion time is increased, such behavior of the dopant impurity yields changes in the energy structure of GaAs, a decrease in the absorption at free chargemore » carriers, and a lowering of the surface recombination rate. As a result, the photoluminescence signal from the sample is significantly enhanced.« less

Assessment of diffuser pressure loss on WWTPs in Baden-Württemberg.

PubMed

Krampe, J

2011-01-01

Aeration of activated sludge is a critical treatment step for the operation of activated sludge plants. To achieve a cost effective treatment process, assessing and benchmarking of aeration system performance are important measures. A simple means of gauging the relative condition of a fine bubble diffused aeration system is to evaluate the pressure loss of the diffusers as oxygen transfer tests are rarely applied during the lifetime of an aeration system. This paper shows an assessment of fine bubble diffuser systems in Baden-Württemberg, Germany, based on the results of a questionnaire sent to 941 WWTPs. Apart from the results with regards to the diffuser pressure loss, this paper also presents information on the current state of diffuser technology such as types and materials as well as the diffuser cleaning methods used in Baden-Württemberg. The majority of the WWTPs were equipped with tube diffusers (71%) with 50% of all plants having EPDM membranes installed. Regular mechanical cleaning is the most common cleaning method followed by regular pressure release/air-bumping programs during operations. With regard to the diffuser pressure loss it was found that 50% of the evaluated plants had a diffuser pressure loss that was twice as high as measured for new diffusers.

A first-passage scheme for determination of overall rate constants for non-diffusion-limited suspensions

NASA Astrophysics Data System (ADS)

Lu, Shih-Yuan; Yen, Yi-Ming

2002-02-01

A first-passage scheme is devised to determine the overall rate constant of suspensions under the non-diffusion-limited condition. The original first-passage scheme developed for diffusion-limited processes is modified to account for the finite incorporation rate at the inclusion surface by using a concept of the nonzero survival probability of the diffusing entity at entity-inclusion encounters. This nonzero survival probability is obtained from solving a relevant boundary value problem. The new first-passage scheme is validated by an excellent agreement between overall rate constant results from the present development and from an accurate boundary collocation calculation for the three common spherical arrays [J. Chem. Phys. 109, 4985 (1998)], namely simple cubic, body-centered cubic, and face-centered cubic arrays, for a wide range of P and f. Here, P is a dimensionless quantity characterizing the relative rate of diffusion versus surface incorporation, and f is the volume fraction of the inclusion. The scheme is further applied to random spherical suspensions and to investigate the effect of inclusion coagulation on overall rate constants. It is found that randomness in inclusion arrangement tends to lower the overall rate constant for f up to the near close-packing value of the regular arrays because of the inclusion screening effect. This screening effect turns stronger for regular arrays when f is near and above the close-packing value of the regular arrays, and consequently the overall rate constant of the random array exceeds that of the regular array. Inclusion coagulation too induces the inclusion screening effect, and leads to lower overall rate constants.

Multicriteria relocation analysis of an off-site radioactive monitoring network for a nuclear power plant.

PubMed

Chang, Ni-Bin; Ning, Shu-Kuang; Chen, Jen-Chang

2006-08-01

Due to increasing environmental consciousness in most countries, every utility that owns a commercial nuclear power plant has been required to have both an on-site and off-site emergency response plan since the 1980s. A radiation monitoring network, viewed as part of the emergency response plan, can provide information regarding the radiation dosage emitted from a nuclear power plant in a regular operational period and/or abnormal measurements in an emergency event. Such monitoring information might help field operators and decision-makers to provide accurate responses or make decisions to protect the public health and safety. This study aims to conduct an integrated simulation and optimization analysis looking for the relocation strategy of a long-term regular off-site monitoring network at a nuclear power plant. The planning goal is to downsize the current monitoring network but maintain its monitoring capacity as much as possible. The monitoring sensors considered in this study include the thermoluminescence dosimetry (TLD) and air sampling system (AP) simultaneously. It is designed for detecting the radionuclide accumulative concentration, the frequency of violation, and the possible population affected by a long-term impact in the surrounding area regularly while it can also be used in an accidental release event. With the aid of the calibrated Industrial Source Complex-Plume Rise Model Enhancements (ISC-PRIME) simulation model to track down the possible radionuclide diffusion, dispersion, transport, and transformation process in the atmospheric environment, a multiobjective evaluation process can be applied to achieve the screening of monitoring stations for the nuclear power plant located at Hengchun Peninsula, South Taiwan. To account for multiple objectives, this study calculated preference weights to linearly combine objective functions leading to decision-making with exposure assessment in an optimization context. Final suggestions should be useful for narrowing the set of scenarios that decision-makers need to consider in this relocation process.

Classification Order of Surface-Confined Intermixing at Epitaxial Interface

NASA Astrophysics Data System (ADS)

Michailov, M.

The self-organization phenomena at epitaxial interface hold special attention in contemporary material science. Being relevant to the fundamental physical problem of competing, long-range and short-range atomic interactions in systems with reduced dimensionality, these phenomena have found exacting academic interest. They are also of great technological importance for their ability to bring spontaneous formation of regular nanoscale surface patterns and superlattices with exotic properties. The basic phenomenon involved in this process is surface diffusion. That is the motivation behind the present study which deals with important details of diffusion scenarios that control the fine atomic structure of epitaxial interface. Consisting surface imperfections (terraces, steps, kinks, and vacancies), the interface offers variety of barriers for surface diffusion. Therefore, the adatoms and clusters need a certain critical energy to overcome the corresponding diffusion barriers. In the most general case the critical energies can be attained by variation of the system temperature. Hence, their values define temperature limits of system energy gaps associated with different diffusion scenarios. This systematization imply classification order of surface alloying: blocked, incomplete, and complete. On that background, two diffusion problems, related to the atomic-scale surface morphology, will be discussed. The first problem deals with diffusion of atomic clusters on atomically smooth interface. On flat domains, far from terraces and steps, we analyzed the impact of size, shape, and cluster/substrate lattice misfit on the diffusion behavior of atomic clusters (islands). We found that the lattice constant of small clusters depends on the number N of building atoms at 1 < N ≤ 10. In heteroepitaxy, this effect of variable lattice constant originates from the enhanced charge transfer and the strong influence of the surface potential on cluster atomic arrangement. At constant temperature, the variation of the lattice constant leads to variable misfit which affects the island migration. The cluster/substrate commensurability influences the oscillation behavior of the diffusion coefficient caused by variation in the cluster shape. We discuss the results in a physical model that implies cluster diffusion with size-dependent cluster/substrate misfit. The second problem is devoted to diffusion phenomena in the vicinity of atomic terraces on stepped or vicinal surfaces. Here, we develop a computational model that refines important details of diffusion behavior of adatoms accounting for the energy barriers at specific atomic sites (smooth domains, terraces, and steps) located on the crystal surface. The dynamic competition between energy gained by mixing and substrate strain energy results in diffusion scenario where adatoms form alloyed islands and alloyed stripes in the vicinity of terrace edges. Being in agreement with recent experimental findings, the observed effect of stripe and island alloy formation opens up a way regular surface patterns to be configured at different atomic levels on the crystal surface. The complete surface alloying of the entire interface layer is also briefly discussed with critical analysis and classification of experimental findings and simulation data.

Retaining both discrete and smooth features in 1D and 2D NMR relaxation and diffusion experiments

NASA Astrophysics Data System (ADS)

Reci, A.; Sederman, A. J.; Gladden, L. F.

2017-11-01

A new method of regularization of 1D and 2D NMR relaxation and diffusion experiments is proposed and a robust algorithm for its implementation is introduced. The new form of regularization, termed the Modified Total Generalized Variation (MTGV) regularization, offers a compromise between distinguishing discrete and smooth features in the reconstructed distributions. The method is compared to the conventional method of Tikhonov regularization and the recently proposed method of L1 regularization, when applied to simulated data of 1D spin-lattice relaxation, T1, 1D spin-spin relaxation, T2, and 2D T1-T2 NMR experiments. A range of simulated distributions composed of two lognormally distributed peaks were studied. The distributions differed with regard to the variance of the peaks, which were designed to investigate a range of distributions containing only discrete, only smooth or both features in the same distribution. Three different signal-to-noise ratios were studied: 2000, 200 and 20. A new metric is proposed to compare the distributions reconstructed from the different regularization methods with the true distributions. The metric is designed to penalise reconstructed distributions which show artefact peaks. Based on this metric, MTGV regularization performs better than Tikhonov and L1 regularization in all cases except when the distribution is known to only comprise of discrete peaks, in which case L1 regularization is slightly more accurate than MTGV regularization.

On the local well-posedness and a Prodi-Serrin-type regularity criterion of the three-dimensional MHD-Boussinesq system without thermal diffusion

NASA Astrophysics Data System (ADS)

Larios, Adam; Pei, Yuan

2017-07-01

We prove a Prodi-Serrin-type global regularity condition for the three-dimensional Magnetohydrodynamic-Boussinesq system (3D MHD-Boussinesq) without thermal diffusion, in terms of only two velocity and two magnetic components. To the best of our knowledge, this is the first Prodi-Serrin-type criterion for such a 3D hydrodynamic system which is not fully dissipative, and indicates that such an approach may be successful on other systems. In addition, we provide a constructive proof of the local well-posedness of solutions to the fully dissipative 3D MHD-Boussinesq system, and also the fully inviscid, irresistive, non-diffusive MHD-Boussinesq equations. We note that, as a special case, these results include the 3D non-diffusive Boussinesq system and the 3D MHD equations. Moreover, they can be extended without difficulty to include the case of a Coriolis rotational term.

Combined Tensor Fitting and TV Regularization in Diffusion Tensor Imaging Based on a Riemannian Manifold Approach.

PubMed

Baust, Maximilian; Weinmann, Andreas; Wieczorek, Matthias; Lasser, Tobias; Storath, Martin; Navab, Nassir

2016-08-01

In this paper, we consider combined TV denoising and diffusion tensor fitting in DTI using the affine-invariant Riemannian metric on the space of diffusion tensors. Instead of first fitting the diffusion tensors, and then denoising them, we define a suitable TV type energy functional which incorporates the measured DWIs (using an inverse problem setup) and which measures the nearness of neighboring tensors in the manifold. To approach this functional, we propose generalized forward- backward splitting algorithms which combine an explicit and several implicit steps performed on a decomposition of the functional. We validate the performance of the derived algorithms on synthetic and real DTI data. In particular, we work on real 3D data. To our knowledge, the present paper describes the first approach to TV regularization in a combined manifold and inverse problem setup.

Turing instability in reaction-diffusion models on complex networks

NASA Astrophysics Data System (ADS)

Ide, Yusuke; Izuhara, Hirofumi; Machida, Takuya

2016-09-01

In this paper, the Turing instability in reaction-diffusion models defined on complex networks is studied. Here, we focus on three types of models which generate complex networks, i.e. the Erdős-Rényi, the Watts-Strogatz, and the threshold network models. From analysis of the Laplacian matrices of graphs generated by these models, we numerically reveal that stable and unstable regions of a homogeneous steady state on the parameter space of two diffusion coefficients completely differ, depending on the network architecture. In addition, we theoretically discuss the stable and unstable regions in the cases of regular enhanced ring lattices which include regular circles, and networks generated by the threshold network model when the number of vertices is large enough.

Fem Simulation of Triple Diffusive Natural Convection Along Inclined Plate in Porous Medium: Prescribed Surface Heat, Solute and Nanoparticles Flux

NASA Astrophysics Data System (ADS)

Goyal, M.; Goyal, R.; Bhargava, R.

2017-12-01

In this paper, triple diffusive natural convection under Darcy flow over an inclined plate embedded in a porous medium saturated with a binary base fluid containing nanoparticles and two salts is studied. The model used for the nanofluid is the one which incorporates the effects of Brownian motion and thermophoresis. In addition, the thermal energy equations include regular diffusion and cross-diffusion terms. The vertical surface has the heat, mass and nanoparticle fluxes each prescribed as a power law function of the distance along the wall. The boundary layer equations are transformed into a set of ordinary differential equations with the help of group theory transformations. A wide range of parameter values are chosen to bring out the effect of buoyancy ratio, regular Lewis number and modified Dufour parameters of both salts and nanofluid parameters with varying angle of inclinations. The effects of parameters on the velocity, temperature, solutal and nanoparticles volume fraction profiles, as well as on the important parameters of heat and mass transfer, i.e., the reduced Nusselt, regular and nanofluid Sherwood numbers, are discussed. Such problems find application in extrusion of metals, polymers and ceramics, production of plastic films, insulation of wires and liquid packaging.

On regularization and error estimates for the backward heat conduction problem with time-dependent thermal diffusivity factor

NASA Astrophysics Data System (ADS)

Karimi, Milad; Moradlou, Fridoun; Hajipour, Mojtaba

2018-10-01

This paper is concerned with a backward heat conduction problem with time-dependent thermal diffusivity factor in an infinite "strip". This problem is drastically ill-posed which is caused by the amplified infinitely growth in the frequency components. A new regularization method based on the Meyer wavelet technique is developed to solve the considered problem. Using the Meyer wavelet technique, some new stable estimates are proposed in the Hölder and Logarithmic types which are optimal in the sense of given by Tautenhahn. The stability and convergence rate of the proposed regularization technique are proved. The good performance and the high-accuracy of this technique is demonstrated through various one and two dimensional examples. Numerical simulations and some comparative results are presented.

Facilitated diffusion in chromatin lattices: mechanistic diversity and regulatory potential.

PubMed

Kampmann, Martin

2005-08-01

The interaction between a protein and a specific DNA site is the molecular basis for vital processes in all organisms. Location of the DNA target site by the protein commonly involves facilitated diffusion. Mechanisms of facilitated diffusion vary among proteins; they include one- and two-dimensional sliding along DNA, direct transfer between uncorrelated sites, as well as combinations of these mechanisms. Facilitated diffusion has almost exclusively been studied in vitro. This review discusses facilitated diffusion in the context of the living cell and proposes a theoretical model for facilitated diffusion in chromatin lattices. Chromatin structure differentially affects proteins in different modes of diffusion. The interplay of facilitated diffusion and chromatin structure can determine the rate of protein association with the target site, the frequency of association-dissociation events at the target site, and, under particular conditions, the occupancy of the target site. Facilitated diffusion is required in vivo for efficient DNA repair and bacteriophage restriction and has potential roles in fine-tuning gene regulatory networks and kinetically compartmentalizing the eukaryotic nucleus.

Diffusion of Sites versus Polymers in Polyelectrolyte Complexes and Multilayers.

PubMed

Fares, Hadi M; Schlenoff, Joseph B

2017-10-18

It has long been assumed that the spontaneous formation of materials such as complexes and multilayers from charged polymers depends on (inter)diffusion of these polyelectrolytes. Here, we separately examine the mass transport of polymer molecules and extrinsic sites-charged polyelectrolyte repeat units balanced by counterions-within thin films of polyelectrolyte complex, PEC, using sensitive isotopic labeling techniques. The apparent diffusion coefficients of these sites within PEC films of poly(diallyldimethylammonium), PDADMA, and poly(styrenesulfonate), PSS, are at least 2 orders of magnitude faster than the diffusion of polyelectrolytes themselves. This is because site diffusion requires only local rearrangements of polyelectrolyte repeat units, placing far fewer kinetic limitations on the assembly of polyelectrolyte complexes in all of their forms. Site diffusion strongly depends on the salt concentration (ionic strength) of the environment, and diffusion of PDADMA sites is faster than that of PSS sites, accounting for the asymmetric nature of multilayer growth. Site diffusion is responsible for multilayer growth in the linear and into the exponential regimes, which explains how PDADMA can mysteriously "pass through" layers of PSS. Using quantitative relationships between site diffusion coefficient and salt concentration, conditions were identified that allowed the diffusion length to always exceed the film thickness, leading to full exponential growth over 3 orders of magnitude thickness. Both site and polymer diffusion were independent of molecular weight, suggesting that ion pairing density is a limiting factor. Polyelectrolyte complexes are examples of a broader class of dynamic bulk polymeric materials that (self-) assemble via the transport of cross-links or defects rather than actual molecules.

Detecting chaos, determining the dimensions of tori and predicting slow diffusion in Fermi-Pasta-Ulam lattices by the Generalized Alignment Index method

NASA Astrophysics Data System (ADS)

Skokos, C.; Bountis, T.; Antonopoulos, C.

2008-12-01

The recently introduced GALI method is used for rapidly detecting chaos, determining the dimensionality of regular motion and predicting slow diffusion in multi-dimensional Hamiltonian systems. We propose an efficient computation of the GALIk indices, which represent volume elements of k randomly chosen deviation vectors from a given orbit, based on the Singular Value Decomposition (SVD) algorithm. We obtain theoretically and verify numerically asymptotic estimates of GALIs long-time behavior in the case of regular orbits lying on low-dimensional tori. The GALIk indices are applied to rapidly detect chaotic oscillations, identify low-dimensional tori of Fermi-Pasta-Ulam (FPU) lattices at low energies and predict weak diffusion away from quasiperiodic motion, long before it is actually observed in the oscillations.

A combined reconstruction-classification method for diffuse optical tomography.

PubMed

Hiltunen, P; Prince, S J D; Arridge, S

2009-11-07

We present a combined classification and reconstruction algorithm for diffuse optical tomography (DOT). DOT is a nonlinear ill-posed inverse problem. Therefore, some regularization is needed. We present a mixture of Gaussians prior, which regularizes the DOT reconstruction step. During each iteration, the parameters of a mixture model are estimated. These associate each reconstructed pixel with one of several classes based on the current estimate of the optical parameters. This classification is exploited to form a new prior distribution to regularize the reconstruction step and update the optical parameters. The algorithm can be described as an iteration between an optimization scheme with zeroth-order variable mean and variance Tikhonov regularization and an expectation-maximization scheme for estimation of the model parameters. We describe the algorithm in a general Bayesian framework. Results from simulated test cases and phantom measurements show that the algorithm enhances the contrast of the reconstructed images with good spatial accuracy. The probabilistic classifications of each image contain only a few misclassified pixels.

Effects of high-frequency damping on iterative convergence of implicit viscous solver

NASA Astrophysics Data System (ADS)

Nishikawa, Hiroaki; Nakashima, Yoshitaka; Watanabe, Norihiko

2017-11-01

This paper discusses effects of high-frequency damping on iterative convergence of an implicit defect-correction solver for viscous problems. The study targets a finite-volume discretization with a one parameter family of damped viscous schemes. The parameter α controls high-frequency damping: zero damping with α = 0, and larger damping for larger α (> 0). Convergence rates are predicted for a model diffusion equation by a Fourier analysis over a practical range of α. It is shown that the convergence rate attains its minimum at α = 1 on regular quadrilateral grids, and deteriorates for larger values of α. A similar behavior is observed for regular triangular grids. In both quadrilateral and triangular grids, the solver is predicted to diverge for α smaller than approximately 0.5. Numerical results are shown for the diffusion equation and the Navier-Stokes equations on regular and irregular grids. The study suggests that α = 1 and 4/3 are suitable values for robust and efficient computations, and α = 4 / 3 is recommended for the diffusion equation, which achieves higher-order accuracy on regular quadrilateral grids. Finally, a Jacobian-Free Newton-Krylov solver with the implicit solver (a low-order Jacobian approximately inverted by a multi-color Gauss-Seidel relaxation scheme) used as a variable preconditioner is recommended for practical computations, which provides robust and efficient convergence for a wide range of α.

Determination of heat transfer parameters by use of finite integral transform and experimental data for regular geometric shapes

NASA Astrophysics Data System (ADS)

Talaghat, Mohammad Reza; Jokar, Seyyed Mohammad

2017-12-01

This article offers a study on estimation of heat transfer parameters (coefficient and thermal diffusivity) using analytical solutions and experimental data for regular geometric shapes (such as infinite slab, infinite cylinder, and sphere). Analytical solutions have a broad use in experimentally determining these parameters. Here, the method of Finite Integral Transform (FIT) was used for solutions of governing differential equations. The temperature change at centerline location of regular shapes was recorded to determine both the thermal diffusivity and heat transfer coefficient. Aluminum and brass were used for testing. Experiments were performed for different conditions such as in a highly agitated water medium ( T = 52 °C) and in air medium ( T = 25 °C). Then, with the known slope of the temperature ratio vs. time curve and thickness of slab or radius of the cylindrical or spherical materials, thermal diffusivity value and heat transfer coefficient may be determined. According to the method presented in this study, the estimated of thermal diffusivity of aluminum and brass is 8.395 × 10-5 and 3.42 × 10-5 for a slab, 8.367 × 10-5 and 3.41 × 10-5 for a cylindrical rod and 8.385 × 10-5 and 3.40 × 10-5 m2/s for a spherical shape, respectively. The results showed there is close agreement between the values estimated here and those already published in the literature. The TAAD% is 0.42 and 0.39 for thermal diffusivity of aluminum and brass, respectively.

The IMADA-AVER Boundary Layer Experiment in the Mexico City Area.

NASA Astrophysics Data System (ADS)

Doran, J. C.; Bian, X.; de Wekker, S. F. J.; Edgerton, S.; Fast, J. D.; Hubbe, J. M.; Shaw, W. J.; Whiteman, C. D.; Abbott, S.; King, C.; Leach, J.; Mulhearn, M.; Russell, C.; Templeman, B.; Wolfe, D.; Archuleta, J.; Elliott, S.; Fernandez, A.; Langley, D.; Lee, J. T.; Porch, W.; Tellier, L.; Chow, J.; Watson, J. G.; Coulter, R. L.; Martin, T. J.; Shannon, J. D.; White, R.; Martinez, D.; Martinez, J. L.; Mora, V.; Sosa, G.; Mercado, G.; Pena, J. L.; Salas, R.; Petty, R.

1998-11-01

A boundary layer field experiment in the Mexico City basin during the period 24 February-22 March 1997 is described. A total of six sites were instrumented. At four of the sites, 915-MHz radar wind profilers were deployed and radiosondes were released five times per day. Two of these sites also had sodars collocated with the profilers. Radiosondes were released twice per day at a fifth site to the south of the basin, and rawinsondes were flown from another location to the northeast of the city three times per day. Mixed layers grew to depths of 2500-3500 m, with a rapid period of growth beginning shortly before noon and lasting for several hours. Significant differences between the mixed-layer temperatures in the basin and outside the basin were observed. Three thermally and topographically driven flow patterns were observed that are consistent with previously hypothesized topographical and thermal forcing mechanisms. Despite these features, the circulation patterns in the basin important for the transport and diffusion of air pollutants show less day-to-day regularity than had been anticipated on the basis of Mexico City's tropical location, high altitude and strong insolation, and topographical setting.

UV spectral irradiance measurements in New Zealand: Effects of Pinatubo volcanic aerosol

NASA Technical Reports Server (NTRS)

Mckenzie, Richard L.

1994-01-01

Since late 1989, regular UV spectral irradiance measurements have been made at Lauder, New Zealand (45 deg S, 170 deg E), whenever weather permits. Here, the instrumentation and measurement strategy are outlined, and early results are discussed. Following the eruption of Mt Pinatubo in June 1991, large amounts of volcanic aerosol were injected into the stratosphere and were subsequently transported to New Zealand's latitudes in the latter half of 1991. This provides an opportunity to investigate the effects of volcanic aerosols on UV irradiances measured at this clean-air site. Although changes in global (sum of diffuse plus direct) irradiances were below the detection threshold, there were significant changes in the partitioning of radiation between the direct beam and diffuse skylight. Decreases by nearly a factor of two in the direct/diffuse ratio were observed at longer wavelengths, and at smaller solar zenith angles (sza's). The aerosol optical depth due to volcanic aerosol over Lauder in December 1991 was 0.15 plus or minus 0.02 at 450 nm, with lower values at shorter wavelengths. Although effects were relatively small in the UVB region, an implication of the changes is that the contrast between shade and direct sun is reduced, so that shaded areas received relatively more radiation in the summer of 1991/92 in New Zealand.

Lung injury induced by secondhand smoke exposure detected with hyperpolarized helium-3 diffusion MR.

PubMed

Wang, Chengbo; Mugler, John P; de Lange, Eduard E; Patrie, James T; Mata, Jaime F; Altes, Talissa A

2014-01-01

To determine whether helium-3 diffusion MR can detect the changes in the lungs of healthy nonsmoking individuals who were regularly exposed to secondhand smoke. Three groups were studied (age: 59 ± 9 years): 23 smokers, 37 exposure-to-secondhand-smoke subjects, and 29 control subjects. We measured helium-3 diffusion values at diffusion times from 0.23 to 1.97 s. One-way analysis of variance revealed that the mean area under the helium-3 diffusion curves (ADC AUC) of the smokers was significantly elevated compared with the controls and to the exposure-to-secondhand-smoke subjects (P < 0.001 both). No difference between the mean ADC AUC of the exposure-to-secondhand-smoke subjects and that of the controls was found (P = 0.115). However, application of a receiver operator characteristic-derived rule to classify subjects as either a "control" or a "smoker," based on ADC AUC, revealed that 30% (11/37) of the exposure-to-secondhand subjects were classified as "smokers" indicating an elevation of the ADC AUC. Using helium-3 diffusion MR, elevated ADC values were detected in 30% of nonsmoking healthy subjects who had been regularly exposed to secondhand smoke, supporting the concept that, in susceptible individuals, secondhand smoke causes mild lung damage. Copyright © 2013 Wiley Periodicals, Inc.

Template-directed fabrication of porous gas diffusion layer for magnesium air batteries

NASA Astrophysics Data System (ADS)

Xue, Yejian; Miao, He; Sun, Shanshan; Wang, Qin; Li, Shihua; Liu, Zhaoping

2015-11-01

The uniform micropore distribution in the gas diffusion layers (GDLs) of the air-breathing cathode is very important for the metal air batteries. In this work, the super-hydrophobic GDL with the interconnected regular pores is prepared by a facile silica template method, and then the electrochemical properties of the Mg air batteries containing these GDLs are investigated. The results indicate that the interconnected and uniform pore structure, the available water-breakout pressure and the high gas permeability coefficient of the GDL can be obtained by the application of 30% silica template. The maximum power density of the Mg air battery containing the GDL with 30% regular pores reaches 88.9 mW cm-2 which is about 1.2 times that containing the pristine GDL. Furthermore, the GDL with 30% regular pores exhibits the improved the long term hydrophobic stability.

q-Space Upsampling Using x-q Space Regularization.

PubMed

Chen, Geng; Dong, Bin; Zhang, Yong; Shen, Dinggang; Yap, Pew-Thian

2017-09-01

Acquisition time in diffusion MRI increases with the number of diffusion-weighted images that need to be acquired. Particularly in clinical settings, scan time is limited and only a sparse coverage of the vast q -space is possible. In this paper, we show how non-local self-similar information in the x - q space of diffusion MRI data can be harnessed for q -space upsampling. More specifically, we establish the relationships between signal measurements in x - q space using a patch matching mechanism that caters to unstructured data. We then encode these relationships in a graph and use it to regularize an inverse problem associated with recovering a high q -space resolution dataset from its low-resolution counterpart. Experimental results indicate that the high-resolution datasets reconstructed using the proposed method exhibit greater quality, both quantitatively and qualitatively, than those obtained using conventional methods, such as interpolation using spherical radial basis functions (SRBFs).

Obtaining sparse distributions in 2D inverse problems.

PubMed

Reci, A; Sederman, A J; Gladden, L F

2017-08-01

The mathematics of inverse problems has relevance across numerous estimation problems in science and engineering. L 1 regularization has attracted recent attention in reconstructing the system properties in the case of sparse inverse problems; i.e., when the true property sought is not adequately described by a continuous distribution, in particular in Compressed Sensing image reconstruction. In this work, we focus on the application of L 1 regularization to a class of inverse problems; relaxation-relaxation, T 1 -T 2 , and diffusion-relaxation, D-T 2 , correlation experiments in NMR, which have found widespread applications in a number of areas including probing surface interactions in catalysis and characterizing fluid composition and pore structures in rocks. We introduce a robust algorithm for solving the L 1 regularization problem and provide a guide to implementing it, including the choice of the amount of regularization used and the assignment of error estimates. We then show experimentally that L 1 regularization has significant advantages over both the Non-Negative Least Squares (NNLS) algorithm and Tikhonov regularization. It is shown that the L 1 regularization algorithm stably recovers a distribution at a signal to noise ratio<20 and that it resolves relaxation time constants and diffusion coefficients differing by as little as 10%. The enhanced resolving capability is used to measure the inter and intra particle concentrations of a mixture of hexane and dodecane present within porous silica beads immersed within a bulk liquid phase; neither NNLS nor Tikhonov regularization are able to provide this resolution. This experimental study shows that the approach enables discrimination between different chemical species when direct spectroscopic discrimination is impossible, and hence measurement of chemical composition within porous media, such as catalysts or rocks, is possible while still being stable to high levels of noise. Copyright © 2017. Published by Elsevier Inc.

Obtaining sparse distributions in 2D inverse problems

NASA Astrophysics Data System (ADS)

Reci, A.; Sederman, A. J.; Gladden, L. F.

2017-08-01

The mathematics of inverse problems has relevance across numerous estimation problems in science and engineering. L1 regularization has attracted recent attention in reconstructing the system properties in the case of sparse inverse problems; i.e., when the true property sought is not adequately described by a continuous distribution, in particular in Compressed Sensing image reconstruction. In this work, we focus on the application of L1 regularization to a class of inverse problems; relaxation-relaxation, T1-T2, and diffusion-relaxation, D-T2, correlation experiments in NMR, which have found widespread applications in a number of areas including probing surface interactions in catalysis and characterizing fluid composition and pore structures in rocks. We introduce a robust algorithm for solving the L1 regularization problem and provide a guide to implementing it, including the choice of the amount of regularization used and the assignment of error estimates. We then show experimentally that L1 regularization has significant advantages over both the Non-Negative Least Squares (NNLS) algorithm and Tikhonov regularization. It is shown that the L1 regularization algorithm stably recovers a distribution at a signal to noise ratio < 20 and that it resolves relaxation time constants and diffusion coefficients differing by as little as 10%. The enhanced resolving capability is used to measure the inter and intra particle concentrations of a mixture of hexane and dodecane present within porous silica beads immersed within a bulk liquid phase; neither NNLS nor Tikhonov regularization are able to provide this resolution. This experimental study shows that the approach enables discrimination between different chemical species when direct spectroscopic discrimination is impossible, and hence measurement of chemical composition within porous media, such as catalysts or rocks, is possible while still being stable to high levels of noise.

Radon Diffusion Measurement in Polyethylene based on Alpha Detection

NASA Astrophysics Data System (ADS)

Rau, Wolfgang

2011-04-01

We present a method to measure the diffusion of Radon in solid materials based on the alpha decay of the radon daughter products. In contrast to usual diffusion measurements which detect the radon that penetrates a thin barrier, we let the radon diffuse into the material and then measure the alpha decays of the radon daughter products in the material. We applied this method to regular and ultra high molecular weight poly ethylene and find diffusion lengths of order of mm as expected. However, the preliminary analysis shows significant differences between two different approaches we have chosen. These differences may be explained by the different experimental conditions.

Efficiency analysis of diffusion on T-fractals in the sense of random walks.

PubMed

Peng, Junhao; Xu, Guoai

2014-04-07

Efficiently controlling the diffusion process is crucial in the study of diffusion problem in complex systems. In the sense of random walks with a single trap, mean trapping time (MTT) and mean diffusing time (MDT) are good measures of trapping efficiency and diffusion efficiency, respectively. They both vary with the location of the node. In this paper, we analyze the effects of node's location on trapping efficiency and diffusion efficiency of T-fractals measured by MTT and MDT. First, we provide methods to calculate the MTT for any target node and the MDT for any source node of T-fractals. The methods can also be used to calculate the mean first-passage time between any pair of nodes. Then, using the MTT and the MDT as the measure of trapping efficiency and diffusion efficiency, respectively, we compare the trapping efficiency and diffusion efficiency among all nodes of T-fractal and find the best (or worst) trapping sites and the best (or worst) diffusing sites. Our results show that the hub node of T-fractal is the best trapping site, but it is also the worst diffusing site; and that the three boundary nodes are the worst trapping sites, but they are also the best diffusing sites. Comparing the maximum of MTT and MDT with their minimums, we find that the maximum of MTT is almost 6 times of the minimum of MTT and the maximum of MDT is almost equal to the minimum for MDT. Thus, the location of target node has large effect on the trapping efficiency, but the location of source node almost has no effect on diffusion efficiency. We also simulate random walks on T-fractals, whose results are consistent with the derived results.

Statistical analysis of nonlinearly reconstructed near-infrared tomographic images: Part I--Theory and simulations.

PubMed

Pogue, Brian W; Song, Xiaomei; Tosteson, Tor D; McBride, Troy O; Jiang, Shudong; Paulsen, Keith D

2002-07-01

Near-infrared (NIR) diffuse tomography is an emerging method for imaging the interior of tissues to quantify concentrations of hemoglobin and exogenous chromophores non-invasively in vivo. It often exploits an optical diffusion model-based image reconstruction algorithm to estimate spatial property values from measurements of the light flux at the surface of the tissue. In this study, mean-squared error (MSE) over the image is used to evaluate methods for regularizing the ill-posed inverse image reconstruction problem in NIR tomography. Estimates of image bias and image standard deviation were calculated based upon 100 repeated reconstructions of a test image with randomly distributed noise added to the light flux measurements. It was observed that the bias error dominates at high regularization parameter values while variance dominates as the algorithm is allowed to approach the optimal solution. This optimum does not necessarily correspond to the minimum projection error solution, but typically requires further iteration with a decreasing regularization parameter to reach the lowest image error. Increasing measurement noise causes a need to constrain the minimum regularization parameter to higher values in order to achieve a minimum in the overall image MSE.

A Comparison Between Modeled and Measured Clear-Sky Radiative Shortwave Fluxes in Arctic Environments, with Special Emphasis on Diffuse Radiation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Barnard, James C.; Flynn, Donna M.

2002-10-08

The ability of the SBDART radiative transfer model to predict clear-sky diffuse and direct normal broadband shortwave irradiances is investigated. Model calculations of these quantities are compared with data from the Atmospheric Radiation Measurement (ARM) program’s Southern Great Plains (SGP) and North Slope of Alaska (NSA) sites. The model tends to consistently underestimate the direct normal irradiances at both sites by about 1%. In regards to clear-sky diffuse irradiance, the model overestimates this quantity at the SGP site in a manner similar to what has been observed in other studies (Halthore and Schwartz, 2000). The difference between the diffuse SBDARTmore » calculations and Halthore and Schwartz’s MODTRAN calculations is very small, thus demonstrating that SBDART performs similarly to MODTRAN. SBDART is then applied to the NSA site, and here it is found that the discrepancy between the model calculations and corrected diffuse measurements (corrected for daytime offsets, Dutton et al., 2001) is 0.4 W/m2 when averaged over the 12 cases considered here. Two cases of diffuse measurements from a shaded “black and white” pyranometer are also compared with the calculations and the discrepancy is again minimal. Thus, it appears as if the “diffuse discrepancy” that exists at the SGP site does not exist at the NSA sites. We cannot yet explain why the model predicts diffuse radiation well at one site but not at the other.« less

Electrochemical Impedance Imaging via the Distribution of Diffusion Times

NASA Astrophysics Data System (ADS)

Song, Juhyun; Bazant, Martin Z.

2018-03-01

We develop a mathematical framework to analyze electrochemical impedance spectra in terms of a distribution of diffusion times (DDT) for a parallel array of random finite-length Warburg (diffusion) or Gerischer (reaction-diffusion) circuit elements. A robust DDT inversion method is presented based on complex nonlinear least squares regression with Tikhonov regularization and illustrated for three cases of nanostructured electrodes for energy conversion: (i) a carbon nanotube supercapacitor, (ii) a silicon nanowire Li-ion battery, and (iii) a porous-carbon vanadium flow battery. The results demonstrate the feasibility of nondestructive "impedance imaging" to infer microstructural statistics of random, heterogeneous materials.

Diffuse optical correlation tomography of cerebral blood flow during cortical spreading depression in rat brain

NASA Astrophysics Data System (ADS)

Zhou, Chao; Yu, Guoqiang; Furuya, Daisuke; Greenberg, Joel; Yodh, Arjun; Durduran, Turgut

2006-02-01

Diffuse optical correlation methods were adapted for three-dimensional (3D) tomography of cerebral blood flow (CBF) in small animal models. The image reconstruction was optimized using a noise model for diffuse correlation tomography which enabled better data selection and regularization. The tomographic approach was demonstrated with simulated data and during in-vivo cortical spreading depression (CSD) in rat brain. Three-dimensional images of CBF were obtained through intact skull in tissues(~4mm) deep below the cortex.

Multi-site study of diffusion metric variability: effects of site, vendor, field strength, and echo time on regions-of-interest and histogram-bin analyses.

PubMed

Helmer, K G; Chou, M-C; Preciado, R I; Gimi, B; Rollins, N K; Song, A; Turner, J; Mori, S

2016-02-27

It is now common for magnetic-resonance-imaging (MRI) based multi-site trials to include diffusion-weighted imaging (DWI) as part of the protocol. It is also common for these sites to possess MR scanners of different manufacturers, different software and hardware, and different software licenses. These differences mean that scanners may not be able to acquire data with the same number of gradient amplitude values and number of available gradient directions. Variability can also occur in achievable b-values and minimum echo times. The challenge of a multi-site study then, is to create a common protocol by understanding and then minimizing the effects of scanner variability and identifying reliable and accurate diffusion metrics. This study describes the effect of site, scanner vendor, field strength, and TE on two diffusion metrics: the first moment of the diffusion tensor field (mean diffusivity, MD), and the fractional anisotropy (FA) using two common analyses (region-of-interest and mean-bin value of whole brain histograms). The goal of the study was to identify sources of variability in diffusion-sensitized imaging and their influence on commonly reported metrics. The results demonstrate that the site, vendor, field strength, and echo time all contribute to variability in FA and MD, though to different extent. We conclude that characterization of the variability of DTI metrics due to site, vendor, field strength, and echo time is a worthwhile step in the construction of multi-center trials.

Novel harmonic regularization approach for variable selection in Cox's proportional hazards model.

PubMed

Chu, Ge-Jin; Liang, Yong; Wang, Jia-Xuan

2014-01-01

Variable selection is an important issue in regression and a number of variable selection methods have been proposed involving nonconvex penalty functions. In this paper, we investigate a novel harmonic regularization method, which can approximate nonconvex Lq  (1/2 < q < 1) regularizations, to select key risk factors in the Cox's proportional hazards model using microarray gene expression data. The harmonic regularization method can be efficiently solved using our proposed direct path seeking approach, which can produce solutions that closely approximate those for the convex loss function and the nonconvex regularization. Simulation results based on the artificial datasets and four real microarray gene expression datasets, such as real diffuse large B-cell lymphoma (DCBCL), the lung cancer, and the AML datasets, show that the harmonic regularization method can be more accurate for variable selection than existing Lasso series methods.

Regularity of random attractors for fractional stochastic reaction-diffusion equations on Rn

NASA Astrophysics Data System (ADS)

Gu, Anhui; Li, Dingshi; Wang, Bixiang; Yang, Han

2018-06-01

We investigate the regularity of random attractors for the non-autonomous non-local fractional stochastic reaction-diffusion equations in Hs (Rn) with s ∈ (0 , 1). We prove the existence and uniqueness of the tempered random attractor that is compact in Hs (Rn) and attracts all tempered random subsets of L2 (Rn) with respect to the norm of Hs (Rn). The main difficulty is to show the pullback asymptotic compactness of solutions in Hs (Rn) due to the noncompactness of Sobolev embeddings on unbounded domains and the almost sure nondifferentiability of the sample paths of the Wiener process. We establish such compactness by the ideas of uniform tail-estimates and the spectral decomposition of solutions in bounded domains.

Field-scale forward and back diffusion through low-permeability zones

NASA Astrophysics Data System (ADS)

Yang, Minjune; Annable, Michael D.; Jawitz, James W.

2017-07-01

Understanding the effects of back diffusion of groundwater contaminants from low-permeability zones to aquifers is critical to making site management decisions related to remedial actions. Here, we combine aquifer and aquitard data to develop recommended site characterization strategies using a three-stage classification of plume life cycle based on the solute origins: aquifer source zone dissolution, source zone dissolution combined with back diffusion from an aquitard, and only back diffusion. We use measured aquitard concentration profile data from three field sites to identify signature shapes that are characteristic of these three stages. We find good fits to the measured data with analytical solutions that include the effects of advection and forward and back diffusion through low-permeability zones, and linearly and exponentially decreasing flux resulting from source dissolution in the aquifer. Aquifer contaminant time series data at monitoring wells from a mature site were well described using analytical solutions representing the combined case of source zone and back diffusion, while data from a site where the source had been isolated were well described solely by back diffusion. The modeling approach presented in this study is designed to enable site managers to implement appropriate remediation technologies at a proper timing for high- and low-permeability zones, considering estimated plume life cycle.

Field-scale forward and back diffusion through low-permeability zones.

PubMed

Yang, Minjune; Annable, Michael D; Jawitz, James W

2017-07-01

Understanding the effects of back diffusion of groundwater contaminants from low-permeability zones to aquifers is critical to making site management decisions related to remedial actions. Here, we combine aquifer and aquitard data to develop recommended site characterization strategies using a three-stage classification of plume life cycle based on the solute origins: aquifer source zone dissolution, source zone dissolution combined with back diffusion from an aquitard, and only back diffusion. We use measured aquitard concentration profile data from three field sites to identify signature shapes that are characteristic of these three stages. We find good fits to the measured data with analytical solutions that include the effects of advection and forward and back diffusion through low-permeability zones, and linearly and exponentially decreasing flux resulting from source dissolution in the aquifer. Aquifer contaminant time series data at monitoring wells from a mature site were well described using analytical solutions representing the combined case of source zone and back diffusion, while data from a site where the source had been isolated were well described solely by back diffusion. The modeling approach presented in this study is designed to enable site managers to implement appropriate remediation technologies at a proper timing for high- and low-permeability zones, considering estimated plume life cycle. Copyright © 2017 Elsevier B.V. All rights reserved.

Cooperation among cancer cells as public goods games on Voronoi networks.

PubMed

Archetti, Marco

2016-05-07

Cancer cells produce growth factors that diffuse and sustain tumour proliferation, a form of cooperation that can be studied using mathematical models of public goods in the framework of evolutionary game theory. Cell populations, however, form heterogeneous networks that cannot be described by regular lattices or scale-free networks, the types of graphs generally used in the study of cooperation. To describe the dynamics of growth factor production in populations of cancer cells, I study public goods games on Voronoi networks, using a range of non-linear benefits that account for the known properties of growth factors, and different types of diffusion gradients. The results are surprisingly similar to those obtained on regular graphs and different from results on scale-free networks, revealing that network heterogeneity per se does not promote cooperation when public goods diffuse beyond one-step neighbours. The exact shape of the diffusion gradient is not crucial, however, whereas the type of non-linear benefit is an essential determinant of the dynamics. Public goods games on Voronoi networks can shed light on intra-tumour heterogeneity, the evolution of resistance to therapies that target growth factors, and new types of cell therapy. Copyright © 2016 Elsevier Ltd. All rights reserved.

Massively Parallel Simulations of Diffusion in Dense Polymeric Structures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Faulon, Jean-Loup, Wilcox, R.T.

1997-11-01

An original computational technique to generate close-to-equilibrium dense polymeric structures is proposed. Diffusion of small gases are studied on the equilibrated structures using massively parallel molecular dynamics simulations running on the Intel Teraflops (9216 Pentium Pro processors) and Intel Paragon(1840 processors). Compared to the current state-of-the-art equilibration methods this new technique appears to be faster by some orders of magnitude.The main advantage of the technique is that one can circumvent the bottlenecks in configuration space that inhibit relaxation in molecular dynamics simulations. The technique is based on the fact that tetravalent atoms (such as carbon and silicon) fit in themore » center of a regular tetrahedron and that regular tetrahedrons can be used to mesh the three-dimensional space. Thus, the problem of polymer equilibration described by continuous equations in molecular dynamics is reduced to a discrete problem where solutions are approximated by simple algorithms. Practical modeling applications include the constructing of butyl rubber and ethylene-propylene-dimer-monomer (EPDM) models for oxygen and water diffusion calculations. Butyl and EPDM are used in O-ring systems and serve as sealing joints in many manufactured objects. Diffusion coefficients of small gases have been measured experimentally on both polymeric systems, and in general the diffusion coefficients in EPDM are an order of magnitude larger than in butyl. In order to better understand the diffusion phenomena, 10, 000 atoms models were generated and equilibrated for butyl and EPDM. The models were submitted to a massively parallel molecular dynamics simulation to monitor the trajectories of the diffusing species.« less

An MRI study of white matter tract integrity in regular cannabis users: effects of cannabis use and age.

PubMed

Jakabek, David; Yücel, Murat; Lorenzetti, Valentina; Solowij, Nadia

2016-10-01

Conflicting evidence exists on the effects of cannabis use on brain white matter integrity. The extant literature has exclusively focused on younger cannabis users, with no studies sampling older cannabis users. We recruited a sample with a broad age range to examine the integrity of major white matter tracts in association with cannabis use parameters and neurodevelopmental stage. Regular cannabis users (n = 56) and non-users (n = 20) with a mean age of 32 (range 18-55 years) underwent structural and diffusion MRI scans. White matter was examined using voxel-based statistics and via probabilistic tract reconstruction. The integrity of tracts was assessed using average fractional anisotropy, axial diffusivity and radial diffusivity. Diffusion measures were compared between users and non-users and as group-by-age interactions. Correlations between diffusion measures and age of onset, duration, frequency and dose of current cannabis use were examined. Cannabis users overall had lower fractional anisotropy than healthy non-users in the forceps minor tract only (p = .015, partial eta = 0.07), with no voxel-wise differences observed. Younger users showed predominantly reduced axial diffusivity, whereas older users had higher radial diffusivity in widespread tracts. Higher axial diffusivity was associated with duration of cannabis use in the cingulum angular bundle (beta = 5.00 × 10(-5), p = .003). Isolated higher AD in older cannabis users was also observed. The findings suggest that exogenous cannabinoids alter normal brain maturation, with differing effects at various neurodevelopmental stages of life. These age-related differences are posited to account for the disparate results described in the literature.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sarman, Sten, E-mail: sarman@ownit.nu; Wang, Yong-Lei; Laaksonen, Aatto

The self-diffusion coefficients of nematic phases of various model systems consisting of regular convex calamitic and discotic ellipsoids and non-convex bodies such as bent-core molecules and soft ellipsoid strings have been obtained as functions of the shear rate in a shear flow. Then the self-diffusion coefficient is a second rank tensor with three different diagonal components and two off-diagonal components. These coefficients were found to be determined by a combination of two mechanisms, which previously have been found to govern the self-diffusion of shearing isotropic liquids, namely, (i) shear alignment enhancing the diffusion in the direction parallel to the streamlinesmore » and hindering the diffusion in the perpendicular directions and (ii) the distortion of the shell structure in the liquid whereby a molecule more readily can escape from a surrounding shell of nearest neighbors, so that the mobility increases in every direction. Thus, the diffusion parallel to the streamlines always increases with the shear rate since these mechanisms cooperate in this direction. In the perpendicular directions, these mechanisms counteract each other so that the behaviour becomes less regular. In the case of the nematic phases of the calamitic and discotic ellipsoids and of the bent core molecules, mechanism (ii) prevails so that the diffusion coefficients increase. However, the diffusion coefficients of the soft ellipsoid strings decrease in the direction of the velocity gradient because the broadsides of these molecules are oriented perpendicularly to this direction due the shear alignment (i). The cross coupling coefficient relating a gradient of tracer particles in the direction of the velocity gradient and their flow in the direction of the streamlines is negative and rather large, whereas the other coupling coefficient relating a gradient in the direction of the streamlines and a flow in the direction of the velocity gradient is very small.« less

Data mining of molecular dynamics data reveals Li diffusion characteristics in garnet Li7La3Zr2O12

PubMed Central

Chen, Chi; Lu, Ziheng; Ciucci, Francesco

2017-01-01

Understanding Li diffusion in solid conductors is essential for the next generation Li batteries. Here we show that density-based clustering of the trajectories computed using molecular dynamics simulations helps elucidate the Li diffusion mechanism within the Li7La3Zr2O12 (LLZO) crystal lattice. This unsupervised learning method recognizes lattice sites, is able to give the site type, and can identify Li hopping events. Results show that, while the cubic LLZO has a much higher hopping rate compared to its tetragonal counterpart, most of the Li hops in the cubic LLZO do not contribute to the diffusivity due to the dominance of back-and-forth type jumps. The hopping analysis and local Li configuration statistics give evidence that Li diffusivity in cubic LLZO is limited by the low vacancy concentration. The hopping statistics also shows uncorrelated Poisson-like diffusion for Li in the cubic LLZO, and correlated diffusion for Li in the tetragonal LLZO in the temporal scale. Further analysis of the spatio-temporal correlation using site-to-site mutual information confirms the weak site dependence of Li diffusion in the cubic LLZO as the origin for the uncorrelated diffusion. This work puts forward a perspective on combining machine learning and information theory to interpret results of molecular dynamics simulations. PMID:28094317

Data mining of molecular dynamics data reveals Li diffusion characteristics in garnet Li7La3Zr2O12

NASA Astrophysics Data System (ADS)

Chen, Chi; Lu, Ziheng; Ciucci, Francesco

2017-01-01

Understanding Li diffusion in solid conductors is essential for the next generation Li batteries. Here we show that density-based clustering of the trajectories computed using molecular dynamics simulations helps elucidate the Li diffusion mechanism within the Li7La3Zr2O12 (LLZO) crystal lattice. This unsupervised learning method recognizes lattice sites, is able to give the site type, and can identify Li hopping events. Results show that, while the cubic LLZO has a much higher hopping rate compared to its tetragonal counterpart, most of the Li hops in the cubic LLZO do not contribute to the diffusivity due to the dominance of back-and-forth type jumps. The hopping analysis and local Li configuration statistics give evidence that Li diffusivity in cubic LLZO is limited by the low vacancy concentration. The hopping statistics also shows uncorrelated Poisson-like diffusion for Li in the cubic LLZO, and correlated diffusion for Li in the tetragonal LLZO in the temporal scale. Further analysis of the spatio-temporal correlation using site-to-site mutual information confirms the weak site dependence of Li diffusion in the cubic LLZO as the origin for the uncorrelated diffusion. This work puts forward a perspective on combining machine learning and information theory to interpret results of molecular dynamics simulations.

Social Media Use and Access to Digital Technology in US Young Adults in 2016.

PubMed

Villanti, Andrea C; Johnson, Amanda L; Ilakkuvan, Vinu; Jacobs, Megan A; Graham, Amanda L; Rath, Jessica M

2017-06-07

In 2015, 90% of US young adults with Internet access used social media. Digital and social media are highly prevalent modalities through which young adults explore identity formation, and by extension, learn and transmit norms about health and risk behaviors during this developmental life stage. The purpose of this study was to provide updated estimates of social media use from 2014 to 2016 and correlates of social media use and access to digital technology in data collected from a national sample of US young adults in 2016. Young adult participants aged 18-24 years in Wave 7 (October 2014, N=1259) and Wave 9 (February 2016, N=989) of the Truth Initiative Young Adult Cohort Study were asked about use frequency for 11 social media sites and access to digital devices, in addition to sociodemographic characteristics. Regular use was defined as using a given social media site at least weekly. Weighted analyses estimated the prevalence of use of each social media site, overlap between regular use of specific sites, and correlates of using a greater number of social media sites regularly. Bivariate analyses identified sociodemographic correlates of access to specific digital devices. In 2014, 89.42% (weighted n, 1126/1298) of young adults reported regular use of at least one social media site. This increased to 97.5% (weighted n, 965/989) of young adults in 2016. Among regular users of social media sites in 2016, the top five sites were Tumblr (85.5%), Vine (84.7%), Snapchat (81.7%), Instagram (80.7%), and LinkedIn (78.9%). Respondents reported regularly using an average of 7.6 social media sites, with 85% using 6 or more sites regularly. Overall, 87% of young adults reported access or use of a smartphone with Internet access, 74% a desktop or laptop computer with Internet access, 41% a tablet with Internet access, 29% a smart TV or video game console with Internet access, 11% a cell phone without Internet access, and 3% none of these. Access to all digital devices with Internet was lower in those reporting a lower subjective financial situation; there were also significant differences in access to specific digital devices with Internet by race, ethnicity, and education. The high mean number of social media sites used regularly and the substantial overlap in use of multiple social media sites reflect the rapidly changing social media environment. Mobile devices are a primary channel for social media, and our study highlights disparities in access to digital technologies with Internet access among US young adults by race/ethnicity, education, and subjective financial status. Findings from this study may guide the development and implementation of future health interventions for young adults delivered via the Internet or social media sites. ©Andrea C Villanti, Amanda L Johnson, Vinu Ilakkuvan, Megan A Jacobs, Amanda L Graham, Jessica M Rath. Originally published in the Journal of Medical Internet Research (http://www.jmir.org), 07.06.2017.

Social Media Use and Access to Digital Technology in US Young Adults in 2016

PubMed Central

Johnson, Amanda L; Ilakkuvan, Vinu; Jacobs, Megan A; Graham, Amanda L; Rath, Jessica M

2017-01-01

Background In 2015, 90% of US young adults with Internet access used social media. Digital and social media are highly prevalent modalities through which young adults explore identity formation, and by extension, learn and transmit norms about health and risk behaviors during this developmental life stage. Objective The purpose of this study was to provide updated estimates of social media use from 2014 to 2016 and correlates of social media use and access to digital technology in data collected from a national sample of US young adults in 2016. Methods Young adult participants aged 18-24 years in Wave 7 (October 2014, N=1259) and Wave 9 (February 2016, N=989) of the Truth Initiative Young Adult Cohort Study were asked about use frequency for 11 social media sites and access to digital devices, in addition to sociodemographic characteristics. Regular use was defined as using a given social media site at least weekly. Weighted analyses estimated the prevalence of use of each social media site, overlap between regular use of specific sites, and correlates of using a greater number of social media sites regularly. Bivariate analyses identified sociodemographic correlates of access to specific digital devices. Results In 2014, 89.42% (weighted n, 1126/1298) of young adults reported regular use of at least one social media site. This increased to 97.5% (weighted n, 965/989) of young adults in 2016. Among regular users of social media sites in 2016, the top five sites were Tumblr (85.5%), Vine (84.7%), Snapchat (81.7%), Instagram (80.7%), and LinkedIn (78.9%). Respondents reported regularly using an average of 7.6 social media sites, with 85% using 6 or more sites regularly. Overall, 87% of young adults reported access or use of a smartphone with Internet access, 74% a desktop or laptop computer with Internet access, 41% a tablet with Internet access, 29% a smart TV or video game console with Internet access, 11% a cell phone without Internet access, and 3% none of these. Access to all digital devices with Internet was lower in those reporting a lower subjective financial situation; there were also significant differences in access to specific digital devices with Internet by race, ethnicity, and education. Conclusions The high mean number of social media sites used regularly and the substantial overlap in use of multiple social media sites reflect the rapidly changing social media environment. Mobile devices are a primary channel for social media, and our study highlights disparities in access to digital technologies with Internet access among US young adults by race/ethnicity, education, and subjective financial status. Findings from this study may guide the development and implementation of future health interventions for young adults delivered via the Internet or social media sites. PMID:28592394

Inverse Diffusion Curves Using Shape Optimization.

PubMed

Zhao, Shuang; Durand, Fredo; Zheng, Changxi

2018-07-01

The inverse diffusion curve problem focuses on automatic creation of diffusion curve images that resemble user provided color fields. This problem is challenging since the 1D curves have a nonlinear and global impact on resulting color fields via a partial differential equation (PDE). We introduce a new approach complementary to previous methods by optimizing curve geometry. In particular, we propose a novel iterative algorithm based on the theory of shape derivatives. The resulting diffusion curves are clean and well-shaped, and the final image closely approximates the input. Our method provides a user-controlled parameter to regularize curve complexity, and generalizes to handle input color fields represented in a variety of formats.

Electrochemical Impedance Imaging via the Distribution of Diffusion Times.

PubMed

Song, Juhyun; Bazant, Martin Z

2018-03-16

We develop a mathematical framework to analyze electrochemical impedance spectra in terms of a distribution of diffusion times (DDT) for a parallel array of random finite-length Warburg (diffusion) or Gerischer (reaction-diffusion) circuit elements. A robust DDT inversion method is presented based on complex nonlinear least squares regression with Tikhonov regularization and illustrated for three cases of nanostructured electrodes for energy conversion: (i) a carbon nanotube supercapacitor, (ii) a silicon nanowire Li-ion battery, and (iii) a porous-carbon vanadium flow battery. The results demonstrate the feasibility of nondestructive "impedance imaging" to infer microstructural statistics of random, heterogeneous materials.

Novel Harmonic Regularization Approach for Variable Selection in Cox's Proportional Hazards Model

PubMed Central

Chu, Ge-Jin; Liang, Yong; Wang, Jia-Xuan

2014-01-01

Variable selection is an important issue in regression and a number of variable selection methods have been proposed involving nonconvex penalty functions. In this paper, we investigate a novel harmonic regularization method, which can approximate nonconvex Lq  (1/2 < q < 1) regularizations, to select key risk factors in the Cox's proportional hazards model using microarray gene expression data. The harmonic regularization method can be efficiently solved using our proposed direct path seeking approach, which can produce solutions that closely approximate those for the convex loss function and the nonconvex regularization. Simulation results based on the artificial datasets and four real microarray gene expression datasets, such as real diffuse large B-cell lymphoma (DCBCL), the lung cancer, and the AML datasets, show that the harmonic regularization method can be more accurate for variable selection than existing Lasso series methods. PMID:25506389

Thermogravimetric analysis for rapid assessment of moisture diffusivity in polydisperse powder and thin film matrices.

PubMed

Thirunathan, Praveena; Arnz, Patrik; Husny, Joeska; Gianfrancesco, Alessandro; Perdana, Jimmy

2018-03-01

Accurate description of moisture diffusivity is key to precisely understand and predict moisture transfer behaviour in a matrix. Unfortunately, measuring moisture diffusivity is not trivial, especially at low moisture values and/or elevated temperatures. This paper presents a novel experimental procedure to accurately measure moisture diffusivity based on thermogravimetric approach. The procedure is capable to measure diffusivity even at elevated temperatures (>70°C) and low moisture values (>1%). Diffusivity was extracted from experimental data based on "regular regime approach". The approach was tailored to determine diffusivity from thin film and from poly-dispersed powdered samples. Subsequently, measured diffusivity was validated by comparing to available literature data, showing good agreement. Ability of this approach to accurately measure diffusivity at a wider range of temperatures provides better insight on temperature dependency of diffusivity. Thus, this approach can be crucial to ensure good accuracy of moisture transfer description/prediction especially when involving elevated temperatures. Copyright © 2017 Elsevier Ltd. All rights reserved.

Nonequilibrium diffusive gas dynamics: Poiseuille microflow

NASA Astrophysics Data System (ADS)

Abramov, Rafail V.; Otto, Jasmine T.

2018-05-01

We test the recently developed hierarchy of diffusive moment closures for gas dynamics together with the near-wall viscosity scaling on the Poiseuille flow of argon and nitrogen in a one micrometer wide channel, and compare it against the corresponding Direct Simulation Monte Carlo computations. We find that the diffusive regularized Grad equations with viscosity scaling provide the most accurate approximation to the benchmark DSMC results. At the same time, the conventional Navier-Stokes equations without the near-wall viscosity scaling are found to be the least accurate among the tested closures.

From localization to anomalous diffusion in the dynamics of coupled kicked rotors

NASA Astrophysics Data System (ADS)

Notarnicola, Simone; Iemini, Fernando; Rossini, Davide; Fazio, Rosario; Silva, Alessandro; Russomanno, Angelo

2018-02-01

We study the effect of many-body quantum interference on the dynamics of coupled periodically kicked systems whose classical dynamics is chaotic and shows an unbounded energy increase. We specifically focus on an N -coupled kicked rotors model: We find that the interplay of quantumness and interactions dramatically modifies the system dynamics, inducing a transition between energy saturation and unbounded energy increase. We discuss this phenomenon both numerically and analytically through a mapping onto an N -dimensional Anderson model. The thermodynamic limit N →∞ , in particular, always shows unbounded energy growth. This dynamical delocalization is genuinely quantum and very different from the classical one: Using a mean-field approximation, we see that the system self-organizes so that the energy per site increases in time as a power law with exponent smaller than 1. This wealth of phenomena is a genuine effect of quantum interference: The classical system for N ≥2 always behaves ergodically with an energy per site linearly increasing in time. Our results show that quantum mechanics can deeply alter the regularity or ergodicity properties of a many-body-driven system.

Development and application of an exchange model for anisotropic water diffusion in the microporous MOF aluminum fumarate

NASA Astrophysics Data System (ADS)

Splith, Tobias; Fröhlich, Dominik; Henninger, Stefan K.; Stallmach, Frank

2018-06-01

Diffusion of water in aluminum fumarate was studied by means of pulsed field gradient (PFG) nuclear magnetic resonance (NMR). Due to water molecules exchanging between the intracrystalline anisotropic pore space and the isotropic intercrystalline void space the model of intracrystalline anisotropic diffusion fails to describe the experimental PFG NMR data at high observation times. Therefore, the two-site exchange model developed by Kärger is extended to the case of exchange between an anisotropic and an isotropic site. This extended exchange model is solved by numerical integration. It describes the experimental data very well and yields values for the intracrystalline diffusion coefficient and the mean residence times of the respective sites. Further PFG NMR studies were performed with coatings consisting of small aluminum fumarate crystals, which are used in adsorptive heat transformation applications. The diffusion coefficients of water in the small crystal coating are compared to the values expected from the extended two-site exchange model and from the model of long-range diffusion.

Obstructed metabolite diffusion within skeletal muscle cells in silico.

PubMed

Aliev, Mayis K; Tikhonov, Alexander N

2011-12-01

Using a Monte Carlo simulation technique, we have modeled 3D diffusion of low molecular weight metabolites inside a skeletal muscle cell. The following structural elements are considered: (i) a regular lattice of actin and myosin filaments inside a myofibril, (ii) the membranes of sarcoplasmic reticulum and mitochondria surrounding the myofibrils, (iii) a set of myofibrils inside a skeletal muscle cell encircled by the outer cell membrane, and (iv) an additional set of regular intracellular structures ("macrocompartments") embedded into the cell interior. The macrocompartments are considered to simulate diffusion restrictions because of hypothetical cylindrical structures (16-22 μm in diameter) suggested earlier (de Graaf et al. Biophys J 78: 1657-1664, 2000). This model allowed us to calculate the apparent coefficients of particle diffusion in the radial and axial directions, D(app)(⊥) and D(app)(II), respectively. Particle movements in the axial direction are considered, at first approximation, as unrestricted diffusion (D(app)(II) = const). The apparent coefficient of radial diffusion, D(app)(⊥), decreases with time because of particle collisions with myofilaments and other rigid obstacles. Results of our random walk simulations are in fairly good agreement with experimental data on NMR measurements of restricted radial diffusion of phosphocreatine in white and red skeletal muscles of goldfish (Kinsey et al. NMR Biomed 12:1-7, 1999). Particle reflections from the low-permeable borders of macrocompartments (efficient diameter, D(eff)(MC) ≈ 9.2-10.4 μm) are the prerequisite for agreeing theoretical and experimental data. The low-permeable coverage of hypothetical macrocompartments (99.8% of coverage) provides the main contribution to time-dependent decrease in D(app)(⊥).

Self-diffusion in the non-Newtonian regime of shearing liquid crystal model systems based on the Gay-Berne potential

NASA Astrophysics Data System (ADS)

Sarman, Sten; Wang, Yong-Lei; Laaksonen, Aatto

2016-02-01

The self-diffusion coefficients of nematic phases of various model systems consisting of regular convex calamitic and discotic ellipsoids and non-convex bodies such as bent-core molecules and soft ellipsoid strings have been obtained as functions of the shear rate in a shear flow. Then the self-diffusion coefficient is a second rank tensor with three different diagonal components and two off-diagonal components. These coefficients were found to be determined by a combination of two mechanisms, which previously have been found to govern the self-diffusion of shearing isotropic liquids, namely, (i) shear alignment enhancing the diffusion in the direction parallel to the streamlines and hindering the diffusion in the perpendicular directions and (ii) the distortion of the shell structure in the liquid whereby a molecule more readily can escape from a surrounding shell of nearest neighbors, so that the mobility increases in every direction. Thus, the diffusion parallel to the streamlines always increases with the shear rate since these mechanisms cooperate in this direction. In the perpendicular directions, these mechanisms counteract each other so that the behaviour becomes less regular. In the case of the nematic phases of the calamitic and discotic ellipsoids and of the bent core molecules, mechanism (ii) prevails so that the diffusion coefficients increase. However, the diffusion coefficients of the soft ellipsoid strings decrease in the direction of the velocity gradient because the broadsides of these molecules are oriented perpendicularly to this direction due the shear alignment (i). The cross coupling coefficient relating a gradient of tracer particles in the direction of the velocity gradient and their flow in the direction of the streamlines is negative and rather large, whereas the other coupling coefficient relating a gradient in the direction of the streamlines and a flow in the direction of the velocity gradient is very small.

An efficient method for model refinement in diffuse optical tomography

NASA Astrophysics Data System (ADS)

Zirak, A. R.; Khademi, M.

2007-11-01

Diffuse optical tomography (DOT) is a non-linear, ill-posed, boundary value and optimization problem which necessitates regularization. Also, Bayesian methods are suitable owing to measurements data are sparse and correlated. In such problems which are solved with iterative methods, for stabilization and better convergence, the solution space must be small. These constraints subject to extensive and overdetermined system of equations which model retrieving criteria specially total least squares (TLS) must to refine model error. Using TLS is limited to linear systems which is not achievable when applying traditional Bayesian methods. This paper presents an efficient method for model refinement using regularized total least squares (RTLS) for treating on linearized DOT problem, having maximum a posteriori (MAP) estimator and Tikhonov regulator. This is done with combination Bayesian and regularization tools as preconditioner matrices, applying them to equations and then using RTLS to the resulting linear equations. The preconditioning matrixes are guided by patient specific information as well as a priori knowledge gained from the training set. Simulation results illustrate that proposed method improves the image reconstruction performance and localize the abnormally well.

Multi-site Study of Diffusion Metric Variability: Characterizing the Effects of Site, Vendor, Field Strength, and Echo Time using the Histogram Distance.

PubMed

Helmer, K G; Chou, M-C; Preciado, R I; Gimi, B; Rollins, N K; Song, A; Turner, J; Mori, S

2016-02-27

MRI-based multi-site trials now routinely include some form of diffusion-weighted imaging (DWI) in their protocol. These studies can include data originating from scanners built by different vendors, each with their own set of unique protocol restrictions, including restrictions on the number of available gradient directions, whether an externally-generated list of gradient directions can be used, and restrictions on the echo time (TE). One challenge of multi-site studies is to create a common imaging protocol that will result in a reliable and accurate set of diffusion metrics. The present study describes the effect of site, scanner vendor, field strength, and TE on two common metrics: the first moment of the diffusion tensor field (mean diffusivity, MD), and the fractional anisotropy (FA). We have shown in earlier work that ROI metrics and the mean of MD and FA histograms are not sufficiently sensitive for use in site characterization. Here we use the distance between whole brain histograms of FA and MD to investigate within- and between-site effects. We concluded that the variability of DTI metrics due to site, vendor, field strength, and echo time could influence the results in multi-center trials and that histogram distance is sensitive metrics for each of these variables.

Insights into new-onset atrial fibrillation following open heart surgery and implications for type II atrial flutter.

PubMed

Sadrpour, Shervin A; Srinivasan, Deepa; Bhimani, Ashish A; Lee, Seungyup; Ryu, Kyungmoo; Cakulev, Ivan; Khrestian, Celeen M; Markowitz, Alan H; Waldo, Albert L; Sahadevan, Jayakumar

2015-12-01

Postoperative atrial fibrillation (POAF), new-onset AF after open heart surgery (OHS), is thought to be related to pericarditis. Based on AF studies in the canine sterile pericarditis model, we hypothesized that POAF in patients after OHS may be associated with a rapid, regular rhythm in the left atrium (LA), suggestive of an LA driver maintaining AF. The aim of this study was to test the hypothesis that in patients with POAF, atrial electrograms (AEGs) recorded from at least one of the two carefully selected LA sites would manifest a rapid, regular rhythm with AEGs of short cycle length (CL) and constant morphology, but a selected right atrial (RA) site would manifest AEGs with irregular CLs and variable morphology. In 44 patients undergoing OHS, AEGs recorded from the epicardial surface of the RA, the LA portion of Bachmann's bundle, and the posterior LA during sustained AF were analysed for regularity of CL and morphology. Sustained AF occurred in 15 of 44 patients. Atrial electrograms were recorded in 11 of 15 patients; 8 of 11 had rapid, regular activation with constant morphology recorded from at least one LA site; no regular AEG sites were present in 3 of 11 patients. Atrial electrograms recorded during sustained POAF frequently demonstrated rapid, regular activation in at least one LA site, consistent with a driver maintaining AF. Published by Oxford University Press on behalf of the European Society of Cardiology 2015. This work is written by (a) US Government employee(s) and is in the public domain in the US.

Threshold-based epidemic dynamics in systems with memory

NASA Astrophysics Data System (ADS)

Bodych, Marcin; Ganguly, Niloy; Krueger, Tyll; Mukherjee, Animesh; Siegmund-Schultze, Rainer; Sikdar, Sandipan

2016-11-01

In this article we analyze an epidemic dynamics model (SI) where we assume that there are k susceptible states, that is a node would require multiple (k) contacts before it gets infected. In specific, we provide a theoretical framework for studying diffusion rate in complete graphs and d-regular trees with extensions to dense random graphs. We observe that irrespective of the topology, the diffusion process could be divided into two distinct phases: i) the initial phase, where the diffusion process is slow, followed by ii) the residual phase where the diffusion rate increases manifold. In fact, the initial phase acts as an indicator for the total diffusion time in dense graphs. The most remarkable lesson from this investigation is that such a diffusion process could be controlled and even contained if acted upon within its initial phase.

Atomistic Modeling of Cation Diffusion in Transition Metal Perovskites La1-xSrxMnO3+/-δfor Solid Oxide Fuel Cell Cathodes Applications

NASA Astrophysics Data System (ADS)

Lee, Yueh-Lin; Duan, Yuhua; Morgan, Dane; Sorescu, Dan; Abernathy, Harry

Cation diffusion in La1-xSrxMnO3+/-δ (LSM) and in related perovskite materials play an important role in controlling long term performance and stability of solid oxide fuel cell (SOFCs) cathodes. Due to sluggish rates of cation diffusion and complex coupling between defect chemistry and cation diffusion pathways, currently there is still lack of quantitative theoretical model predictions on cation diffusivity vs. T and P(O2) to describe experimental cation tracer diffusivities. In this work, based on ab initio modeling of LSM defect chemistry and migration barriers of the possible cation diffusion pathways, we assess the rates of A-site and B-site cation diffusion in a wide range of T and P(O2) at x =0.0 and 0.2 for SOFC applications. We demonstrate the active cation diffusion pathways in LSM involve cation defect clusters as cation transport carriers, where reduction in the cation migration barriers, which are governed by the steric effect associated with the metal-oxygen cage in the perovskite lattice, is much greater than the penalty of repulsive interaction in the A-site and B-site cation vacancy clusters, leading to higher cation diffusion rates as compared to those of single cation vacancy hopping mechanisms. The predicted Mn and La/Sr cation self-diffusion coefficients of LSM at at x =0.0 and 0.2 along with their 1/T and P(O2) dependences, are in good agreement with the experimental tracer diffusion coefficients.

Multisite longitudinal reliability of tract-based spatial statistics in diffusion tensor imaging of healthy elderly subjects.

PubMed

Jovicich, Jorge; Marizzoni, Moira; Bosch, Beatriz; Bartrés-Faz, David; Arnold, Jennifer; Benninghoff, Jens; Wiltfang, Jens; Roccatagliata, Luca; Picco, Agnese; Nobili, Flavio; Blin, Oliver; Bombois, Stephanie; Lopes, Renaud; Bordet, Régis; Chanoine, Valérie; Ranjeva, Jean-Philippe; Didic, Mira; Gros-Dagnac, Hélène; Payoux, Pierre; Zoccatelli, Giada; Alessandrini, Franco; Beltramello, Alberto; Bargalló, Núria; Ferretti, Antonio; Caulo, Massimo; Aiello, Marco; Ragucci, Monica; Soricelli, Andrea; Salvadori, Nicola; Tarducci, Roberto; Floridi, Piero; Tsolaki, Magda; Constantinidis, Manos; Drevelegas, Antonios; Rossini, Paolo Maria; Marra, Camillo; Otto, Josephin; Reiss-Zimmermann, Martin; Hoffmann, Karl-Titus; Galluzzi, Samantha; Frisoni, Giovanni B

2014-11-01

Large-scale longitudinal neuroimaging studies with diffusion imaging techniques are necessary to test and validate models of white matter neurophysiological processes that change in time, both in healthy and diseased brains. The predictive power of such longitudinal models will always be limited by the reproducibility of repeated measures acquired during different sessions. At present, there is limited quantitative knowledge about the across-session reproducibility of standard diffusion metrics in 3T multi-centric studies on subjects in stable conditions, in particular when using tract based spatial statistics and with elderly people. In this study we implemented a multi-site brain diffusion protocol in 10 clinical 3T MRI sites distributed across 4 countries in Europe (Italy, Germany, France and Greece) using vendor provided sequences from Siemens (Allegra, Trio Tim, Verio, Skyra, Biograph mMR), Philips (Achieva) and GE (HDxt) scanners. We acquired DTI data (2 × 2 × 2 mm(3), b = 700 s/mm(2), 5 b0 and 30 diffusion weighted volumes) of a group of healthy stable elderly subjects (5 subjects per site) in two separate sessions at least a week apart. For each subject and session four scalar diffusion metrics were considered: fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD) and axial (AD) diffusivity. The diffusion metrics from multiple subjects and sessions at each site were aligned to their common white matter skeleton using tract-based spatial statistics. The reproducibility at each MRI site was examined by looking at group averages of absolute changes relative to the mean (%) on various parameters: i) reproducibility of the signal-to-noise ratio (SNR) of the b0 images in centrum semiovale, ii) full brain test-retest differences of the diffusion metric maps on the white matter skeleton, iii) reproducibility of the diffusion metrics on atlas-based white matter ROIs on the white matter skeleton. Despite the differences of MRI scanner configurations across sites (vendors, models, RF coils and acquisition sequences) we found good and consistent test-retest reproducibility. White matter b0 SNR reproducibility was on average 7 ± 1% with no significant MRI site effects. Whole brain analysis resulted in no significant test-retest differences at any of the sites with any of the DTI metrics. The atlas-based ROI analysis showed that the mean reproducibility errors largely remained in the 2-4% range for FA and AD and 2-6% for MD and RD, averaged across ROIs. Our results show reproducibility values comparable to those reported in studies using a smaller number of MRI scanners, slightly different DTI protocols and mostly younger populations. We therefore show that the acquisition and analysis protocols used are appropriate for multi-site experimental scenarios. Copyright © 2014 Elsevier Inc. All rights reserved.

Modification of the flow pass method as applied to problems of chemistry of planet atmospheres

NASA Technical Reports Server (NTRS)

Parshev, V. A.

1980-01-01

It was shown that the modified flow pass method possesses considerable effectiveness, both in the case when the coefficient of diffusion changes severely in the examined region and in the case when diffusion is the prevalent process, as compared with chemical reactions. The case when a regular pass proves inapplicable, or applicable in a limited interval of the decisive parameters, was examined.

Solving regularly and singularly perturbed reaction-diffusion equations in three space dimensions

NASA Astrophysics Data System (ADS)

Moore, Peter K.

2007-06-01

In [P.K. Moore, Effects of basis selection and h-refinement on error estimator reliability and solution efficiency for higher-order methods in three space dimensions, Int. J. Numer. Anal. Mod. 3 (2006) 21-51] a fixed, high-order h-refinement finite element algorithm, Href, was introduced for solving reaction-diffusion equations in three space dimensions. In this paper Href is coupled with continuation creating an automatic method for solving regularly and singularly perturbed reaction-diffusion equations. The simple quasilinear Newton solver of Moore, (2006) is replaced by the nonlinear solver NITSOL [M. Pernice, H.F. Walker, NITSOL: a Newton iterative solver for nonlinear systems, SIAM J. Sci. Comput. 19 (1998) 302-318]. Good initial guesses for the nonlinear solver are obtained using continuation in the small parameter ɛ. Two strategies allow adaptive selection of ɛ. The first depends on the rate of convergence of the nonlinear solver and the second implements backtracking in ɛ. Finally a simple method is used to select the initial ɛ. Several examples illustrate the effectiveness of the algorithm.

MODIS Solar Diffuser: Modelled and Actual Performance

NASA Technical Reports Server (NTRS)

Waluschka, Eugene; Xiong, Xiao-Xiong; Esposito, Joe; Wang, Xin-Dong; Krebs, Carolyn (Technical Monitor)

2001-01-01

The Moderate Resolution Imaging Spectroradiometer (MODIS) instrument's solar diffuser is used in its radiometric calibration for the reflective solar bands (VIS, NTR, and SWIR) ranging from 0.41 to 2.1 micron. The sun illuminates the solar diffuser either directly or through a attenuation screen. The attenuation screen consists of a regular array of pin holes. The attenuated illumination pattern on the solar diffuser is not uniform, but consists of a multitude of pin-hole images of the sun. This non-uniform illumination produces small, but noticeable radiometric effects. A description of the computer model used to simulate the effects of the attenuation screen is given and the predictions of the model are compared with actual, on-orbit, calibration measurements.

Evolution of the concentration PDF in random environments modeled by global random walk

NASA Astrophysics Data System (ADS)

Suciu, Nicolae; Vamos, Calin; Attinger, Sabine; Knabner, Peter

2013-04-01

The evolution of the probability density function (PDF) of concentrations of chemical species transported in random environments is often modeled by ensembles of notional particles. The particles move in physical space along stochastic-Lagrangian trajectories governed by Ito equations, with drift coefficients given by the local values of the resolved velocity field and diffusion coefficients obtained by stochastic or space-filtering upscaling procedures. A general model for the sub-grid mixing also can be formulated as a system of Ito equations solving for trajectories in the composition space. The PDF is finally estimated by the number of particles in space-concentration control volumes. In spite of their efficiency, Lagrangian approaches suffer from two severe limitations. Since the particle trajectories are constructed sequentially, the demanded computing resources increase linearly with the number of particles. Moreover, the need to gather particles at the center of computational cells to perform the mixing step and to estimate statistical parameters, as well as the interpolation of various terms to particle positions, inevitably produce numerical diffusion in either particle-mesh or grid-free particle methods. To overcome these limitations, we introduce a global random walk method to solve the system of Ito equations in physical and composition spaces, which models the evolution of the random concentration's PDF. The algorithm consists of a superposition on a regular lattice of many weak Euler schemes for the set of Ito equations. Since all particles starting from a site of the space-concentration lattice are spread in a single numerical procedure, one obtains PDF estimates at the lattice sites at computational costs comparable with those for solving the system of Ito equations associated to a single particle. The new method avoids the limitations concerning the number of particles in Lagrangian approaches, completely removes the numerical diffusion, and speeds up the computation by orders of magnitude. The approach is illustrated for the transport of passive scalars in heterogeneous aquifers, with hydraulic conductivity modeled as a random field.

Correlation factors for atomic diffusion in nondilute multicomponent alloys with arbitrary vacancy concentration

NASA Astrophysics Data System (ADS)

Tahir-Kheli, R. A.

1983-09-01

Vacancy-assisted tracer diffusion in a multicomponent kinetic alloy consisting of xλ N atoms with hopping rate Jλ (where λ≡A, B, C, etc.) and υN vacancies (where υ=1-λxλ) distributed randomly over a regular d-dimensional (where d>=2) hypercubic, or close-packed, lattice of N sites is analyzed through a self-consistent renormalization of a recent theory of Tahir-Kheli and Elliott combined with a generalization of concepts introduced by Manning. The result for the tracer-diffusion correlation factor is the following: ftr=H''(tr)[H''(tr)+2J0], where J0 is the tracer-hopping rate, H''(tr) is a generalized effective vacancy escape frequency, H''(tr)=[M(1-υ)[J0υftr+Jeff], where Jeff is an effective hopping rate of the background atoms averaged with a weighting factor proportional to xλ and fλ, i.e., Jeff=λ(Jλxλfλ)λ(xλfλ) and M=-(1+<θ>)<θ>. For a single-component alloy, with particle concentration x, Jλ=J, and vacancy concentration υ=1-x our theory provides an excellent overall description of the correlation factor as long as JJ0>~z-2. Indeed, even for J-->0, the calculated results agree with the Monte Carlo estimates, except in the immediate vicinity of the percolation threshold, υp, which is located self-consistently to an accuracy of the order 1z.

Multiple binding sites for transcriptional repressors can produce regular bursting and enhance noise suppression

NASA Astrophysics Data System (ADS)

Lengyel, Iván M.; Morelli, Luis G.

2017-04-01

Cells may control fluctuations in protein levels by means of negative autoregulation, where transcription factors bind DNA sites to repress their own production. Theoretical studies have assumed a single binding site for the repressor, while in most species it is found that multiple binding sites are arranged in clusters. We study a stochastic description of negative autoregulation with multiple binding sites for the repressor. We find that increasing the number of binding sites induces regular bursting of gene products. By tuning the threshold for repression, we show that multiple binding sites can also suppress fluctuations. Our results highlight possible roles for the presence of multiple binding sites of negative autoregulators.

Scale-free network provides an optimal pattern for knowledge transfer

NASA Astrophysics Data System (ADS)

Lin, Min; Li, Nan

2010-02-01

We study numerically the knowledge innovation and diffusion process on four representative network models, such as regular networks, small-world networks, random networks and scale-free networks. The average knowledge stock level as a function of time is measured and the corresponding growth diffusion time, τ is defined and computed. On the four types of networks, the growth diffusion times all depend linearly on the network size N as τ∼N, while the slope for scale-free network is minimal indicating the fastest growth and diffusion of knowledge. The calculated variance and spatial distribution of knowledge stock illustrate that optimal knowledge transfer performance is obtained on scale-free networks. We also investigate the transient pattern of knowledge diffusion on the four networks, and a qualitative explanation of this finding is proposed.

Sinc-Galerkin estimation of diffusivity in parabolic problems

NASA Technical Reports Server (NTRS)

Smith, Ralph C.; Bowers, Kenneth L.

1991-01-01

A fully Sinc-Galerkin method for the numerical recovery of spatially varying diffusion coefficients in linear partial differential equations is presented. Because the parameter recovery problems are inherently ill-posed, an output error criterion in conjunction with Tikhonov regularization is used to formulate them as infinite-dimensional minimization problems. The forward problems are discretized with a sinc basis in both the spatial and temporal domains thus yielding an approximate solution which displays an exponential convergence rate and is valid on the infinite time interval. The minimization problems are then solved via a quasi-Newton/trust region algorithm. The L-curve technique for determining an approximate value of the regularization parameter is briefly discussed, and numerical examples are given which show the applicability of the method both for problems with noise-free data as well as for those whose data contains white noise.

Fluid Registration of Diffusion Tensor Images Using Information Theory

PubMed Central

Chiang, Ming-Chang; Leow, Alex D.; Klunder, Andrea D.; Dutton, Rebecca A.; Barysheva, Marina; Rose, Stephen E.; McMahon, Katie L.; de Zubicaray, Greig I.; Toga, Arthur W.; Thompson, Paul M.

2008-01-01

We apply an information-theoretic cost metric, the symmetrized Kullback-Leibler (sKL) divergence, or J-divergence, to fluid registration of diffusion tensor images. The difference between diffusion tensors is quantified based on the sKL-divergence of their associated probability density functions (PDFs). Three-dimensional DTI data from 34 subjects were fluidly registered to an optimized target image. To allow large image deformations but preserve image topology, we regularized the flow with a large-deformation diffeomorphic mapping based on the kinematics of a Navier-Stokes fluid. A driving force was developed to minimize the J-divergence between the deforming source and target diffusion functions, while reorienting the flowing tensors to preserve fiber topography. In initial experiments, we showed that the sKL-divergence based on full diffusion PDFs is adaptable to higher-order diffusion models, such as high angular resolution diffusion imaging (HARDI). The sKL-divergence was sensitive to subtle differences between two diffusivity profiles, showing promise for nonlinear registration applications and multisubject statistical analysis of HARDI data. PMID:18390342

Assessment of using ultrasound images as prior for diffuse optical tomography regularization matrix

NASA Astrophysics Data System (ADS)

Althobaiti, Murad; Vavadi, Hamed; Zhu, Quing

2017-02-01

Imaging of tissue with Ultrasound-guided diffuse optical tomography (DOT) is a rising imaging technique to map hemoglobin concentrations within tissue for breast cancer detection and diagnosis. Near-infrared optical imaging received a lot of attention in research as a possible technique to be used for such purpose especially for breast tumors. Since DOT images contrast is closely related to oxygenation and deoxygenating of the hemoglobin, which is an important factor in differentiating malignant and benign tumors. One of the optical imaging modalities used is the diffused optical tomography (DOT); which probes deep scattering tissue (1-5cm) by NIR optical source-detector probe and detects NIR photons in the diffusive regime. The photons in the diffusive regime usually reach the detector without significant information about their source direction and the propagation path. Because of that, the optical reconstruction problem of the medium characteristics is ill-posed even with the tomography and Back-projection techniques. The accurate recovery of images requires an effective image reconstruction method. Here, we illustrate a method in which ultrasound images are encoded as prior for regularization of the inversion matrix. Results were evaluated using phantom experiments of low and high absorption contrasts. This method improves differentiation between the low and the high contrasts targets. Ultimately, this method could improve malignant and benign cases by increasing reconstructed absorption ratio of malignant to benign. Besides that, the phantom results show improvements in target shape as well as the spatial resolution of the DOT reconstructed images.

Measurement of radon diffusion in polyethylene based on alpha detection

NASA Astrophysics Data System (ADS)

Rau, Wolfgang

2012-02-01

Radon diffusion in different materials has been measured in the past. Usually the diffusion measurements are based on a direct determination of the amount of radon that diffuses through a thin layer of material. Here we present a method based on the measurement of the radon daughter products which are deposited inside the material. Looking at the decay of 210Po allows us to directly measure the exponential diffusion profile characterized by the diffusion length. In addition we can determine the solubility of radon in PE. We also describe a second method to determine the diffusion constant based on the short-lived radon daughter products 218Po and 214Po, using the identical experimental setup. Measurements for regular polyethylene (PE) and High Molecular Weight Polyethylene (HMWPE) yielded diffusion lengths of (1.3±0.3) mm and (0.8±0.2) mm and solubilities of 0.5±0.1 and 0.7±0.2, respectively, for the first method; the diffusion lengths extracted from the second method are noticeably larger which may be caused by different experimental conditions during diffusion.

Multi-site Study of Diffusion Metric Variability: Characterizing the Effects of Site, Vendor, Field Strength, and Echo Time using the Histogram Distance

PubMed Central

Helmer, K. G.; Chou, M-C.; Preciado, R. I.; Gimi, B.; Rollins, N. K.; Song, A.; Turner, J.; Mori, S.

2016-01-01

MRI-based multi-site trials now routinely include some form of diffusion-weighted imaging (DWI) in their protocol. These studies can include data originating from scanners built by different vendors, each with their own set of unique protocol restrictions, including restrictions on the number of available gradient directions, whether an externally-generated list of gradient directions can be used, and restrictions on the echo time (TE). One challenge of multi-site studies is to create a common imaging protocol that will result in a reliable and accurate set of diffusion metrics. The present study describes the effect of site, scanner vendor, field strength, and TE on two common metrics: the first moment of the diffusion tensor field (mean diffusivity, MD), and the fractional anisotropy (FA). We have shown in earlier work that ROI metrics and the mean of MD and FA histograms are not sufficiently sensitive for use in site characterization. Here we use the distance between whole brain histograms of FA and MD to investigate within- and between-site effects. We concluded that the variability of DTI metrics due to site, vendor, field strength, and echo time could influence the results in multi-center trials and that histogram distance is sensitive metrics for each of these variables. PMID:27350723

Mathematical Modeling the Geometric Regularity in Proteus Mirabilis Colonies

NASA Astrophysics Data System (ADS)

Zhang, Bin; Jiang, Yi; Minsu Kim Collaboration

Proteus Mirabilis colony exhibits striking spatiotemporal regularity, with concentric ring patterns with alternative high and low bacteria density in space, and periodicity for repetition process of growth and swarm in time. We present a simple mathematical model to explain the spatiotemporal regularity of P. Mirabilis colonies. We study a one-dimensional system. Using a reaction-diffusion model with thresholds in cell density and nutrient concentration, we recreated periodic growth and spread patterns, suggesting that the nutrient constraint and cell density regulation might be sufficient to explain the spatiotemporal periodicity in P. Mirabilis colonies. We further verify this result using a cell based model.

Solving ill-posed control problems by stabilized finite element methods: an alternative to Tikhonov regularization

NASA Astrophysics Data System (ADS)

Burman, Erik; Hansbo, Peter; Larson, Mats G.

2018-03-01

Tikhonov regularization is one of the most commonly used methods for the regularization of ill-posed problems. In the setting of finite element solutions of elliptic partial differential control problems, Tikhonov regularization amounts to adding suitably weighted least squares terms of the control variable, or derivatives thereof, to the Lagrangian determining the optimality system. In this note we show that the stabilization methods for discretely ill-posed problems developed in the setting of convection-dominated convection-diffusion problems, can be highly suitable for stabilizing optimal control problems, and that Tikhonov regularization will lead to less accurate discrete solutions. We consider some inverse problems for Poisson’s equation as an illustration and derive new error estimates both for the reconstruction of the solution from the measured data and reconstruction of the source term from the measured data. These estimates include both the effect of the discretization error and error in the measurements.

Product of the SNPP VIIRS SD Screen Transmittance and the SD BRDF (RSB) From Both Yaw Maneuver and Regular On-Orbit Data

NASA Technical Reports Server (NTRS)

Lei, Ning; Xiong, Xiaoxiong

2016-01-01

To assure data quality, the Earth-observing Visible Infrared Imaging Radiometer Suite (VIIRS) regularly performs on-orbit radiometric calibrations of its 22 spectral bands. The primary calibration radiance source for the reflective solar bands (RSBs) is a sunlit solar diffuser (SD). During the calibration process, sunlight goes through a perforated plate (the SD screen) and then strikes the SD. The SD scattered sunlight is used for the calibration, with the spectral radiance proportional to the product of the SD screen transmittance and the SD bidirectional reflectance distribution function (BRDF). The BRDF is decomposed to the product of its value at launch and a numerical factor quantifying its change since launch. Therefore, the RSB calibration requires accurate knowledge of the product of the SD screen transmittance and the BRDF (RSB; launch time). Previously, we calculated the product with yaw maneuver data and found that the product had improved accuracy over the prelaunch one. With both yaw maneuver and regular on orbit data, we were able to improve the accuracy of the SDSM screen transmittance and the product for the solar diffuser stability monitor SD view. In this study, we use both yaw maneuver and a small portion of regular on-orbit data to determine the product for the RSB SD view.

Extending the Diffuse Layer Model of Surface Acidity Behavior: III. Estimating Bound Site Activity Coefficients

EPA Science Inventory

Although detailed thermodynamic analyses of the 2-pK diffuse layer surface complexation model generally specify bound site activity coefficients for the purpose of accounting for those non-ideal excess free energies contributing to bound site electrochemical potentials, in applic...

Diffusion coefficient and shear viscosity of rigid water models.

PubMed

Tazi, Sami; Boţan, Alexandru; Salanne, Mathieu; Marry, Virginie; Turq, Pierre; Rotenberg, Benjamin

2012-07-18

We report the diffusion coefficient and viscosity of popular rigid water models: two non-polarizable ones (SPC/E with three sites, and TIP4P/2005 with four sites) and a polarizable one (Dang-Chang, four sites). We exploit the dependence of the diffusion coefficient on the system size (Yeh and Hummer 2004 J. Phys. Chem. B 108 15873) to obtain the size-independent value. This also provides an estimate of the viscosity of all water models, which we compare to the Green-Kubo result. In all cases, a good agreement is found. The TIP4P/2005 model is in better agreement with the experimental data for both diffusion and viscosity. The SPC/E and Dang-Chang models overestimate the diffusion coefficient and underestimate the viscosity.

Bone marrow metastasis of malignant melanoma in childhood arising within a congenital melanocytic nevus.

PubMed

Volejnikova, Jana; Bajciova, Viera; Sulovska, Lucie; Geierova, Marie; Buriankova, Eva; Jarosova, Marie; Hajduch, Marian; Sterba, Jaroslav; Mihal, Vladimir

2016-09-01

Malignant melanoma in childhood is infrequent and can arise within congenital melanocytic nevi. Spread of malignant melanoma to the bone marrow, especially in children, is extremely rare. Reported is a case of a 5-year-old boy with a congenital large melanocytic nevus of the head and neck who presented with a short history of low back and leg pain, fever and cervical lymphadenopathy. Despite regular follow-up by a dermatologist and plastic surgeon and repeatedly negative histology of previous partial excisions, diffuse bone marrow infiltration with malignant melanoma was diagnosed. The primary site was identified in the post-excision area. The disease progressed rapidly on ipilimumab immunotherapy and led to death at four months from the diagnosis. Surveillance is indispensable in children with a predisposition to melanoma and nonspecific symptoms such as bone pain, gait impairment or cytopenia, should always be taken into account.

Turbulence in the presence of internal waves in the bottom boundary layer of the California inner shelf

NASA Astrophysics Data System (ADS)

Allen, Rachel M.; Simeonov, Julian A.; Calantoni, Joseph; Stacey, Mark T.; Variano, Evan A.

2018-05-01

Turbulence measurements were collected in the bottom boundary layer of the California inner shelf near Point Sal, CA, for 2 months during summer 2015. The water column at Point Sal is stratified by temperature, and internal bores propagate through the region regularly. We collected velocity, temperature, and turbulence data on the inner shelf at a 30-m deep site. We estimated the turbulent shear production ( P), turbulent dissipation rate ( ɛ), and vertical diffusive transport ( T), to investigate the near-bed local turbulent kinetic energy (TKE) budget. We observed that the local TKE budget showed an approximate balance ( P ≈ ɛ) during the observational period, and that buoyancy generally did not affect the TKE balance. On a finer resolution timescale, we explored the balance between dissipation and models for production and observed that internal waves did not affect the balance in TKE at this depth.

Retention and diffusion of H, He, O, C impurities in Be

NASA Astrophysics Data System (ADS)

Zhang, Pengbo; Zhao, Jijun; Wen, Bin

2012-04-01

We report the energetics and diffusion behavior of H, He, O, and C impurities in beryllium as fusion materials from first-principles calculations. Among the six interstitial sites in Be, the basal tetrahedral one is most stable for H, He, O, while C prefers to occupy an octahedral site. Solution of O impurity in Be is an exothermic process with solution energy of -2.37 eV, whereas solution of H, C and He is an endothermic process (solution energy: 1.55 eV, 2.46 eV, and 5.70 eV, respectively). Overall speaking, these impurities prefer to diffuse along longer paths. The H and O impurities share the same out-of-plane diffusion path via basal tetrahedral sites, while the He and C impurities in Be mainly diffuse via basal tetrahedral and octahedral sites along the (0 0 1) plane. Diffusion of He in Be is easiest with a lowest barrier of 0.14 eV; whereas H diffusion in Be is also rather fast with migration energies of 0.4 eV. On the contrary, diffusion of C and O impurities is more difficult because of strong bonding with lattice atoms and high energy barriers of 0.42 and 1.63 eV, respectively. Our theoretical results provide the fundamental parameters for understanding the impurity aggregation and bubble formation in early stage of irradiation damage.

A Novel Diffuse Fraction-Based Two-Leaf Light Use Efficiency Model: An Application Quantifying Photosynthetic Seasonality across 20 AmeriFlux Flux Tower Sites

NASA Astrophysics Data System (ADS)

Yan, Hao; Wang, Shao-Qiang; Yu, Kai-Liang; Wang, Bin; Yu, Qin; Bohrer, Gil; Billesbach, Dave; Bracho, Rosvel; Rahman, Faiz; Shugart, Herman H.

2017-10-01

Diffuse radiation can increase canopy light use efficiency (LUE). This creates the need to differentiate the effects of direct and diffuse radiation when simulating terrestrial gross primary production (GPP). Here, we present a novel GPP model, the diffuse-fraction-based two-leaf model (DTEC), which includes the leaf response to direct and diffuse radiation, and treats maximum LUE for shaded leaves (ɛmsh defined as a power function of the diffuse fraction (Df)) and sunlit leaves (ɛmsu defined as a constant) separately. An Amazonian rainforest site (KM67) was used to calibrate the model by simulating the linear relationship between monthly canopy LUE and Df. This showed a positive response of forest GPP to atmospheric diffuse radiation, and suggested that diffuse radiation was more limiting than global radiation and water availability for Amazon rainforest GPP on a monthly scale. Further evaluation at 20 independent AmeriFlux sites showed that the DTEC model, when driven by monthly meteorological data and MODIS leaf area index (LAI) products, explained 70% of the variability observed in monthly flux tower GPP. This exceeded the 51% accounted for by the MODIS 17A2 big-leaf GPP product. The DTEC model's explicit accounting for the impacts of diffuse radiation and soil water stress along with its parameterization for C4 and C3 plants was responsible for this difference. The evaluation of DTEC at Amazon rainforest sites demonstrated its potential to capture the unique seasonality of higher GPP during the diffuse radiation-dominated wet season. Our results highlight the importance of diffuse radiation in seasonal GPP simulation.