Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

Gradient-based reinforcement learning techniques for underwater robotics behavior learning

Darrerament, l'interès pel desenvolupament d'aplicacions amb robots submarins autònoms (AUV) ha crescut de forma considerable. Els AUVs són atractius gràcies al seu tamany i el fet que no necessiten un operador humà per pilotar-los. Tot i això, és impossible comparar, en termes d'eficiència i flexibilitat, l'habilitat d'un pilot humà amb les escasses capacitats operatives que ofereixen els AUVs actuals. L'utilització de AUVs per cobrir grans àrees implica resoldre problemes complexos, especialment si es desitja que el nostre robot reaccioni en temps real a canvis sobtats en les condicions de treball. Per aquestes raons, el desenvolupament de sistemes de control autònom amb l'objectiu de millorar aquestes capacitats ha esdevingut una prioritat. Aquesta tesi tracta sobre el problema de la presa de decisions utilizant AUVs. El treball presentat es centra en l'estudi, disseny i aplicació de comportaments per a AUVs utilitzant tècniques d'aprenentatge per reforç (RL). La contribució principal d'aquesta tesi consisteix en l'aplicació de diverses tècniques de RL per tal de millorar l'autonomia dels robots submarins, amb l'objectiu final de demostrar la viabilitat d'aquests algoritmes per aprendre tasques submarines autònomes en temps real. En RL, el robot intenta maximitzar un reforç escalar obtingut com a conseqüència de la seva interacció amb l'entorn. L'objectiu és trobar una política òptima que relaciona tots els estats possibles amb les accions a executar per a cada estat que maximitzen la suma de reforços totals. Així, aquesta tesi investiga principalment dues tipologies d'algoritmes basats en RL: mètodes basats en funcions de valor (VF) i mètodes basats en el gradient (PG). Els resultats experimentals finals mostren el robot submarí Ictineu en una tasca autònoma real de seguiment de cables submarins. Per portar-la a terme, s'ha dissenyat un algoritme anomenat mètode d'Actor i Crític (AC), fruit de la fusió de mètodes VF amb tècniques de PG.

Procés de destintatge del paper per flotació. Avaluació de l'eficàcia de l'eliminació de tinta

La tesi doctoral presentada té com a objectius principals l'estudi de les etapes fonamentals de desintegració i flotació en un procés de destintatge de papers vell de qualitats elevades per a poder millorar l'eficàcia d'aquestes etapes clau. Conté una revisió teòrica completa i molt actualitzada del procés de desintegració i flotació tant a nivell macroscòpic com microscòpic. La metodologia de treball en el laboratori, la posada a punt dels aparells, així com les anàlisis efectuades per a valorar la resposta del procés (anàlisi de blancor, anàlisi d'imatge i anàlisi de la concentració efectiva de tinta residual) estan descrites en el capítol de material i mètodes. La posada en marxa permet obtenir unes primeres conclusions respecte la necessitat de treballar amb una matèria primera homogènia i respecte la no significació de la temperatura de desintegració dins l'interval de treball permès al laboratori (20-50°C). L'anàlisi de les variables mecàniques de desintegració: consistència de desintegració (c), velocitat d'agitació en la desintegració (N) i temps de desintegració (t), permet de discernir que la consistència de desintegració és una variable fonamental. El valor de consistència igual al 10% marca el límit d'existència de les forces d'impacte mecànic en la suspensió fibrosa. A consistències superiors, les forces viscoses i d'acceleració dominen l'etapa de desintegració. Existeix una interacció entre la consistència i el temps de desintegració, optimitzant-se aquesta darrera variable en funció del valor de la consistència. La velocitat d'agitació és significativa només per a valors de consistència de desintegració inferiors al 10%. En aquests casos, incrementar el valor de N de 800 a 1400 rpm representa una disminució de 14 punts en el factor de destintabilitat. L'anàlisi de les variables químiques de desintegració: concentració de silicat sòdic (% Na2SiO3), peròxid d'hidrogen (% H2O2) i hidròxid sòdic (% Na2OH), proporciona resultats força significatius. El silicat sòdic presenta un efecte altament dispersant corroborat per les corbes de distribució dels diàmetres de partícula de tinta obtingudes mitjançant anàlisi d'imatges. L'hidròxid sòdic també presenta un efecte dispersant tot i que no és tant important com el del silicat sòdic. Aquests efectes dispersants són deguts principalment a l'increment de les repulsions electrostàtiques que aporten a la suspensió fibrosa aquests reactius químics fent disminuir l'eficàcia d'eliminació de l'etapa de flotació. El peròxid d'hidrogen utilitzat generalment com agent blanquejant, actua en aquests casos com a neutralitzador dels grups hidroxil provinents tant del silicat sòdic com de l'hidròxid sòdic, disminuint la repulsió electrostàtica dins la suspensió. Amb l'anàlisi de les variables hidrodinàmiques de flotació: consistència de flotació (c), velocitat d'agitació durant la flotació (N) i cabal d'aire aplicat (q), s'aconsegueix la seva optimització dins el rang de treball permès al laboratori. Valors elevats tant de la velocitat d'agitació com del cabal d'aire aplicat durant la flotació permeten eliminar majors quantitats de tinta. La consistència de flotació assoleix valors òptims depenent de les condicions de flux dins la cel·la de flotació. Les metodologies d'anàlisi emprades permeten obtenir diferents factors de destintabilitat. Entre aquests factors existeix una correlació important (determinada pels coeficients de correlació de Pearson) que permet assegurar la utilització de la blancor com a paràmetre fonamental en l'anàlisi del destintatge de papers vells, sempre i quan es complementi amb anàlisis d'imatge o bé amb anàlisi de la concentració efectiva de tinta residual. S'aconsegueixen expressions empíriques tipus exponencial que relacionen aquests factors de destintabilitat amb les variables experimentals. L' estudi de les cinètiques de flotació permet calcular les constants cinètiques (kBlancor, kERIC, kSupimp) en funció de les variables experimentals, obtenint un model empíric de flotació que relacionant-lo amb els paràmetres microscòpics que afecten realment l'eliminació de partícules de tinta, deriva en un model fonamental molt més difícil d'interpretar. Mitjançant l'estudi d'aquestes cinètiques separades per classes, també s'aconsegueix determinar que l'eficàcia d'eliminació de partícules de tinta és màxima si el seu diàmetre equivalent és superior a 50 μm. Les partícules amb diàmetres equivalents inferiors a 15 μm no són eliminades en les condicions de flotació analitzades. Es pot dir que és físicament impossible eliminar partícules de tinta de diàmetres molt diferents amb la mateixa eficiència i sota les mateixes condicions de treball. El rendiment del procés analitzat en funció de l'eliminació de sòlids per l'etapa de flotació no ha presentat relacions significatives amb cap de les variables experimentals analitzades. Únicament es pot concloure que addicionar quantitats elevades de silicat sòdic provoca una disminució tant de sòlids com de matèria inorgànica presents en les escumes de flotació.