Effect of charge asymmetry and charge screening on structure of superlattices formed by oppositely charged colloidal particles

Colloidal suspensions made up of oppositely charged particles have been shown to self-assemble into substitutionally ordered superlattices. For a given colloidal suspension, the structure of the superlattice formed from self-assembly depends on its composition, charges on the particles, and charge screening. In this study we have computed the pressure-composition phase diagrams of colloidal suspensions made up of binary mixtures of equal sized and oppositely charged particles interacting via hard core Yukawa potential for varying values of charge screening and charge asymmetry. The systems are studied under conditions where the thermal energy is equal or greater in magnitude to the contact energy of the particles and the Debye screening length is smaller than the size of the particles. Our studies show that charge asymmetry has a significant effect on the ability of colloidal suspensions to form substitutionally ordered superlattices. Slight deviations of the charges from the stoichiometric ratio are found to drastically reduce the thermodynamic stability of substitutionally ordered superlattices. These studies also show that for equal-sized particles, there is an optimum amount of charge screening that favors the formation of substitutionally ordered superlattices. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.3700226]

Experimental and simulation studies on the performance of standing wave thermoacoustic prime mover for pulse tube refrigerator

The thermoacoustic prime mover (TAPM) has gained considerable attention as a pressure wave generator to drive pulse tube refrigerator (PTR) due to no moving parts, reasonable efficiency, use of environmental friendly working fluids etc. To drive PTCs, lower frequencies (f) with larger pressure amplitudes (Delta P) are essential, which are affected by geometric and operating parameters of TAPM as well as working fluids. For driving PTRs, a twin standing wave TAPM is built and studied by using different working fluids such as helium, argon, nitrogen and their binary mixtures. Simulation results of DeltaEc are compared with experimental data wherever possible. DeltaEc predicts slightly increased resonance frequencies, but gives larger Delta P and lower temperature difference Delta T across stack. High mass number working fluid leads to lower frequency with larger Delta P, but higher Delta T. Studies indicate that the binary gas mixture of right composition with lower Delta T can be arrived at to drive TAPM of given geometry. (C) 2013 Elsevier Ltd and IIR. All rights reserved.

Combined toxicity of two crystalline phases (anatase and rutile) of Titania nanoparticles towards freshwater microalgae: Chlorella sp

In view of the increasing usage of anatase and rutile crystalline phases of titania NPs in the consumer products, their entry into the aquatic environment may pose a serious risk to the ecosystem. In the present study, the possible toxic impact of anatase and rutile nanoparticles (individually and in binary mixture) was investigated using freshwater microalgae, Chlorella sp. at low exposure concentrations (0.25, 0.5 and 1 mg/L) in freshwater medium under UV irradiation. Reduction of cell viability as well as a reduction in chlorophyll content were observed due to the presence of NPs. An antagonistic effect was noted at certain concentrations of binary mixture such as (0.25, 0.25), (0.25, 0.5), and (0.5, 0.5) mg/L, and an additive effect for the other combinations, (0.25, 1), (0.5, 0.25), (0.5, 1), (1, 0.25), (1, 0.5), and (1, 1) mg/L. The hydrodynamic size analyses in the test medium revealed that rutile NPs were more stable in lake water than the anatase and binary mixtures at 6 h, the sizes of anatase (1 mg/L), rutile NPs (1 mg/L), and binary mixture (1, 1 mg/L) were 948.83 +/- 35.01 nm, 555.74 +/- 19.93 nm, and 1620.24 +/- 237.87 nm, respectively]. The generation of oxidative stress was found to be strongly dependent on the crystallinity of the nanoparticles. The transmission electron microscopic images revealed damages in the nucleus and cell membrane of algal cells due to the interaction of anatase NPs, whereas rutile NPs were found to cause chloroplast and internal organelle damages. Mis-shaped chloroplasts, lack of nucleus, and starch-pyrenoid complex were noted in binary-treated cells. The findings from the current study may facilitate the environmental risk assessment of titania NPs in an aquatic ecosystem. (C) 2015 Elsevier B.V. All rights reserved.

Spatio-temporal correlations in aqueous systems: computational studies of static and dynamic heterogeneity by 2D-IR spectroscopy

Local heterogeneity is ubiquitous in natural aqueous systems. It can be caused locally by external biomolecular subsystems like proteins, DNA, micelles and reverse micelles, nanoscopic materials etc., but can also be intrinsic to the thermodynamic nature of the aqueous solution itself (like binary mixtures or at the gas-liquid interface). The altered dynamics of water in the presence of such diverse surfaces has attracted considerable attention in recent years. As these interfaces are quite narrow, only a few molecular layers thick, they are hard to study by conventional methods. The recent development of two dimensional infra-red (2D-IR) spectroscopy allows us to estimate length and time scales of such dynamics fairly accurately. In this work, we present a series of interesting studies employing two dimensional infra-red spectroscopy (2D-IR) to investigate (i) the heterogeneous dynamics of water inside reverse micelles of varying sizes, (ii) supercritical water near the Widom line that is known to exhibit pronounced density fluctuations and also study (iii) the collective and local polarization fluctuation of water molecules in the presence of several different proteins. The spatio-temporal correlation of confined water molecules inside reverse micelles of varying sizes is well captured through the spectral diffusion of corresponding 2D-IR spectra. In the case of supercritical water also, we observe a strong signature of dynamic heterogeneity from the elongated nature of the 2D-IR spectra. In this case the relaxation is ultrafast. We find remarkable agreement between the different tools employed to study the relaxation of density heterogeneity. For aqueous protein solutions, we find that the calculated dielectric constant of the respective systems unanimously shows a noticeable increment compared to that of neat water. However, the `effective' dielectric constant for successive layers shows significant variation, with the layer adjacent to the protein having a much lower value. Relaxation is also slowest at the surface. We find that the dielectric constant achieves the bulk value at distances more than 3 nm from the surface of the protein.

Suspensions of polymer-grafted nanoparticles with added polymers-Structure and effective pair-interactions

We present the results of combined experimental and theoretical (molecular dynamics simulations and integral equation theory) studies of the structure and effective interactions of suspensions of polymer grafted nanoparticles (PGNPs) in the presence of linear polymers. Due to the absence of systematic experimental and theoretical studies of PGNPs, it is widely believed that the structure and effective interactions in such binary mixtures would be very similar to those of an analogous soft colloidal material-star polymers. In our study, polystyrene-grafted gold nanoparticles with functionality f = 70 were mixed with linear polystyrene (PS) of two different molecular weights for obtaining two PGNP: PS size ratios, xi = 0.14 and 2.76 (where, xi = M-g/M-m, M-g and M-m being the molecular weights of grafting and matrix polymers, respectively). The experimental structure factor of PGNPs could be modeled with an effective potential (Model-X), which has been found to be widely applicable for star polymers. Similarly, the structure factor of the blends with xi = 0.14 could be modeled reasonably well, while the structure of blends with xi = 2.76 could not be captured, especially for high density of added polymers. A model (Model-Y) for effective interactions between PGNPs in a melt of matrix polymers also failed to provide good agreement with the experimental data for samples with xi = 2.76 and high density of added polymers. We tentatively attribute this anomaly in modeling the structure factor of blends with xi = 2.76 to the questionable assumption of Model-X in describing the added polymers as star polymers with functionality 2, which gets manifested in both polymer-polymer and polymer-PGNP interactions especially at higher fractions of added polymers. The failure of Model-Y may be due to the neglect of possible many-body interactions among PGNPs mediated by matrix polymers when the fraction of added polymers is high. These observations point to the need for a new framework to understand not only the structural behavior of PGNPs but also possibly their dynamics and thermo-mechanical properties as well. (C) 2015 AIP Publishing LLC.

Influence Of Gravity On Structure Of Colloidal Crystal Using Simulated Microgravity

Liquid mixtures of water and deuterium oxide as the liquid phase, were used to match the density of charged colloidal particles. Kossel diffraction method was used to detect the crystal structures. The experiments under the density-matched (g=0) and unmatched (g=1) conditions are compared to examine the influence of gravity on the crystal structures formed by self-assembly of 110 nm (in diameter) polystyrene microspheres. The result shows that die gravity tends to make the lattice constants of colloidal crystals smaller at lower positions, which indicates that the effect of gravity should be taken into account in the study of the colloidal crystals.

Engineering thermostable fungal cellobiohydrolases

Meeting the world's growing energy demands while protecting our fragile environment is a challenging issue. Second generation biofuels are liquid fuels like long-chain alcohols produced from lignocellulosic biomass. To reduce the cost of biofuel production, we engineered fungal family 6 cellobiohydrolases (Cel6A) for enhanced thermostability using random mutagenesis and recombination of beneficial mutations. During long-time hydrolysis, engineered thermostable cellulases hydrolyze more sugars than wild-type Cel6A as single enzymes and binary mixtures at their respective optimum temperatures. Engineered thermostable cellulases exhibit synergy in binary mixtures similar to wild-type cellulases, demonstrating the utility of engineering individual cellulases to produce novel thermostable mixtures. Crystal structures of the engineered thermostable cellulases indicate that the stabilization comes from improved hydrophobic interactions and restricted loop conformations by proline substitutions. At high temperature, free cysteines contribute to irreversible thermal inactivation in engineered thermostable Cel6A and wild-type Cel6A. The mechanism of thermal inactivation in this cellulase family is consistent with disulfide bond degradation and thiol-disulfide exchange. Enhancing the thermostability of Cel6A also increases tolerance to pretreatment chemicals, demonstrated by the strong correlation between thermostability and tolerance to 1-ethyl-3-methylimidazolium acetate. Several semi-rational protein engineering approaches – on the basis of consensus sequence analysis, proline stabilization, FoldX energy calculation, and high B-factors – were evaluated to further enhance the thermostability of Cel6A.

Modelagem do Equilíbrio de Fases entre Hidrocarbonetos Leves e Pesados

O comportamento de fases para sistemas binários com um hidrocarboneto leve e um pesado é muito importante tanto para o projeto real de um processo quanto para o desenvolvimento de modelos teóricos. Para atender a crescente demanda por informação experimental de equilíbrio de fases a altas pressões, o objetivo deste estudo é obter uma metodologia que substitua parcialmente ou maximize a pouca informação experimental disponível. Para isto propõe-se a modelagem do equilíbrio de fases em misturas de hidrocarboneto leve com um pesado, sem o conhecimento da estrutura molecular do pesado, inferindo-se os parâmetros do modelo a partir da modelagem de dados de ponto de bolha obtidos na literatura. Esta metodologia implica não só na descrição do equilíbrio de fases de um sistema como na estimação das propriedades críticas do pesado, de difícil obtenção devido ao craqueamento destes a altas temperaturas. Neste contexto, este estudo apresenta uma estratégia que estima indiretamente as propriedades críticas dos compostos pesados. Para isto, foram correlacionados dados experimentais de ponto de bolha de misturas binárias contendo um hidrocarboneto leve e um pesado, usando-se dois modelos: o de Peng-Robinson e o TPT1M (Teoria da Polimerização Termodinâmica de primeira ordem de Wertheim modificada). Os parâmetros ajustados com o modelo de Peng-Robinson correspondem diretamente às propriedades críticas do composto pesado, enquanto os ajustados com o modelo TPT1M foram usados para obtê-las. Esta estratégia fornece parâmetros dependentes do modelo, porém permite o cálculo de outras propriedades termodinâmicas, como a extrapolação da temperatura dos dados estudados. Além disso, acredita-se que a correlação dos parâmetros obtidos com as propriedades críticas disponíveis ajudará na caracterização de frações pesadas de composição desconhecida

I. Electronic states of heavily-doped molecular crystals--naphthalene. I. Theoretical. II. Experimental. II. Lowest singlet and triplet excited states of pyrazine

PART I

The energy spectrum of heavily-doped molecular crystals was treated in the Green’s function formulation. The mixed crystal Green’s function was obtained by averaging over all possible impurity distributions. The resulting Green’s function, which takes the form of an infinite perturbation expansion, was further approximated by a closed form suitable for numerical calculations. The density-of-states functions and optical spectra for binary mixtures of normal naphthalene and deuterated naphthalene were calculated using the pure crystal density-of-state functions. The results showed that when the trap depth is large, two separate energy bands persist, but when the trap depth is small only a single band exists. Furthermore, in the former case it was found that the intensities of the outer Davydov bands are enhanced whereas the inner bands are weakened. Comparisons with previous theoretical calculations and experimental results are also made.

PART II

The energy states and optical spectra of heavily-doped mixed crystals are investigated. Studies are made for the following binary systems: (1) naphthalene-h₈ and d₈, (2) naphthalene--h₈ and αd₄, and (3) naphthalene--h₈ and βd₁, corresponding to strong, medium and weak perturbations. In addition to ordinary absorption spectra at 4˚K, band-to-band transitions at both 4˚K and 77˚K are also analyzed with emphasis on their relations to cooperative excitation and overall density-of-states functions for mixed crystals. It is found that the theoretical calculations presented in a previous paper agree generally with experiments except for cluster states observed in system (1) at lower guest concentrations. These features are discussed semi-quantitatively. As to the intermolecular interaction parameters, it is found that experimental results compare favorably with calculations based on experimental density-of-states functions but not with those based on octopole interactions or charge-transfer interactions. Previous experimental results of Sheka and the theoretical model of Broude and Rashba are also compared with present investigations.

PART III

The phosphorescence, fluorescence and absorption spectra of pyrazine-h₄ and d₄ have been obtained at 4˚K in a benzene matrix. For comparison, those of the isotopically mixed crystal pyrazine-h₄ in d₄ were also taken. All these spectra show extremely sharp and well-resolved lines and reveal detailed vibronic structure.

The analysis of the weak fluorescence spectrum resolves the long-disputed question of whether one or two transitions are involved in the near-ultraviolet absorption of pyrazine. The “mirror-image relationship” between absorption and emission shows that the lowest singlet state is an allowed transition, properly designated as ¹B_3u ← ¹A_1g. The forbidden component ¹B_2g, predicted by both “exciton” and MO theories to be below the allowed component, must lie higher. Its exact location still remains uncertain.

The phosphorescence spectrum when compared with the excitation phosphorescence spectra, indicates that the lowest triplet state is also symmetry allowed, showing a strong 0-0 band and a “mirror-image relationship” between absorption and emission. In accordance with previous work, the triplet state is designated as ³B_3u.

The vibronic structure of the phosphorescence spectrum is very complicated. Previous work on the analysis of this spectrum all concluded that a long progression of v_6a exists. Under the high resolution attainable in our work, the supposed v_6a progression proves to have a composite triplet structure, starting from the second member of the progression. Not only is the v_9a hydrogen-bending mode present as shown by the appearance of the C-D bending mode in the d₄ spectrum, but a band of 1207 cm^-1 in the pyrazine in benzene system and 1231 cm^-1 in the mixed crystal system is also observed. This band is assigned as 2v_6b and of a_1g symmetry. Its anonymously strong intensity in the phosphorescence spectrum is interpreted as due to the Fermi resonance with the 2v_6a and v_9a band.

To help resolve the present controversy over the crystal phosphorescence spectrum of pyrazine, detailed vibrational analyses of the emission spectra were made. The fluorescence spectrum has essentially the same vibronic structure as the phosphorescence spectrum.

The separation of atomic and molecular gases from helium by cataphoresis

The cataphoretic purification of helium was investigated for binary mixtures of He with Ar, Ne, N₂, O₂, CO, and CO₂ in DC glow discharge. An experimental technique was developed to continuously measure the composition in the anode end-bulb without sample withdrawal. Discharge currents ranged from 10 ma to 100 ma. Total gas pressure ranged from 2 torr to 9 torr. Initial compositions of the minority component in He ranged from 1.2 mole percent to 7.5 mole percent.

The cataphoretic separation of Ar and Ne from He was found to be in agreement with previous investigators. The cataphoretic separation of N₂, O₂, and CO from He was found to be similar to noble gas systems in that the steady-state separation improved with (1) increasing discharge current, (2) increasing gas pressure, and (3) decreasing initial composition of the minority component. In the He-CO₂ mixture, the CO₂ dissociated to CO plus O₂. The fraction of CO₂ dissociated was directly proportional to the current and pressure and independent of initial composition.

The experimental results for the separation of Ar, Ne, N₂, O₂, and CO from He were interpreted in the framework of a recently proposed theoretical model involving an electrostatic Peclet number. In the model the electric field was assumed to be constant. This assumption was checked experimentally and the maximum variation in electric field was 35% in time and 30% in position. Consequently, the assumption of constant electric field introduced no more than 55% variation in the electrostatic Peclet number during a separation.

To aid in the design of new cataphoretic systems, the following design criteria were developed and tested in detail: (1) electric field independent of discharge current, (2) electric field directly proportional to total pressure, (3) ion fraction of impurity directly proportional to discharge current, and (4) ion fraction of impurity independent of total pressure. Although these assumptions are approximate, they enabled the steady-state concentration profile to be predicted to within 25% for 75% of the data. The theoretical model was also tested with respect to the characteristic time associated with transient cataphoresis. Over 80% of the data was within a factor of two of the calculated characteristic times.

The electrostatic Peclet number ranged in value from 0.13 to 4.33. Back-calculated ion fractions of the impurity component ranged in value from 4.8x10^-6 to 178x10^-6.

Modelagem da tensão superficial de misturas de hidrocarbonetos

Esta dissertação apresenta um novo modelo criado para a previsão do valor da tensão superficial de misturas de hidrocarbonetos. Pode ser utilizado na forma preditiva ou correlativa. Este modelo foi baseado no modelo preditivo para soluções ideais de HildebrandScott (1950) e no modelo correlativo de JouybanAcree (2004). Outras regras de mistura foram avaliadas e comparadas com o modelo proposto, incluindo abordagens clássicas como de Sprow e Prausnitz, parácoro e estados correspondentes. Foram avaliados os desempenhos dos modelos escolhidos na previsão e na correlação da tensão superficial de misturas de hidrocarbonetos em diferentes temperaturas que continham diferentes assimetrias na mistura, tanto em natureza química (e.g., parafínicos, naftênicos e aromáticos) quanto em tamanho de cadeia. No total foram empregados 773 dados experimentais. Dentre os modelos preditivos, o modelo proposto foi o de melhor desempenho. Para sistemas com comportamento próximo da idealidade, esse modelo apresentou desvios em relação aos dados experimentais compatíveis com o erro experimental. Todos os modelos correlativos estudados apresentaram desvios essencialmente da ordem do erro experimental. O modelo proposto é importante, pois é capaz de prever a tensão superficial de misturas binárias aproximadamente ideais de hidrocarbonetos com elevada acurácia, além de permitir a correlação eficiente dos dados de todas as misturas estudadas

Inversão de funções do plano no plano aplicada ao cálculo de azeótropos

Azeotropia é um fenômeno termodinâmico onde um líquido em ebulição produz um vapor com composição idêntica. Esta situação é um desafio para a Engenharia de Separação, já que os processos de destilação exploram as diferenças entre as volatilidades relativas e, portanto, um azeótropo pode ser uma barreira para a separação. Em misturas binárias, o cálculo da azeotropia é caracterizado por um sistema não-linear do tipo 2 × 2. Um interessante e raro caso é o denominado azeotropia dupla, que pode ser verificado quando este sistema não-linear tem duas soluções, correspondendo a dois azeótropos distintos. Diferentes métodos tem sido utilizados na resolução de problemas desta natureza, como métodos estocásticos de otimização e as técnicas intervalares (do tipo Newton intervalar/bisseção generalizada). Nesta tese apresentamos a formulação do problema de azeotropia dupla e uma nova e robusta abordagem para a resolução dos sistemas não-lineares do tipo 2 × 2, que é a inversão de funções do plano no plano (MALTA; SALDANHA; TOMEI, 1996). No método proposto, as soluções são obtidas através de um conjunto de ações: obtenção de curvas críticas e de pré-imagens de pontos arbritários, inversão da função e por fim, as soluções esperadas para o problema de azeotropia. Esta metodologia foi desenvolvida para resolver sistemas não-lineares do tipo 2 × 2, tendo como objetivo dar uma visão global da função que modela o fenômeno em questão, além, é claro, de gerar as soluções esperadas. Serão apresentados resultados numéricos para o cálculo dos azeótropos no sistema benzeno + hexafluorobenzeno a baixas pressões por este método de inversão. Como ferramentas auxiliares, serão também apresentados aspectos numéricos usando aproximações clássicas, tais como métodos de Newton com técnicas de globalização e o algorítmo de otimização não-linear C-GRASP, para efeito de comparação.

Evaluation of estrogenic activities and mechanism of action of perfluorinated chemicals determined by vitellogenin induction in primary cultured tilapia hepatocytes

Perfluorochemicals (PFCs) are emerging persistent organic pollutants (POPs) and are widely present in the environment, wildlife and humans. Recently, reports have suggested that PFCs may have endocrine-disrupting activities. In the present study, we have developed a non-competitive enzyme-linked immunosorbent assay (ELISA) method to investigate estrogenic activities of selected PFCs using vitellogenin (VTG) induction in primary cultured hepatocytes of freshwater male tilapia (Oreochromis niloticus). Cultured hepatocytes were exposed to various concentrations of perfluorooctanyl sulfonate (PFOS), pentadecafluorooctanoic acid (PFOA), 1H, 1H, 2H, 2H-nonafluoro-1-hexanol (4:2 FTOH), 1H, 1H, 2H, 2H-perfluorooctanol (6:2 FTOH) and 1H, 1H, 2H, 2H-perfluoro-1-decanol (8:2 FTOH) for 48h, while 17 beta-estradiol (E2) and 4-nonylphenol (4-NP) were used as positive controls. A dose-dependent induction of VTG was observed in E2-, 4-NP-, PFOS-, PFOA- and 6:2 FrOH-treated cells, whereas VTG levels remained unchanged in the 4:2 FTOH and 8:2 FTOH exposure groups at the concentrations tested. The estimated 48-h EC50 values for E2,4-NP, PFOS, PFOA and 6:2 FTOH were 4.7 x 10(-7), 7.1 x 10(-6), 1.5 x 10(-5), 2.9 x 10(-5) and 2.8 x 10(-5) M, respectively. In the time-course study, significant VTG induction took place at 24 h (E2), 6 It (4-NP), 48 It (PFOS), 48 It (PFOA), 72 It (4:2 FTOH), 12 h (6:2 FTOH), 72 h (8:2 FTOH), and increased further after 96 It of exposure. Co-exposure to binary mixtures of individual PFCs and E2 for 48 It significantly inhibited E2-induced hepatocellular VTG production in a dose-dependent manner except for 4:2 FTOH. The estimated 48-h IC50 (concentration of a compound that elicits 50% inhibition of maximally E2-induced VTG) values for PFOS, PFOA, 6:2 FTOH and 8:2 FTOH were 3.1 x 10(-7), 5.1 X 10(-7), 1.1 X 10(-6) and 7.5 x 10(-7) M, respectively. In order to further investigate the estrogenic mechanism of PFCs, the hepatocytes were co-exposed to binary mixtures of individual chemicals (E2,4-NP, PFOS, PFOA and 6:2 FTOH) and the known estrogen receptor inhibitor tamoxifen for 48 h; tamoxifen significantly inhibited the ability of these chemicals to stimulate vitellogenesis. The overall results demonstrated that PFOS, PFOA and FTOHs have estrogenic activities and that exposure to a combination of E2 and PFCs produced anti-estrogenic effects. The results of the estrogen receptor inhibition assay further suggested that the estrogenic effect of PFCs may be mediated by the estrogen receptor pathway in primary cultured tilapia hepatocytes. (c) 2007 Elsevier B.V. All rights reserved.

基于Lattice Boltzmann模型的液-液混合流模拟

引入了一种二元Lattice Boltzmann Model(LBM),实现了两种液体组成的混合流的模拟.不同于其它的类似模型,它区分考虑了流体的粘性和扩散特性,可以很容易地模拟各种互溶或者不互溶的混合流现象.此外,由于LBM的运算大都是线性的局部运算,这使得它很容易在可编程图形处理器(Graphics Process Unit,GPU)上进行加速,从而进行实时模拟.给出了若干二元混合流的模拟结果.

simulation and interaction of fluid dynamics

In the fluid simulation, the fluids and their surroundings may greatly change properties such as shape and temperature simultaneously, and different surroundings would characterize different interactions, which would change the shape and motion of the fluids in different ways. On the other hand, interactions among fluid mixtures of different kinds would generate more comprehensive behavior. To investigate the interaction behavior in physically based simulation of fluids, it is of importance to build physically correct models to represent the varying interactions between fluids and the environments, as well as interactions among the mixtures. In this paper, we will make a simple review of the interactions, and focus on those most interesting to us, and model them with various physical solutions. In particular, more detail will be given on the simulation of miscible and immiscible binary mixtures. In some of the methods, it is advantageous to be taken with the graphics processing unit (GPU) to achieve real-time computation for middle-scale simulation.