A molecular dynamics coarse-grained block copolymer (BCP) model was developed and used to studied directed self-assembly (DSA), especially in regards to applications for semiconductor manufacturing. Most of the thesis is spent investigating the effect that guiding layer properties and block copolymer properties have on line roughness and defect density in a BCP-DSA process. These two effects are perhaps the most critical in making BCP-DSA a cost efficient industrial process. It is found that guiding patterns have little effect on line roughness and in fact that the BCP heals the majority of roughness in the underlying pattern. BCP properties have a larger effect on line roughness. Segregation strength (as measured by χN, where χ is the Flory- Huggins interaction parameter and N is the degree of polymerization) resulted in a larger than expected increase in line roughness when χN was low. Polydispersity resulted in a moderate increase in line roughness. In regards to equilibrium defect density, free energy calculations showed that χ was the primary determining factor, not χN as many expected. Equilibrium defect density was found to decrease exponentially with increasing χ. Defect density is also found to scale exponentially with polydispersity. Concerning defect heal rate, which can increase the real defect rate of a process if said rate is too low, it is found that increasing χN linearly increased the barrier to defect healing, which means that the defect heal rate decreases exponentially. However, for thin films this is only true for χN > ~ 50. Below χN ~ 50, the barrier is approximately constant. These results give excellent guidance to the type of materials and processes necessary to optimize a BCP-DSA process. A simulation technique designed to more efficiently sample over energy barriers called protracted noise dynamics for polymer systems was developed and studied. It was found that a decrease in simulation time of up to 4 orders of magnitude was achieved. The effect of box size on allowable pitches for a lamellar forming BCP was derived and demonstrated. It was found that more elongated boxes yielded more possible pitches and more accurate results. A short study on the effect of multiblock copolymers on the location of the order-disorder transition was also carried out and it was found that multiblock copolymers had small effect on the ODT. The distribution of chain conformations was also calculated.

The coarse-grained configurational-bias Monte Carlo (CG-CBMC) simulation was employed to study the structural behaviors of methyl methacrylate (MMA) oligomers adsorbed on grooved substrate due to molecular dynamics (MD) simulation is probably trapped at some local energy minima and difficult to carry out over a long enough time to allow relaxation of chain motion for an enormous polymeric system. Therefore, the CG-CBMC simulation was adopted in the present study. In this study, three types of chains are classified according to their positions relative to the groove. Type 1, Type 2, and Type 3 represent the entire MMA-oligomer within the groove, the MMA-oligomer partially within the groove, and the oligomer outside the groove, respectively. The orientational order parameters of Type 1 and Type 2 oligomers decrease with the increase of groove width, but the orientational order parameter of Type 3 oligomers is approximately equal to 0.1. In addition, observation of the orientational order parameters of Type 1 oligomers interacting with the grooved substrate at different interaction strengths decrease with increasing the groove width. Furthermore, the orientational order parameters of Type 1 oligomers within the narrowest (20 Å) and the widest (35 Å) groove with different depths were determined. For the narrowest groove, the arrangement of Type 1 oligomers will be influenced by the groove width. However, in the case of the widest groove, the orientational order parameter of Type 1 oligomers is approximately equal to 0.2. This study can help engineers clarify the characteristics and phenomena of physical adsorption of the molecules, as well as contributing to the application of recent technology.

No description available.

In the first part of the thesis, the molecular structure and electronic properties of phospholipids at the single molecule level and also for a monolayer structure are investigated via ab initio calculations under different degrees of hydration. The focus of the study is on phosphatidylcholines, in particular dipalmitoylphosphatidylcholine (DPPC), which are the most abundant phospholipids in biological membranes. Upon hydration, the phospholipid shape into a sickle-like structure. The hydration dramatically alters the surface potential, dipole and quadrupole moments of the lipids, and probably guides the interactions of the lipids with other molecules and the communication between cells. The vibrational spectrum of DPPC and DPPC-water complexes are completely assigned and it is shown that water hydrating the lipid head groups enables efficient energy transfer across membrane leaflets on sub-picosecond time scales. Moreover, the vibrational modes and lifetimes of pure and hydrated DPPC lipids, at human body temperature, are estimated by performing ab initio molecular dynamics simulations. The vibrational modes of the water molecules close to the head group of DPPC are active in the frequency range between 0.5 - 55 THz, with a peak at 2.80 THz in the energy spectrum. The computed lifetimes for the high-frequency modes agree well with recent data measured at room temperature, where high-order phonon scattering is not negligible. The structure and auto-ionization of water at the water-phospholipid interface are investigated by ab initio molecular dynamics and ab initio Monte Carlo simulations using local density approximation and generalized gradient approximation for the exchange-correlation energy functional. Depending on the lipid head group, strongly enhanced ionization is observed, leading to dissociation of several water molecules into H+ and OH- per lipid. The results can shed light on the phenomena of the high proton conductivity along membranes that has been reported experimentally. In the second part of the thesis, Monte Carlo simulations of the lipid bilayer, on the basis of a coarse grained model, are performed to gain insight into the mechanical properties of planar lipid bilayers. By using a rescaling method, the Poisson's ratio is calculated for different phases. Additional information on the bending rigidity, determined from height fluctuations on the basis of the Helfrich Hamiltonian, allows for calculation of the Young's modulus for each phase. In addition, the free energy barrier for lipid flip-flop process in the fluid and gel phases are estimated. The main rate-limiting step to complete a flip-flop process is related to a free energy barrier that has to be crossed in order to reach the center of the bilayer. The free energy cost for performing a lipid flip-flop in the gel phase is found to be five times greater than in the fluid phase, demonstrating the rarity of such events in the gel phase. Moreover, an energy barrier is estimated for formation of transient water pores that often precedes lipid translocation events and accounts for the rate-limiting step of these pore-associated lipid translocation processes.