Theoretical insights on electron donor-acceptor interactions involving carbon dioxide

Electron donor-acceptor (EDA) interactions are widely involved in chemistry and their understanding is essential to design new technological applications in a variety of fields ranging from material sciences and chemical engineering to medicine. In this work, we study EDA complexes of carbon dioxide with ketones using several ab initio and Density Functional Theory methods. Energy contributions to the interaction energy have been analyzed in detail using both variational and perturbational treatments. Dispersion energy has been shown to play a key role in explaining the high stability of a non-conventional structure, which can roughly be described by a cooperative EDA interaction.

Analysis of adsorption of gases and vapors on nonporous graphitized thermal carbon black

In this paper we analyzed the adsorption of a large number of gases and vapors on graphitized thermal carbon black. The Henry constant was used to determine the adsorbate-adsorbent interaction energy, which is found to be a modest decreasing function of temperature. Analysis of the complete adsorption isotherm over a wider range of pressure yields information on the monolayer coverage concentration and the adsorbate-adsorbate interaction energy. Among the various equations tested, the Hill-de Boer equation accounting for BET-postulated multilayer formation describes well the adsorption isotherms of all adsorbates. On average, the adsorbate-adsorbate interaction energy in the adsorbed phase is less than that in the bulk phase, suggesting that the distance between adsorbed molecules in the first layer of the adsorbed phase is slightly less than the equilibrium distance between two adsorbate molecules in the bulk phase. This suggests that the first layer is in a compressed state, which is due to the attraction of the adsorbent surface. The monolayer concentration as determined from the fitting of the Hill-de Boer equation with experimental data is slightly larger than the values calculated from the molecular projection area, suggesting that molecules can be oriented such that a larger number of molecules can be accommodated on the carbon black surface. This further supports the shorter distance between adsorbate molecules in the adsorbed phase.

Effects of potential models in the vapor-liquid equilibria and adsorption of simple gases on graphitized thermal carbon black

In this paper, we investigate the effects of various potential models in the description of vapor–liquid equilibria (VLE) and adsorption of simple gases on highly graphitized thermal carbon black. It is found that some potential models proposed in the literature are not suitable for the description of VLE (saturated gas and liquid densities and the vapor pressure with temperature). Simple gases, such as neon, argon, krypton, xenon, nitrogen, and methane are studied in this paper. To describe the isotherms on graphitized thermal carbon black correctly, the surface mediation damping factor introduced in our recent publication should be used to calculate correctly the fluid–fluid interaction energy between particles close to the surface. It is found that the damping constant for the noble gases family is linearly dependent on the polarizability, suggesting that the electric field of the graphite surface has a direct induction effect on the induced dipole of these molecules. As a result of this polarization by the graphite surface, the fluid–fluid interaction energy is reduced whenever two particles are near the surface. In the case of methane, we found that the damping constant is less than that of a noble gas having the similar polarizability, while in the case of nitrogen the damping factor is much greater and this could most likely be due to the quadrupolar nature of nitrogen.

Adsorption of Lennard-Jones fluids in carbon slit pores of a finite length. A computer simulation study

The adsorption of simple Lennard-Jones fluids in a carbon slit pore of finite length was studied with Canonical Ensemble (NVT) and Gibbs Ensemble Monte Carlo Simulations (GEMC). The Canonical Ensemble was a collection of cubic simulation boxes in which a finite pore resides, while the Gibbs Ensemble was that of the pore space of the finite pore. Argon was used as a model for Lennard-Jones fluids, while the adsorbent was modelled as a finite carbon slit pore whose two walls were composed of three graphene layers with carbon atoms arranged in a hexagonal pattern. The Lennard-Jones (LJ) 12-6 potential model was used to compute the interaction energy between two fluid particles, and also between a fluid particle and a carbon atom. Argon adsorption isotherms were obtained at 87.3 K for pore widths of 1.0, 1.5 and 2.0 nm using both Canonical and Gibbs Ensembles. These results were compared with isotherms obtained with corresponding infinite pores using Grand Canonical Ensembles. The effects of the number of cycles necessary to reach equilibrium, the initial allocation of particles, the displacement step and the simulation box size were particularly investigated in the Monte Carlo simulation with Canonical Ensembles. Of these parameters, the displacement step had the most significant effect on the performance of the Monte Carlo simulation. The simulation box size was also important, especially at low pressures at which the size must be sufficiently large to have a statistically acceptable number of particles in the bulk phase. Finally, it was found that the Canonical Ensemble and the Gibbs Ensemble both yielded the same isotherm (within statistical error); however, the computation time for GEMC was shorter than that for canonical ensemble simulation. However, the latter method described the proper interface between the reservoir and the adsorbed phase (and hence the meniscus).

Correct procedures for the calculation of heats of adsorption for heterogeneous adsorbents from molecular simulation

Several procedures for calculating the heat of adsorption from Monte Carlo simulations for a heterogeneous adsorbent are presented. Simulations have been performed to generate isotherms for nitrogen at 77 K and methane at 273.15 K in graphitic slit pores of various widths. The procedures were then applied to calculate the heat of adsorption of an activated carbon with an arbitrary pore size distribution. The consistency of the different procedures shows them to be correct in calculating interaction energy contributions to the heat of adsorption. The currently favored procedure for this type of calculation, from the literature, is shown to be incorrect and in serious error when calculating the heat of adsorption of activated carbon.

The charge effect on the hindrance factors for diffusion and convection of a solute in pores II

The diffusion and convection of a solute suspended in a fluid across porous membranes are known to be reduced compared to those in a bulk solution, owing to the fluid mechanical interaction between the solute and the pore wall as well as steric restriction. If the solute and the pore wall are electrically charged, the electrostatic interaction between them could affect the hindrance to diffusion and convection. In this study, the transport of charged spherical solutes through charged circular cylindrical pores filled with an electrolyte solution containing small ions was studied numerically by using a fluid mechanical and electrostatic model. Based on a mean field theory, the electrostatic interaction energy between the solute and the pore wall was estimated from the Poisson-Boltzmann equation, and the charge effect on the solute transport was examined for the solute and pore wall of like charge. The results were compared with those obtained from the linearized form of the Poisson-Boltzmann equation, i.e.the Debye-Hückel equation. © 2012 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.

Charge effects on hindrance factors for diffusion and convection of solute in pores I

The transport of a spherical solute through a long circular cylindrical pore filled with an electrolyte solution is studied numerically, in the presence of constant surface charge on the solute and the pore wall. Fluid dynamic analyses were carried out to calculate the flow field around the solute in the pore to evaluate the drag coefficients exerted on the solute. Electrical potentials around the solute in the electrolyte solution were computed based on a mean-field theory to provide the interaction energy between the charged solute and the pore wall. Combining the results of the fluid dynamic and electrostatic analyses, we estimated the rate of the diffusive and convective transport of the solute across the pore. Although the present estimates of the drag coefficients on the solute suggest more than 10% difference from existing studies, depending on the radius ratio of the solute relative to the pore and the radial position of the solute center in the pore, this difference leads to a minor effect on the hindrance factors. It was found that even at rather large ion concentrations, the repulsive electrostatic interaction between the charged solute and the pore wall of like charge could significantly reduce the transport rate of the solute.

Fluid mechanical and electrostatic study on the osmotic flow through cylindrical pores

When two solutions differing in solute concentration are separated by a porous membrane, the osmotic pressure will generate a net volume flux of the suspending fluid across the membrane; this is termed osmotic flow. We consider the osmotic flow across a membrane with circular cylindrical pores when the solute and the pore walls are electrically charged, and the suspending fluid is an electrolytic solution containing small cations and anions. Under the condition in which the radius of the pores and that of the solute molecules greatly exceed those of the solvent as well as the ions, a fluid mechanical and electrostatic theory is introduced to describe the osmotic flow in the presence of electric charge. The interaction energy, including the electrostatic interaction between the solute and the pore wall, plays a key role in determining the osmotic flow. We examine the electrostatic effect on the osmotic flow and discuss the difference in the interaction energy determined from the nonlinear Poisson-Boltzmann equation and from its linearized equation (the Debye-Hückel equation).

Electrical charge effect on osmotic flow through pores

An electrostatic model for osmotic flow through circular cylindrical pores is developed to describe the reflection coefficient for the membrane transport in the presence of surface charges on the pore wall and the solute. For a spherical solute placed at an arbitrary radial position in the pore, the electrical potential was computed by a spectral element method applied to the Poisson-Boltzmann equation together with the condition of electrical neutrality. The interaction energy between the surface charges was used to estimate the osmotic reflection coefficient. The proposed model predicts that even for a small Debye length compared to the pore radius, the repulsive electrostatic interaction between the surface charges could significantly increase the osmotic flow through the pore.

Estudo in silico das interações da protease NS3-NS2B de DENV-2 com o inibidor peptídico Bz-nKRR-H com finalidades terapêuticas

Dengue virus is an important patogen that causes Dengue desease in all world, and belongs to Flavivirus gender. The virus consists of enveloped RNA with a single strand positive sense, 11Kb genome. The RNA is translated into a polyprotein precursor, wich is cleaved into 3 structural proteins (C, prM e E) and 7 non-structural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B e NS5). The NS3 is a multifunctional protein, that besides to promote the polyprotein precursor cleavage, also have NTPase, helicase and RTPase activity. The NS3 needs a hydrophilic segment of 40 residues from the transmembrane NS2B protein (who acts like cofator) to realize this functions. Actually, there's no vacines available on the market, and the treatment are just symptomatic. The tetrapeptide inhibitor Bz-Nle-Lys-Arg-Arg-H (Ki de 5,8-7,0 M) was showed as a potent inhibitor μ for NS3prot in Dengue virus. That is a inteligent alternative to treat the dengue desease. The present work aimed analyse the interactions of the ligand bounded to the activity site to provid a clear and depth vision of that interaction. For this purpouse, it was conducted an in silico study, by using quantum mechanical calculations based on Density Functional Theory (DFT), with Generalized Gradient approximation (GGA) to describe the effects of exchange and correlation. The interaction energy of each amino acid belonging to the binding site to the ligand was calculated the using the method of molecular fragmentation with conjugated caps (MFCC). Besides energy, we calculated the distances, types of molecular interactions and atomic groups involved. The theoretical models used were satisfactory and show a more accurate description when the dielectric constant = 20 ε and 80 was used. The results demonstrate that the interaction energy of the system reached convergence at 13.5 A. Within a radius of 13,5A the most important residues were identified. Met49, Met84 and Asp81 perform interactions of hydrogen with the ligant. The Asp79 and Asp75 residues present high energy of attraction. Arg54, Arg85 and Lys 131 perform hydrogen interactions with the ligand, however, appear in BIRD graph having high repulsion energy with the inhibitor. The data also emphasizes the importance of residue Tyr161 and the involvement of the catalytic triad composed by Asp75, His51 and Ser135

Generalized atoms in molecules approach (GAIM approach)

This thesis involves two parts. The first is a new-proposed theoretical approach called generalized atoms in molecules (GAIM). The second is a computational study on the deamination reaction of adenine with OH⁻/nH₂O (n=0, 1, 2, 3) and 3H₂O. The GAIM approach aims to solve the energy of each atom variationally in the first step and then to build the energy of a molecule from each atom. Thus the energy of a diatomic molecule (A-B) is formulated as a sum of its atomic energies, EA and EB. Each of these atomic energies is expressed as, EA = Hᴬ + Vₑₑᴬᴬ + 1/2Vₑₑᴬ<>ᴮ EB = Hᴮ + Vₑₑᴮᴮ + 1/2Vₑₑᴬ<>ᴮ where; Hᴬ and Hᴮ are the kinetic and nuclear attraction energy of electrons of atoms A and B, respectively; Vₑₑᴬᴬ and Vₑₑᴮᴮ are the interaction energy between the electrons on atoms A and B, respectively; and Vₑₑᴬ<>ᴮ is the interaction energy between the electrons of atom A with the electrons of atom B. The energy of the molecule is then minimized subject to the following constraint, |ρA(r)dr + |ρB(r)dr = N where ρA(r) and ρB(r) are the electron densities of atoms A and B, respectively, and N is the number of electrons. The initial testing of the performance of GAIM was done through calculating dissociation curves for H₂, LiH, Li₂, BH, HF, HCl, N₂, F₂, and Cl₂. The numerical results show that GAIM performs very well with H₂, LiH, Li₂, BH, HF, and HCl. GAIM shows convergence problems with N₂, F₂, and Cl₂ due to difficulties in reordering the degenerate atomic orbitals Pₓ, Py, and Pz in N, F, and Cl atoms. Further work for the development of GAIM is required. Deamination of adenine results in one of several forms of premutagenic lesions occurring in DNA. In this thesis, mechanisms for the deamination reaction of adenine with OH⁻/nH₂O, (n = 0, 1, 2, 3) and 3H₂O were investigated. HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), and B3LYP/6-31+G(d) levels of theory were employed to optimize all the geometries. Energies were calculated at the G3MP2B3 and CBS-QB3 levels of theory. The effect of solvent (water) was computed using the polarizable continuum model (PCM). Intrinsic reaction coordinate (IRC) calculations were performed for all transition states. Five pathways were investigated for the deamination reaction of adenine with OH⁻/nH₂O and 3H₂O. The first four pathways (A-D) begin with by deprotonation at the amino group of adenine by OH⁻, while pathway E is initiated by tautomerization of adenine. For all pathways, the next two steps involve the formation of a tetrahedral intermediate followed by dissociation to yield products via a 1,3-hydrogen shift. Deamination with a single OH⁻ has a high activation barrier (190 kJ mol⁻¹ using G3MP2B3 level) for the rate-determining step. Addition of one water molecule reduces this barrier by 68 kJ mol⁻¹ calculated at G3MP2B3 level. Adding more water molecules decreases the overall activation energy of the reaction, but the effect becomes smaller with each additional water molecule. The most plausible mechanism is pathway E, the deamination reaction of adenine with 3H₂O, which has an overall G3MP2B3 activation energy of 139 and 137 kJ mol⁻¹ in the gas phase and PCM, respectively. This barrier is lower than that for the deamination with OH⁻/3H₂O by 6 and 2 kJ mol⁻¹ in the gas phase and PCM, respectively.

Analysis of Self-Organization

The dissertation is devoted to the study of problems in calculus of variation, free boundary problems and gradient flows with respect to the Wasserstein metric. More concretely, we consider the problem of characterizing the regularity of minimizers to a certain interaction energy. Minimizers of the interaction energy have a somewhat surprising relationship with solutions to obstacle problems. Here we prove and exploit this relationship to obtain novel regularity results. Another problem we tackle is describing the asymptotic behavior of the Cahn-Hilliard equation with degenerate mobility. By framing the Cahn-Hilliard equation with degenerate mobility as a gradient flow in Wasserstein metric, in one space dimension, we prove its convergence to a degenerate parabolic equation under the framework recently developed by Sandier-Serfaty.

Ação de anestésicos gerais em canais iônicos

Dissertação (mestrado)—Universidade de Brasília, Instituto de Ciências Biológicas, Departamento de Biologia Molecular, 2016.

Study and characterization of polymorphs of BTBT derivatives and their application in organic electronics

The perquisites of organic semiconductors (OSCs) in the field of organic electronics have attracted much attention due to the advantages like cost-effectiveness, solution processibility, etc. A key property in OSCs is charge carrier mobility, which depends on molecular packing, as even the slightest changes in the packing of OSC can significantly impact the mobility. Organic molecules are constructed by weak interactions, which makes the OSCs prone to adopt multiple packing arrangements, thus giving rise to polymorphism. Therefore, polymorph screening in bulk and thin films is crucial for material development. This thesis aims to present a systematic study of polymorphism of [1]benzothieno[3,2-b]benzothiophene (BTBT) derivatives functionalized with different side chains. The role of peripheral side chains has been studied since they can promote different packing arrangements. The bulk polymorph screening of OSCs was approached with conventional solution mediated recrystallization experiments like evaporation, slurry maturation, anti-solvent precipitation, etc. Each of the polymorphs were inspected for their relative stability and the kinetics of transformation was evaluated. Polymorphism in thin films was also investigated for selected OSCs. Non-equilibrium methods like, thermal gradient and solution shearing were employed to examine the nucleation, crystal growth and morphology in controlled crystallization conditions. After careful analysis of crystal phases in bulk and thin films, OFETs have been fabricated by optimizing the manufacturing conditions and the hole mobility values were extracted. The charge transport property of the OSCs tested for OFETs was supported by the ionization potential and transfer integrals calculation. An attempt to correlate the solid-state structure to electronic properties was carried out. For some of the molecules, mechanical properties have been also investigated, as the response to mechanical stress is highly susceptible to packing arrangements and the intermolecular interaction energy contributions. Additionally, collaborative research was carried out by solving and analysing the crystal structures of six oligorylene molecules.

Enhancing predictive capability of models for solubility and permeability in polymers and composites

The interpretation of phase equilibrium and mass transport phenomena in gas/solvent - polymer system at molten or glassy state is relevant in many industrial applications. Among tools available for the prediction of thermodynamics properties in these systems, at molten/rubbery state, is the group contribution lattice-fluid equation of state (GCLF-EoS), developed by Lee and Danner and ultimately based on Panayiotou and Vera LF theory. On the other side, a thermodynamic approach namely non-equilibrium lattice-fluid (NELF) was proposed by Doghieri and Sarti to consistently extend the description of thermodynamic properties of solute polymer systems obtained through a suitable equilibrium model to the case of non-equilibrium conditions below the glass transition temperature. The first objective of this work is to investigate the phase behaviour in solvent/polymer at glassy state by using NELF model and to develop a predictive tool for gas or vapor solubility that could be applied in several different applications: membrane gas separation, barrier materials for food packaging, polymer-based gas sensors and drug delivery devices. Within the efforts to develop a predictive tool of this kind, a revision of the group contribution method developed by High and Danner for the application of LF model by Panayiotou and Vera is considered, with reference to possible alternatives for the mixing rule for characteristic interaction energy between segments. The work also devotes efforts to the analysis of gas permeability in polymer composite materials as formed by a polymer matrix in which domains are dispersed of a second phase and attention is focused on relation for deviation from Maxwell law as function of arrangement, shape of dispersed domains and loading.