As funções de correlação temporal para o momento de dipolo induzido do metano CH4 líquido, foram simuladas pelo método de dinâmica molecular em quatro estados termodinâmicos, a densidades de 28,10; 27,29; 26,41 e 25,34 mol/dm3 e temperaturas de 91,2; 100,7; 110,7 e 122,2 K respectivamente. Foram avaliadas as funções de 2- 3- e 4-corpos, mostraram sua importância no efeito de cancelamento e exibiram no tempo longo decaimentos aproximadamente exponenciais. Incluiu-se os mecanismos de indução emanados pelos campos elétricos dos momentos de octopolo e hexadecapolo, seus campos gradientes, assim como a superposição eletrônica. Con- seguiu-se reproduzir com concordância muito boa os espectros experimentais, no intervalo de freqüências de 0-600 cm-1. O dipolo induzido por octopolo representa a mais importante contribuição ao espectro de absorção no infravermelho longínquo, seguido do dipolo induzido por hexadecapolo. Os momentos espectrais zero estão perto dos valores experimentais, tendo o primeiro estado termodinâmico o melhor valor, havendo uma porcentagem média de erro dos quatro estados termodinâmicos igual ao 6 % a mais do valor experimental. / Time correlations functions for the induced dipole moment for liquid methane were simulated by molecular dynamics in four states of density 28,10, 27,29, 26,41, 25,34 mol/dm3 and temperature 91,2, 100,2, 110,7, 122,2 K respectively. The calculation of n-body terms show their importance on cancellation effects and exhibit in long time aproximately exponential decay. The induction mechanisms included were octupole and hexadecapole moments emanated by the electric fields and their field gradients, as well as electronic overlap. The simulated spectra described agrees very good with the experimental spectra in the frequencies of 0-600 cm−1, the octupole induced dipole represents the most important contribution to the far infrared absorption spectrum, followed by the hexadecapole induced dipole. The zero spectral moments were nearly the experimental values, having the first thermodynamic state the best value, resulting in an average percentual of errors in the four thermodynamics states of 6% over the experimental value.

As funções de correlação temporal para o momento de dipolo induzido do metano CH4 líquido, foram simuladas pelo método de dinâmica molecular em quatro estados termodinâmicos, a densidades de 28,10; 27,29; 26,41 e 25,34 mol/dm3 e temperaturas de 91,2; 100,7; 110,7 e 122,2 K respectivamente. Foram avaliadas as funções de 2- 3- e 4-corpos, mostraram sua importância no efeito de cancelamento e exibiram no tempo longo decaimentos aproximadamente exponenciais. Incluiu-se os mecanismos de indução emanados pelos campos elétricos dos momentos de octopolo e hexadecapolo, seus campos gradientes, assim como a superposição eletrônica. Con- seguiu-se reproduzir com concordância muito boa os espectros experimentais, no intervalo de freqüências de 0-600 cm-1. O dipolo induzido por octopolo representa a mais importante contribuição ao espectro de absorção no infravermelho longínquo, seguido do dipolo induzido por hexadecapolo. Os momentos espectrais zero estão perto dos valores experimentais, tendo o primeiro estado termodinâmico o melhor valor, havendo uma porcentagem média de erro dos quatro estados termodinâmicos igual ao 6 % a mais do valor experimental. / Time correlations functions for the induced dipole moment for liquid methane were simulated by molecular dynamics in four states of density 28,10, 27,29, 26,41, 25,34 mol/dm3 and temperature 91,2, 100,2, 110,7, 122,2 K respectively. The calculation of n-body terms show their importance on cancellation effects and exhibit in long time aproximately exponential decay. The induction mechanisms included were octupole and hexadecapole moments emanated by the electric fields and their field gradients, as well as electronic overlap. The simulated spectra described agrees very good with the experimental spectra in the frequencies of 0-600 cm−1, the octupole induced dipole represents the most important contribution to the far infrared absorption spectrum, followed by the hexadecapole induced dipole. The zero spectral moments were nearly the experimental values, having the first thermodynamic state the best value, resulting in an average percentual of errors in the four thermodynamics states of 6% over the experimental value.

As funções de correlação temporal para o momento de dipolo induzido do metano CH4 líquido, foram simuladas pelo método de dinâmica molecular em quatro estados termodinâmicos, a densidades de 28,10; 27,29; 26,41 e 25,34 mol/dm3 e temperaturas de 91,2; 100,7; 110,7 e 122,2 K respectivamente. Foram avaliadas as funções de 2- 3- e 4-corpos, mostraram sua importância no efeito de cancelamento e exibiram no tempo longo decaimentos aproximadamente exponenciais. Incluiu-se os mecanismos de indução emanados pelos campos elétricos dos momentos de octopolo e hexadecapolo, seus campos gradientes, assim como a superposição eletrônica. Con- seguiu-se reproduzir com concordância muito boa os espectros experimentais, no intervalo de freqüências de 0-600 cm-1. O dipolo induzido por octopolo representa a mais importante contribuição ao espectro de absorção no infravermelho longínquo, seguido do dipolo induzido por hexadecapolo. Os momentos espectrais zero estão perto dos valores experimentais, tendo o primeiro estado termodinâmico o melhor valor, havendo uma porcentagem média de erro dos quatro estados termodinâmicos igual ao 6 % a mais do valor experimental. / Time correlations functions for the induced dipole moment for liquid methane were simulated by molecular dynamics in four states of density 28,10, 27,29, 26,41, 25,34 mol/dm3 and temperature 91,2, 100,2, 110,7, 122,2 K respectively. The calculation of n-body terms show their importance on cancellation effects and exhibit in long time aproximately exponential decay. The induction mechanisms included were octupole and hexadecapole moments emanated by the electric fields and their field gradients, as well as electronic overlap. The simulated spectra described agrees very good with the experimental spectra in the frequencies of 0-600 cm−1, the octupole induced dipole represents the most important contribution to the far infrared absorption spectrum, followed by the hexadecapole induced dipole. The zero spectral moments were nearly the experimental values, having the first thermodynamic state the best value, resulting in an average percentual of errors in the four thermodynamics states of 6% over the experimental value.

Este trabalho investiga a termodifusão em coloides magnéticos através da técnica de varredura Z. O ponto de partida é a generalização do modelo de lente térmica, supondo o surgimento de um gradiente de concentração dos grãos magnéticos devido ao gradiente de temperatura causado pelo feixe de laser sobre a amostra. A partir do uso da técnica de varredura Z foi possível o estudo do coeficiente Soret (S IND.T) em ferrofluidos iônicos, surfactados e citrados, em amostras com baixa concentração de grãos (fração volumétrica de Fe, ø, menor que 1%). Na generalização do modelo de lente térmica que efetuamos, consideramos que a variação no índice de refração da amostra, em uma experiência de varredura Z, depende da variação da intensidade do feixe laser (I), da variação da temperatura (T) e da variação da concentração de grãos magnéticos (ø), onde C IND.N, C IND.T e C IND.S são seus respectivos parâmetros adimensionais no modelo. Uma vez que o tempo característico da termodifusão é da ordem de segundos, uma varredura Z com pulsos da ordem de 20ms é utilizada para a determinação de C IND.N. C IND.T é obtido independentemente por meio de métodos de óptica linear.Após a determinação de C IND.N e C IND.T, uma varredura Z com duração de pulso da ordem de 1 segundo é feita para determinar C IND.S e, posteriormente, o coeficiente Soret. A partir do comportamento da curva de evolução temporal da tranmitância com pulsos de 1 segundo pode-se determinar o sinal do coeficiente Soret. O sinal está relacionado com a tendência dos grãos de migrarem para a região mais fria (termofóbico, S IND.T>0) ou mais quente (termofílico, S IND.T<0) da amostra, dependendo de suas características físico-químicas. Mostramos que o módulo de S IND.T é proporcional a ø, em concordância com resultados obtidos para soluções mais concentradas (ø1) através da técnica de Espalhamento Rayleigh Forçado. Uma possível origem física para os comportamentos termofóbico e termofílico dos fluidos magnéticos poderia estar relacionada a mudanças na intensidade das forças que mantêm o equilíbrio coloidal, por ação da temperatura. / This work explores the thermodiffusion in magnetic colloids through the Z-Scan technique. The starting point is the generalization of the thermal lens model based on the assumption that the concentration gradient of the magnetic grains emerges due to the temperature gradient caused by the laser beam on the sample. By using the Z-Scan technique it was possible to study the Soret coefficient (ST) for ionic, surfacted and citrated ferrofluids in samples with low concentration of grains (Fe volumetric percentage, ø, less than 1%). In this thermal lens model generalization, we have considered that the refraction índex variation in a Z-Scan experiment depends on the laser beam intensity (I), the temperature variation (T) and the variation of the magnetic grains concentration (ø), where CN, CT and CS are their respective dimensionless parameters in the model. As characteristic time of thermodiffusion is of the order of seconds, a Z-Scan with pulses around 20 ms is used in order to determine CN. CT is obtained independently by using lenear optics methods. After the determination of CN and CT, a Z-Scan with pulses around 1 second is made in order to determine CS and, Consequently, the Soret coefficient. Through the behavior of the time dependent transmittance with 1-second pulses we were able to determine the sign of the Soret coefficient. The sign is related to the tendency of the grains to migrate to the colder region (thermophobic, ST>0) or to the warmer region (thermophilic, ST<0) of the sample, depending on its physical-chemical characteristics. We have showed that the ST module is proportional to ø, in agreement with the results for higher concentration solutions (ø1%) obtained through Forced Rayleigh Scattering. A possible physical originfor the thermophobic and thermophilic behavior of magnetic fluids could berelated to changes in the intensity of the forces that keep the colloidal balance, by means of temperature.

Fizemos um estudo teórico sistemático e criterioso para entender o papael das moléculas de água na estabilização da geometria de alfa-ciclodextrina (alfa-CD) em solução aquosa. Adicionalmente, analisamos a confiabilidade dos campos de força clássicos amplamente utilizados e a aplicabilidade de geometrias otimizadas com cálculos quânticos, através de comparação com resultados experimentais. Inicialmente realizamos otimizações de geometria da alfa-CD isolada utilizando os métodos clássicos AMBER, OPLS e MM2 e os métodos quânticos Hartree-Fock (HF), AM1 e PM3. Dentre as 42 geometrias encontradas, identificamos a mais estável obtida com métodos clássicos (AMBER) e com método quântico (HF). Através de comparação de propriedades estruturais dessas duas geometrias com os dados experimentais, observamos que elas são muito próximas e que diferem basicamente pelo tamanho da cavidade, abertura do cone e conformação das hidroxilas primárias. Na geometria AMBER essas hidroxilas são todas (-)gauche enquanto que na HF há duas hidroxilas na conformação trans provocadas pela formação de ligações de hidrogênio intramoleculares e com isso dando uma estabilidade maior para essa geometria de cerca de 2 kcal/mol. Utilizando simulações computacionais com método de Monte Carlo, analisamos a solvatação dessas duas geometrias em água. Verificamos que a geometria AMBER, embora tendo um dipolo menor, interage mais fortemente com a água por cerca de 21 kcal/mol. Isso é decorrência de 6 ligações de hidrogênio a mais que essa geometria faz com a água em comparação com a geometria HF. Analisamos também o caráter hidrofóbico da cavidade e calculamos a densidade interna na geometria HF e 0.33 g/cm POT.3 enquanto que na cavidade da geometria AMBER essa densidade é 0.46 g/cm POT.3. O efeito hidrofóbico provocado pela cavidade em ambas as geometrias é bastante forte, sendo na geometria HF um pouco maior. Para finalizar nossa análise, ) realizamos simulações com Teoria de Perturbação Termodinâmica (PT) e calculamos a variação de energia livre entre essas duas geometrias em solução aquosa. Encontramos que a geometria AMBER é mais estável em água por cerca de 20 kcal/mol. Como a entropia é dada por ( alfaF-alfaE)/T, temos que a variação da entropia entre os dois estados é desprezível. Assim verificamos que dentre os métodos teóricos utilizados para otimização de geometria, o campo de força clássico AMBER foi o que melhor descreveu a geometria da alfa-CD / Theoretical study of the structural stability and hydrophilic solvation of \'alpha\'-cyclodextrin

Fizemos um estudo teórico sistemático e criterioso para entender o papael das moléculas de água na estabilização da geometria de alfa-ciclodextrina (alfa-CD) em solução aquosa. Adicionalmente, analisamos a confiabilidade dos campos de força clássicos amplamente utilizados e a aplicabilidade de geometrias otimizadas com cálculos quânticos, através de comparação com resultados experimentais. Inicialmente realizamos otimizações de geometria da alfa-CD isolada utilizando os métodos clássicos AMBER, OPLS e MM2 e os métodos quânticos Hartree-Fock (HF), AM1 e PM3. Dentre as 42 geometrias encontradas, identificamos a mais estável obtida com métodos clássicos (AMBER) e com método quântico (HF). Através de comparação de propriedades estruturais dessas duas geometrias com os dados experimentais, observamos que elas são muito próximas e que diferem basicamente pelo tamanho da cavidade, abertura do cone e conformação das hidroxilas primárias. Na geometria AMBER essas hidroxilas são todas (-)gauche enquanto que na HF há duas hidroxilas na conformação trans provocadas pela formação de ligações de hidrogênio intramoleculares e com isso dando uma estabilidade maior para essa geometria de cerca de 2 kcal/mol. Utilizando simulações computacionais com método de Monte Carlo, analisamos a solvatação dessas duas geometrias em água. Verificamos que a geometria AMBER, embora tendo um dipolo menor, interage mais fortemente com a água por cerca de 21 kcal/mol. Isso é decorrência de 6 ligações de hidrogênio a mais que essa geometria faz com a água em comparação com a geometria HF. Analisamos também o caráter hidrofóbico da cavidade e calculamos a densidade interna na geometria HF e 0.33 g/cm POT.3 enquanto que na cavidade da geometria AMBER essa densidade é 0.46 g/cm POT.3. O efeito hidrofóbico provocado pela cavidade em ambas as geometrias é bastante forte, sendo na geometria HF um pouco maior. Para finalizar nossa análise, ) realizamos simulações com Teoria de Perturbação Termodinâmica (PT) e calculamos a variação de energia livre entre essas duas geometrias em solução aquosa. Encontramos que a geometria AMBER é mais estável em água por cerca de 20 kcal/mol. Como a entropia é dada por ( alfaF-alfaE)/T, temos que a variação da entropia entre os dois estados é desprezível. Assim verificamos que dentre os métodos teóricos utilizados para otimização de geometria, o campo de força clássico AMBER foi o que melhor descreveu a geometria da alfa-CD / Theoretical study of the structural stability and hydrophilic solvation of \'alpha\'-cyclodextrin

Neste trabalho apresentamos uma revisão de métodos de simulação de energia eletrostática de sistemas de cargas e uma proposta de adaptação de algoritmo ultilizado na literatura de sistemas gravitacionais para estudo das propriedades estatísticas de sistemas coulombianos. Na primeira parte do estudo, revisamos os fundamentos teóricos do método de Ewald e suas condições de aplicabilidade, procurando esclarecer as referências mais importantes no assunto, que são de difícil compreensão, gerando equívocos na utilização do termo de dipolo. Detalhamos o estudo sobre a análise da convergência da série em que a técnica se baseia, bem como sua interpretação física mostrando a equivalência entre as duas abordagens . Na segunda parte do trabalho analisamos os fundamentos do Fast Multipole Method desenvolvido para interação gravitacional, para o qual construímos programas em linguagem C para uma versão na rede. Criamos um algoritmo que denominamos Fast Multipole Monte Carlo (FMMC) e desenvolvemos um programa para cálculo das propriedades termodinâmicas de sistemas coulombianos. Os programas são testados comparando resultados para a energia e propriedades térmicas do modelo LRPM com resultados de simulação através de cálculo direto. / In this work we present a review of methods of simulation for the electrostatic energy of charged systems and an adaptation of an algorithm from the literature on gravitational systems for the study of the statistical properties of Coulomb systems. In the first part of the work, we review the fundamentals for the theoretical method of Ewald and its conditions of applicability, seeking to clarify the most important references on the subject, which because of the involved mathematics, have led to misuse of the so-called dipole correction. We detail the study on the convergence of the series for the electrostatic potential on which the Ewald technique is based, as well as the physical interpretation given elsewhere, showing the equivalence between the two approaches. In the second part of this work, we analyse the foundations of the Fast Multipole Method developed for gravitational interactions, and present programs in C language for a network version of neutral charged systems. Finally, an algorithm, which we name Fast Multipole Monte Carlo, and the corresponding code for calculating the thermodynamic properties of Coulomb systems are presented. The programs are tested by comparing results for the energy and thermal properties of the Lattice Restricted Primitive model with results of simulations based on direct calculations for the Coulomb energies.

Este trabalho investiga a termodifusão em coloides magnéticos através da técnica de varredura Z. O ponto de partida é a generalização do modelo de lente térmica, supondo o surgimento de um gradiente de concentração dos grãos magnéticos devido ao gradiente de temperatura causado pelo feixe de laser sobre a amostra. A partir do uso da técnica de varredura Z foi possível o estudo do coeficiente Soret (S IND.T) em ferrofluidos iônicos, surfactados e citrados, em amostras com baixa concentração de grãos (fração volumétrica de Fe, ø, menor que 1%). Na generalização do modelo de lente térmica que efetuamos, consideramos que a variação no índice de refração da amostra, em uma experiência de varredura Z, depende da variação da intensidade do feixe laser (I), da variação da temperatura (T) e da variação da concentração de grãos magnéticos (ø), onde C IND.N, C IND.T e C IND.S são seus respectivos parâmetros adimensionais no modelo. Uma vez que o tempo característico da termodifusão é da ordem de segundos, uma varredura Z com pulsos da ordem de 20ms é utilizada para a determinação de C IND.N. C IND.T é obtido independentemente por meio de métodos de óptica linear.Após a determinação de C IND.N e C IND.T, uma varredura Z com duração de pulso da ordem de 1 segundo é feita para determinar C IND.S e, posteriormente, o coeficiente Soret. A partir do comportamento da curva de evolução temporal da tranmitância com pulsos de 1 segundo pode-se determinar o sinal do coeficiente Soret. O sinal está relacionado com a tendência dos grãos de migrarem para a região mais fria (termofóbico, S IND.T>0) ou mais quente (termofílico, S IND.T<0) da amostra, dependendo de suas características físico-químicas. Mostramos que o módulo de S IND.T é proporcional a ø, em concordância com resultados obtidos para soluções mais concentradas (ø1) através da técnica de Espalhamento Rayleigh Forçado. Uma possível origem física para os comportamentos termofóbico e termofílico dos fluidos magnéticos poderia estar relacionada a mudanças na intensidade das forças que mantêm o equilíbrio coloidal, por ação da temperatura. / This work explores the thermodiffusion in magnetic colloids through the Z-Scan technique. The starting point is the generalization of the thermal lens model based on the assumption that the concentration gradient of the magnetic grains emerges due to the temperature gradient caused by the laser beam on the sample. By using the Z-Scan technique it was possible to study the Soret coefficient (ST) for ionic, surfacted and citrated ferrofluids in samples with low concentration of grains (Fe volumetric percentage, ø, less than 1%). In this thermal lens model generalization, we have considered that the refraction índex variation in a Z-Scan experiment depends on the laser beam intensity (I), the temperature variation (T) and the variation of the magnetic grains concentration (ø), where CN, CT and CS are their respective dimensionless parameters in the model. As characteristic time of thermodiffusion is of the order of seconds, a Z-Scan with pulses around 20 ms is used in order to determine CN. CT is obtained independently by using lenear optics methods. After the determination of CN and CT, a Z-Scan with pulses around 1 second is made in order to determine CS and, Consequently, the Soret coefficient. Through the behavior of the time dependent transmittance with 1-second pulses we were able to determine the sign of the Soret coefficient. The sign is related to the tendency of the grains to migrate to the colder region (thermophobic, ST>0) or to the warmer region (thermophilic, ST<0) of the sample, depending on its physical-chemical characteristics. We have showed that the ST module is proportional to ø, in agreement with the results for higher concentration solutions (ø1%) obtained through Forced Rayleigh Scattering. A possible physical originfor the thermophobic and thermophilic behavior of magnetic fluids could berelated to changes in the intensity of the forces that keep the colloidal balance, by means of temperature.

Neste trabalho apresentamos uma revisão de métodos de simulação de energia eletrostática de sistemas de cargas e uma proposta de adaptação de algoritmo ultilizado na literatura de sistemas gravitacionais para estudo das propriedades estatísticas de sistemas coulombianos. Na primeira parte do estudo, revisamos os fundamentos teóricos do método de Ewald e suas condições de aplicabilidade, procurando esclarecer as referências mais importantes no assunto, que são de difícil compreensão, gerando equívocos na utilização do termo de dipolo. Detalhamos o estudo sobre a análise da convergência da série em que a técnica se baseia, bem como sua interpretação física mostrando a equivalência entre as duas abordagens . Na segunda parte do trabalho analisamos os fundamentos do Fast Multipole Method desenvolvido para interação gravitacional, para o qual construímos programas em linguagem C para uma versão na rede. Criamos um algoritmo que denominamos Fast Multipole Monte Carlo (FMMC) e desenvolvemos um programa para cálculo das propriedades termodinâmicas de sistemas coulombianos. Os programas são testados comparando resultados para a energia e propriedades térmicas do modelo LRPM com resultados de simulação através de cálculo direto. / In this work we present a review of methods of simulation for the electrostatic energy of charged systems and an adaptation of an algorithm from the literature on gravitational systems for the study of the statistical properties of Coulomb systems. In the first part of the work, we review the fundamentals for the theoretical method of Ewald and its conditions of applicability, seeking to clarify the most important references on the subject, which because of the involved mathematics, have led to misuse of the so-called dipole correction. We detail the study on the convergence of the series for the electrostatic potential on which the Ewald technique is based, as well as the physical interpretation given elsewhere, showing the equivalence between the two approaches. In the second part of this work, we analyse the foundations of the Fast Multipole Method developed for gravitational interactions, and present programs in C language for a network version of neutral charged systems. Finally, an algorithm, which we name Fast Multipole Monte Carlo, and the corresponding code for calculating the thermodynamic properties of Coulomb systems are presented. The programs are tested by comparing results for the energy and thermal properties of the Lattice Restricted Primitive model with results of simulations based on direct calculations for the Coulomb energies.

As mudanças nas propriedades estruturais, eletrônicas, óticas e magnéticas de moléculas em aglomerados e em meio líquido são apresentadas. Para os aglomerados moleculares tais propriedades foram obtidas usando métodos ab initio tradicionais de mecânica quântica em configurações otimizadas de mínima energia. A fase líquida é mais complexa e ainda hoje se mostra como um desafio teórico. Modelos convencionais para a descrição teórica dos efeitos de solvente, como o campo de reação auto-consistente e a aproximação de aglomerado rígido, apesar de largamente aplicados não satisfazem importantes características da fase líquida como, por exemplo, a natureza estatística dos líquidos. Adicionalmente, interações específicas tais como, formação de ligações de hidrogênio e transferência de carga são mais difíceis de incluir. Modelos mais realísticos para tratar a estrutura eletrônica de sistemas líquidos envolvem algum tipo de simulação computacional. Entre estes, os métodos conhecidos como QM/MM são modelos híbridos aplicados com sucesso e que empregam tanto simulações clássicas como métodos de mecânica quântica. Existem hoje dois principais tipos de cálculos QM/MM. Nos cálculos convencionais, o sistema é particionado em duas regiões (clássica e quântica) e0 as interações são calculadas separada e simultaneamente. Nos cálculos seqüenciais QM/MM, as configurações do líquido são geradas via simulações clássicas e posteriormente submetidas a cálculos quânticos. Neste trabalho a metodologia seqüencial (Simulação/MQ) foi empregada no estudo de sistemas em fase líquida, como soluções aquosas de piridina e metanol e líquidos homogêneos como benzeno e água. Empregando simulações Monte Carlo ou Dinâmica Molecular as estruturas do líquido são geradas à temperatura ambiente para subsequentes cálculos quânticos. Os resultados para o líquido são comparados com aqueles obtidos para a molécula isolada e fornecem uma clara descrição dos efeitos de solvente. Entre as propriedades estudadas nesta tese estão energia de interação, momento de dipolo, freqüência vibracional, polarizabilidade, blindagem magnética e propriedades de espalhamento Raman e Rayleigh. / The changes of structural, electronic, optical and magnetic properties of molecules in clusters and in liquid phase are presented. For molecular clusters such properties were obtained using traditional and well defined quantum mechanics of initio methods in the minimum energy geometry - optimized configurations. The liquid phase is more complex and represents a current theoretical challenge. Conventional models for theoretical description of solvent effects, such as self-consistent reaction field and rigid cluster approximation, although widely applied do not satisfy important characteristics of the liquid phase such as the formation of hydrogen bonds and charge transfer are more difficult to include. More realistic models to treat the electronic structure of liquid systems include some sort of computer simulation. Among these the so-called QM/MM methods are successful hybrid models that uses both the classical simulation and quantum mechanics. There are now two main classes of QM/MM calculations. In the conventional QM/MM the system is partitioned into two regions (classical and quantum) and the interactions are calculated separated and simultaneously. In the sequential QM/MM, the configurations of the liquid are generated via classical simulations and subsequently the quantum calculations are performed. In this work the sequential methodology (Simulation/MQ) is used to study molecular systems in the liquid phase as aqueous solutions of pyridine and methanol and homogeneous liquid as water and benzene. Using Monte Carlo or Molecular Dynamics simulations liquid structures at room temperature are generated for subsequent quantum mechanics calculations. The results for the liquid are compared with those obtained for the isolated molecule and give a clear picture of the liquid effects. Among the properties studied in this thesis are interaction energies, dipole moments, vibrational frequencies, polarizabilities, magnetic shielding, chemical shift and optical properties of the Raman and Rayleigh scattering.