← Back to branch summary

~ubuntu-branches/ubuntu/trusty/nwchem/trusty-proposed

« back to all changes in this revision

Viewing changes to src/nwdft/rt_tddft/rt_tddft.F

  • Committer: Package Import Robot
  • Author(s): Michael Banck, Daniel Leidert, Andreas Tille, Michael Banck
  • Date: 2013-07-04 12:14:55 UTC
  • mfrom: (1.1.2)
  • Revision ID: package-import@ubuntu.com-20130704121455-5tvsx2qabor3nrui
Tags: 6.3-1
http://bugs.debian.org/696361
* New upstream release.
* Fixes anisotropic properties (Closes: #696361).
* New features include:
  + Multi-reference coupled cluster (MRCC) approaches
  + Hybrid DFT calculations with short-range HF 
  + New density-functionals including Minnesota (M08, M11) and HSE hybrid
    functionals
  + X-ray absorption spectroscopy (XAS) with TDDFT
  + Analytical gradients for the COSMO solvation model
  + Transition densities from TDDFT 
  + DFT+U and Electron-Transfer (ET) methods for plane wave calculations
  + Exploitation of space group symmetry in plane wave geometry optimizations
  + Local density of states (LDOS) collective variable added to Metadynamics
  + Various new XC functionals added for plane wave calculations, including
    hybrid and range-corrected ones
  + Electric field gradients with relativistic corrections 
  + Nudged Elastic Band optimization method
  + Updated basis sets and ECPs 

[ Daniel Leidert ]
* debian/watch: Fixed.

[ Andreas Tille ]
* debian/upstream: References

[ Michael Banck ]
* debian/upstream (Name): New field.
* debian/patches/02_makefile_flags.patch: Refreshed.
* debian/patches/06_statfs_kfreebsd.patch: Likewise.
* debian/patches/07_ga_target_force_linux.patch: Likewise.
* debian/patches/05_avoid_inline_assembler.patch: Removed, no longer needed.
* debian/patches/09_backported_6.1.1_fixes.patch: Likewise.
* debian/control (Build-Depends): Added gfortran-4.7 and gcc-4.7.
* debian/patches/10_force_gcc-4.7.patch: New patch, explicitly sets
  gfortran-4.7 and gcc-4.7, fixes test suite hang with gcc-4.8 (Closes:
  #701328, #713262).
* debian/testsuite: Added tests for COSMO analytical gradients and MRCC.
* debian/rules (MRCC_METHODS): New variable, required to enable MRCC methods.

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/nwdft/rt_tddft/rt_tddft.F
 
1
C> @file rt_tddft.F
 
2
C> @ingroup rt_tddft
 
3
 
 
4
C--------------------------------------------------------------------
 
5
C
 
6
C  Overall RT-TDDFT documentation:
 
7
C
 
8
C> @defgroup rt_tddft Real-time TDDFT
 
9
C--------------------------------------------------------------------
 
10
 
 
11
C--------------------------------------------------------------------
 
12
C> Driver routine for real-time TDDFT.
 
13
C>       
 
14
C> Call this from main DFT routine after convergence.
 
15
C> It will call the closed-shell, open-shell, spin-orbit, etc versions of RT-TDDFT.
 
16
C--------------------------------------------------------------------
 
17
      logical function rt_tddft (rtdb)
 
18
      implicit none
 
19
 
 
20
#include "errquit.fh"
 
21
#include "mafdecls.fh"
 
22
#include "stdio.fh"
 
23
#include "global.fh"
 
24
#include "msgids.fh"
 
25
#include "util.fh"
 
26
#include "cdft.fh"
 
27
#include "matutils.fh"
 
28
#include "case.fh"
 
29
#include "rtdb.fh"
 
30
#include "headers/rt_tddft.fh"
 
31
 
 
32
C ^^^ strange that I had to do this
 
33
 
 
34
 
 
35
C     == Inputs ==
 
36
      integer, intent(in) :: rtdb !< Runtime database
 
37
 
 
38
      
 
39
C     == Parameters ==
 
40
      character(*), parameter :: pname = "rt_tddft: "
 
41
 
 
42
 
 
43
C     == Variables ==
 
44
      type(rt_params_t) :: params
 
45
      character*80 theory
 
46
      integer me
 
47
 
 
48
 
 
49
      rt_tddft = .false.
 
50
 
 
51
      me = ga_nodeid ()
 
52
 
 
53
 
 
54
 
 
55
C
 
56
C     Store rtdb in params struct, then determine the theory we are using.
 
57
C
 
58
      params%rtdb = rtdb
 
59
 
 
60
      if (.not. rtdb_cget(rtdb, "dft:theory", 1, theory))
 
61
     $     call errquit(pname//"theory not specified",0, RTDB_ERR)
 
62
 
 
63
      if (theory .eq. "dft") then
 
64
         
 
65
         if (ipol.eq.1) then             !closed shell DFT
 
66
            if (me.eq.0) then
 
67
               write (luout, *) " " 
 
68
               write (luout, *) "*************************"
 
69
               write (luout, *) "* Closed-shell RT-TDDFT *"
 
70
               write (luout, *) "*************************"
 
71
               write (luout, *) " "
 
72
               call util_flush (luout)
 
73
            endif
 
74
            call rt_tddft_cs_init (params)
 
75
            call rt_tddft_header (params)
 
76
            call rt_tddft_cs_prop (params)
 
77
            
 
78
         elseif (ipol.eq.2) then         !open shell DFT
 
79
            if (me.eq.0) then
 
80
               write (luout, *) " "
 
81
               write (luout, *) "***********************"
 
82
               write (luout, *) "* Open-shell RT-TDDFT *"
 
83
               write (luout, *) "***********************"
 
84
               write (luout, *) " "
 
85
               call util_flush (luout)
 
86
            endif
 
87
            call rt_tddft_os_init (params)
 
88
            call rt_tddft_header (params)
 
89
            call rt_tddft_os_prop (params)
 
90
         else
 
91
            call errquit (pname//"bad ipol value", 0, 0)
 
92
         endif
 
93
         
 
94
      elseif (theory .eq. "sodft") then  !spin orbit DFT
 
95
         call errquit (pname//"spin-orbit not implemented yet",0,0)
 
96
      else
 
97
         call errquit (pname//"Invalid theory: "//trim(theory), 0, 0)
 
98
      endif
 
99
 
 
100
 
 
101
      
 
102
C
 
103
C     Clean up.
 
104
C
 
105
      call rt_tddft_clean (params)
 
106
 
 
107
      
 
108
 
 
109
C
 
110
C     Finally post-process dens mats to cube files.  Note we explicitly pass subgs option...
 
111
C
 
112
      if (params%viz_dplot) call rt_tddft_snapshot_dplot (params,
 
113
     $     params%viz_subgs)
 
114
 
 
115
 
 
116
      rt_tddft = .true.
 
117
      end function