Proteins in silico-modeling and sampling

Proteins are linear chain molecules made out of amino acids. Only when they fold to their native states, they become functional. This dissertation aims to model the solvent (environment) effect and to develop & implement enhanced sampling methods that enable a reliable study of the protein folding problem in silico. We have developed an enhanced solvation model based on the solution to the Poisson-Boltzmann equation in order to describe the solvent effect. Following the quantum mechanical Polarizable Continuum Model (PCM), we decomposed net solvation free energy into three physical terms– Polarization, Dispersion and Cavitation. All the terms were implemented, analyzed and parametrized individually to obtain a high level of accuracy. In order to describe the thermodynamics of proteins, their conformational space needs to be sampled thoroughly. Simulations of proteins are hampered by slow relaxation due to their rugged free-energy landscape, with the barriers between minima being higher than the thermal energy at physiological temperatures. In order to overcome this problem a number of approaches have been proposed of which replica exchange method (REM) is the most popular. In this dissertation we describe a new variant of canonical replica exchange method in the context of molecular dynamic simulation. The advantage of this new method is the easily tunable high acceptance rate for the replica exchange. We call our method Microcanonical Replica Exchange Molecular Dynamic (MREMD). We have described the theoretical frame work, comment on its actual implementation, and its application to Trp-cage mini-protein in implicit solvent. We have been able to correctly predict the folding thermodynamics of this protein using our approach.

Growth, modification and integration of carbon nanotubes into molecular electronics

Molecules are the smallest possible elements for electronic devices, with active elements for such devices typically a few Angstroms in footprint area. Owing to the possibility of producing ultrahigh density devices, tremendous effort has been invested in producing electronic junctions by using various types of molecules. The major issues for molecular electronics include (1) developing an effective scheme to connect molecules with the present micro- and nano-technology, (2) increasing the lifetime and stabilities of the devices, and (3) increasing their performance in comparison to the state-of-the-art devices. In this work, we attempt to use carbon nanotubes (CNTs) as the interconnecting nanoelectrodes between molecules and microelectrodes. The ultimate goal is to use two individual CNTs to sandwich molecules in a cross-bar configuration while having these CNTs connected with microelectrodes such that the junction displays the electronic character of the molecule chosen. We have successfully developed an effective scheme to connect molecules with CNTs, which is scalable to arrays of molecular electronic devices. To realize this far reaching goal, the following technical topics have been investigated. 1. Synthesis of multi-walled carbon nanotubes (MWCNTs) by thermal chemical vapor deposition (T-CVD) and plasma-enhanced chemical vapor deposition (PECVD) techniques (Chapter 3). We have evaluated the potential use of tubular and bamboo-like MWCNTs grown by T-CVD and PE-CVD in terms of their structural properties. 2. Horizontal dispersion of MWCNTs with and without surfactants, and the integration of MWCNTs to microelectrodes using deposition by dielectrophoresis (DEP) (Chapter 4). We have systematically studied the use of surfactant molecules to disperse and horizontally align MWCNTs on substrates. In addition, DEP is shown to produce impurityfree placement of MWCNTs, forming connections between microelectrodes. We demonstrate the deposition density is tunable by both AC field strength and AC field frequency. 3. Etching of MWCNTs for the impurity-free nanoelectrodes (Chapter 5). We show that the residual Ni catalyst on MWCNTs can be removed by acid etching; the tip removal and collapsing of tubes into pyramids enhances the stability of field emission from the tube arrays. The acid-etching process can be used to functionalize the MWCNTs, which was used to make our initial CNT-nanoelectrode glucose sensors. Finally, lessons learned trying to perform spectroscopic analysis of the functionalized MWCNTs were vital for designing our final devices. 4. Molecular junction design and electrochemical synthesis of biphenyl molecules on carbon microelectrodes for all-carbon molecular devices (Chapter 6). Utilizing the experience gained on the work done so far, our final device design is described. We demonstrate the capability of preparing patterned glassy carbon films to serve as the bottom electrode in the new geometry. However, the molecular switching behavior of biphenyl was not observed by scanning tunneling microscopy (STM), mercury drop or fabricated glassy carbon/biphenyl/MWCNT junctions. Either the density of these molecules is not optimum for effective integration of devices using MWCNTs as the nanoelectrodes, or an electroactive contaminant was reduced instead of the ionic biphenyl species. 5. Self-assembly of octadecanethiol (ODT) molecules on gold microelectrodes for functional molecular devices (Chapter 7). We have realized an effective scheme to produce Au/ODT/MWCNT junctions by spanning MWCNTs across ODT-functionalized microelectrodes. A percentage of the resulting junctions retain the expected character of an ODT monolayer. While the process is not yet optimized, our successful junctions show that molecular electronic devices can be fabricated using simple processes such as photolithography, self-assembled monolayers and dielectrophoresis.

Graphite exfoliation by supercritical carbon dioxide extraction

Supercritical carbon dioxide is used to exfoliate graphite, producing a small, several-layer graphitic flake. The supercritical conditions of 2000, 2500, and 3000 psi and temperatures of 40°, 50°, and 60°C, have been used to study the effect of critical density on the sizes and zeta potentials of the treated flakes. Photon Correlation Spectroscopy (PCS), Brunauer-Emmett-Teller (BET) surface area measurement, field emission scanning electron microscopy (FE-SEM), and atomic force microscopy (AFM) are used to observe the features of the flakes. N-methyl-2-pyrrolidinone (NMP), dimethylformamide (DMF), and isopropanol are used as co-solvents to enhance the supercritical carbon dioxide treatment. As a result, the PCS results show that the flakes obtained from high critical density treatment (low temperature and high pressure) are more stable due to more negative charges of zeta potential, but have smaller sizes than those from low critical density (high temperature and low pressure). However, when an additional 1-hour sonication is applied, the size of the flakes from low critical density treatment becomes smaller than those from high critical density treatment. This is probably due to more CO2 molecules stacked between the layers of the graphitic flakes. The zeta potentials of the sonicated samples were slightly more negative than nonsonicated samples. NMP and DMF co-solvents maintain stability and prevented reaggregation of the flakes better than isopropanol. The flakes tend to be larger and more stable as the treatment time increases since larger flat area of graphite is exfoliated. In these experiments, the temperature has more impact on the flakes than pressure. The BET surface area resultsshow that CO2 penetrates the graphite layers more than N2. Moreover, the negative surface area of the treated graphite indicates that the CO2 molecules may be adsorbed between the graphite layers during supercritical treatment. The FE-SEM and AFM images show that the flakes have various shapes and sizes. The effects of surfactants can be observed on the FE-SEM images of the samples in one percent by weight solution of SDBS in water since the sodium dodecylbenzene sulfonate (SDBS) residue covers all of the remaining flakes. The AFM images show that the vertical thickness of the graphitic flakes can ranges from several nanometers (less than ten layers thick), to more than a hundred nanometers. In conclusion, supercritical carbon dioxide treatment is a promising step compared to mechanical and chemical exfoliation techniques in the large scale production of thin graphitic flake, breaking down the graphite flakes into flakes only a fewer graphene layers thick.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.