Estimating successive pK(a) values of polyprotic acids from ab initio molecular dynamics using metadynamics: the dissociation of phthalic acid and its isomers

Estimation of the dissociation constant, or pK(a), of weak acids continues to be a central goal in theoretical chemistry. Here we show that ab initio Car-Parrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free energy profile of the dissociation reaction can provide reasonable estimates of the successive pK(a) values of polyprotic acids. We use the distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic group as the collective variable to explore the free energy profile of the dissociation process. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. Two distinct minima corresponding to the dissociated and un-dissociated states of the acid are observed and the difference in their free energy values provides the estimate for pK(a), the acid dissociation constant. We show that the method predicts the pK(a) value of benzoic acid in good agreement with experiment and then show using phthalic acid (benzene dicarboxylic acid) as a test system that both the first and second pK(a) values as well, as the subtle difference in their values for different isomers can be predicted in reasonable agreement with experimental data.

Mechanical properties of CNT-Bisphenol E cyanate ester-based CFRP nanocomposite developed through VARTM process

This paper reports on the effect of multiwall carbon nanotubes (CNTs) without and with chemical functionalization on the mechanical properties of Bisphenol E cyanate ester resin (BECy) based carbon fibre reinforced plastic (CFRP) laminated composites. BECy with its low viscosity, low moisture uptake and superior mechanical properties is selected for its application in CFRP laminates through the cost-effective Vacuum Assisted Resin Transfer Moulding (VARTM) process. However, unlike CNT-epoxy-CFRP composites, processing and performance issues such as dispersion of CNTs, chemical bonding with resin, functionalization effects, effects on mechanical properties, etc. for BECy-CNT-CFRP composite system are not well reported. The objective of this study is to improve the mechanical properties of BECy resin with small additions of CNTs and functionalized CNTs in CFRP laminates. CNTs and fCNTs are infused into BECy using ultrasonication and standard mixing methods. Improvements in Young's modulus and strength in tension, compression, shear, flexure and interlaminar shear strength are analysed. It is observed that addition of 0.5wt% CNTs effected for maximum mechanical properties of the resin and 1wt% CNTs for the mechanical properties of CNT-CFRP nanocomposite. Further, improvements obtained with fCNTs are marginal. Dispersion behaviour and effect of CNTs/fCNTs in load transfer corroborated with SEM pictures are presented. The enhanced mechanical properties realized in VARTM processing of BECy-CFRP laminate indicate clear advantage of CNT based modification of the process.

Bio-inspired route for the synthesis of spherical shaped MgO:Fe3+ nanoparticles: Structural, photoluminescence and photocatalytic investigation

MgO:Fe3+ (0.1-5 mol%) nanoparticles (NPs) were synthesized via eco-friendly, inexpensive and simple low temperature solution combustion route using Aloe vera gel as fuel. The final products were characterized by SEM, TEM and HRTEM. PXRD data and Rietveld analysis revealed the formation of cubic system. The influence of Fe3+ ion concentration on the structure morphology, UV absorption, PL emission and photocatalytic activity of MgO:Fe3+ NPs were investigated. The yellow emission with CIE chromaticity coordinates (0.44, 0.52) and average correlated color temperature value was found to be 3540 K which corresponds to warm light of NPs. The control of Fe3+. on MgO matrix influences the photocatalytic decolorization of methylene blue (MB) under UV light. The enhanced photocatalytic activity of MgO:Fe3+ (4 mol%) was attributed to dopant concentration, effective crystallite size, textural properties, decreased band gap and capability for reducing the electron hole pair recombination. Further, the trends of inhibitory effect in the presence of different radical scavengers were explored. These findings open up new avenues for the exploration of Fe-doped MgO in eco-friendly water applications and in the process of display devices. (C) 2015 Elsevier B.V. All rights reserved.

System reliability analysis of granular filter for protection against piping in dams

Granular filters are provided for the safety of water retaining structure for protection against piping failure. The phenomenon of piping triggers when the base soil to be protected starts migrating in the direction of seepage flow under the influence of seepage force. To protect base soil from migration, the voids in the filter media should be small enough but it should not also be too small to block smooth passage of seeping water. Fulfilling these two contradictory design requirements at the same time is a major concern for the successful performance of granular filter media. Since Terzaghi era, conventionally, particle size distribution (PSD) of granular filters is designed based on particle size distribution characteristics of the base soil to be protected. The design approach provides a range of D15f value in which the PSD of granular filter media should fall and there exist infinite possibilities. Further, safety against the two critical design requirements cannot be ensured. Although used successfully for many decades, the existing filter design guidelines are purely empirical in nature accompanied with experience and good engineering judgment. In the present study, analytical solutions for obtaining the factor of safety with respect to base soil particle migration and soil permeability consideration as proposed by the authors are first discussed. The solution takes into consideration the basic geotechnical properties of base soil and filter media as well as existing hydraulic conditions and provides a comprehensive solution to the granular filter design with ability to assess the stability in terms of factor of safety. Considering the fact that geotechnical properties are variable in nature, probabilistic analysis is further suggested to evaluate the system reliability of the filter media that may help in risk assessment and risk management for decision making.

Design for sustainability: case of designing an urban household organic waste management system

Beyond product design, if the notion of product `lifecycle design' enforces the consideration of requirements from all the lifecycle phases of products, design for sustainability enforces the consideration of lifecycle design in the context of the lifecycles of other products, processes, institutions and their design. Consequently, sustainability requirements that need to be met by design are very diverse. In this article, we portray the nature of design process to address sustainability requirements. This is done taking an example of designing a urban household organic waste management system that requires less water and reclaims the nutrients.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.

Dislocation theory of the fracture criterion for anisotropic solids

A dislocation theory of fracture criterion for the mixed dislocation emission and cleavage process in an anisotropic solid is developed in this paper. The complicated cases involving mixed-mode loading are considered here. The explicit formula for dislocations interaction with a semi-infinite crack is obtained. The governing equation for the critical condition of crack cleavage in an anisotropic solid after a number dislocation emissions is established. The effects of elastic anisotropy, crack geometry and load phase angle on the critical energy release rate and the total number of the emitted dislocations at the onset of cleavage are analysed in detail. The analyses revealed that the critical energy release rates can increase to one or two magnitudes larger than the surface energy because of the dislocation emission. It is also found elastic anisotropy and crystal orientation have significant effects on the critical energy release rates. The anisotropic values can be several times the isotropic value in one crack orientation. The values may be as much as 40% less than the isotropic value in another crack orientation and another anisotropy parameter. Then the theory is applied to a fee single crystal. An edge dislocation can emit from the crack tip along the most highly shear stressed slip plane. Crack cleavage can occur along the most highly stressed slip plane after a number of dislocation emissions. Calculation is carried out step by step. Each step we should judge by which slip system is the most highly shear stressed slip system and which slip system has the largest energy release rate. The calculation clearly shows that the crack orientation and the load phase angle have significant effects on the crystal brittle-ductile behaviours.

On nonequilibrium of pyrolysis process in the manufacture of ethene

Nonequilibrium process for cracking ethane and n-buthane in the manufacture of ethene has been analytically and numerically investigated in a Heavi-side function temperature field and through a normal shock wave. The results demonstrate that, while the reaction temperature increases, the maximum value of ethene yield is increased, and the optimal reaction duration is sharply shortened. For the identical initial reaction temperature, the maximum value of ethene yield through a stationary normal shock wave is less than that in a Heavi-side function temperature field. However, the ethene consumption after the maximum value in the former case is less than that in the latter. Higher ethene yield will be obtained by using the gasdynamic heating method than by using the current methods.

Analysis of ingredient and heating value of municipal solid waste

Great differences between municipal solid wastes(MSW)produced at different places and different times in terms of such parameters as physical ingredient and heating value lead to difficulty in effective handling of MSW. In this paper, ingredient, heating value and their temporal varying trends of typical MSW in Beijing were continuously measured and analyzed. With consideration of the process in pyrolysis and incineration, correlation between physical ingredients and heating values was induced, favorable for evaluation of heating value needed in handling of MSW from simple analysis of physical ingredients of it.

Effects of Baffle Design on Fluid Flow and Heat Transfer in Ammonothermal Growth System of Nitrides

Efforts have been made in growing bulk single crystals of GaN front supercritical fluids using the ammonothermal method, which utilizes ammonia as fluid rather than water as in the hydrothermal process. Different mineralizers such as amide or azide and temperatures in the range of 200-600degreesC have been used to increase the solubility. The pressure is from 1 to 4 kbar. Modeling of the ammonothermal growth process has been used to identify factors which may affect the temperature distribution, fluid flow and nutrient transport. The GaN charge is considered as a porous media bed and the flow in the charge is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. The effects of baffle design and opening on flow pattern and temperature distribution in an autoclave are analyzed. Two cases are considered with baffle openings of 15% and 20% in cross-sectional area, respectively.

The ensemble statistics of the energy of a random system subjected to harmonic excitation

This paper is concerned with the ensemble statistics of the response to harmonic excitation of a single dynamic system such as a plate or an acoustic volume. Random point process theory is employed, and various statistical assumptions regarding the system natural frequencies are compared, namely: (i) Poisson natural frequency spacings, (ii) statistically independent Rayleigh natural frequency spacings, and (iii) natural frequency spacings conforming to the Gaussian orthogonal ensemble (GOE). The GOE is found to be the most realistic assumption, and simple formulae are derived for the variance of the energy of the system under either point loading or rain-on-the-roof excitation. The theoretical results are compared favourably with numerical simulations and experimental data for the case of a mass loaded plate. © 2003 Elsevier Ltd. All rights reserved.

A study of the influence of free convection on coagulation process

It was assumed [1, 2] that gravity affects the coagulation process in two ways: free convection, which is hard to be avoided on the ground and sedimentation, which can be greatly reduced by the density-matching method. We present a ground-based experiment set-up to study the influence of convection on the perikinetic coagulation for aqueous polystyrene (PS) dispersions. The turbidity measurement was used to evaluate the relative coagulation rate and convection-driven flows in the solution were checked with a visual-magnification system. The pattern of flow field temperature profile in the sample cell is given. Our experiments show that there was no noticeable difference of coagulation rate observed no matter whether convection flows exist (with the flow speed up to 180 mu m/s) or not.

Determination Ofinterfacial Mechanical Parameters For An Al/Epoxy/Al2O3 System Byusing Peel Test Simulations

Peel test measurements and simulations of the interfacial mechanical parameters for the Al/Epoxy/Al2O3 system are performed in the present investigation. A series of Al film thicknesses between 20 and 250 microns and three peel angles of 90, 135 and 180 degrees are considered. Two types of epoxy adhesives are adopted to obtain both strong and weak interface adhesions. A finite element model with cohesive zone elements is used to identify the interfacial parameters and simulate the peel test process. By simulating and recording normal stress near the crack tip, the separation strength is obtained. Furthermore, the cohesive energy is identified by comparing the simulated steady-state peel force and the experimental result. It is found from the research that both the cohesive energy and the separation strength can be taken as the intrinsic interfacial parameters which are dependent on the thickness of the adhesive layer and independent of the film thickness and peel angle.

Study on flow characteristics of solid/liquid system in lysozyme crystal growth

During the process of lysozyme protein crystallization with batch method, the macroscopic flow field of solid/liquid system was observed by particle image velocimetry (PIV). Furthermore, a normal growth rate of (110) face and local flow field around a single protein crystal were obtained by a long work distance microscope. The experimental results showed that the average velocity, the maximal velocity of macroscopic solid/liquid system and the velocity of local flow field around single protein crystal were fluctuant. The effective boundary layer thickness delta(eff), the concentration at the interface Q and the characteristic velocity V were calculated using a convection-diffusion model. The results showed that the growth of lysozyme crystal in this experiment was dominated by interfacial kinetics rather than bulk transport, and the function of buoyancy-driven flow in bulk transport was small, however, the effect of bulk transport in crystal growth had a tendency to increase with the increase of lysozyme concentration. The calculated results, also showed that the order of magnitude of shear force was about 10(-21) N, which was much less than the bond force between the lysozyme molecules. Therefore the shear force induced by buoyancy-driven flows cannot remove the protein molecules from the interface of crystal.

A systems approach to characterise the dissemination process of household water treatment systems in developing countries

Abstract: Starting in the 1980s, household-level water treatment and safe storage systems (HWTS) have been developed as simple, local, user-friendly, and low cost options to improve drinking water quality at the point of use. However, despite conclusive evidence of the health and economic benefits of HWTS, and promotion efforts in over 50 countries in the past 20 years, implementation outcomes have been slow, reaching only 5-10 million regular users. This study attempts to understand the barriers and drivers affecting HWTS implementation. Although existing literature related to HWTS and innovation diffusion theories proposed ample critical factors and recommendations, there is a lack of holistic and systemic approach to integrate these findings. It is proposed that system dynamics modelling can be a promising tool to map the inter-relationships of different critical factors and to understand the structure of HWTS dissemination process, which may lead to identifying high impact, leveraged mitigation strategies to scale-up HWTS adoption and sustained use.